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^ . Abstract 

. 

This article reviews the current status of precursor superconducting phase fluctu- 
ations as a possible mechanism for pseudogap formation in high-temperature su- 
'"^5 per conductors. In particular we compare this approach which relies on the two- 
q ■ dimensional nature of the superconductivity to the often used T-matrix approach. 
O Starting from simple pairing Hamiltonians we present a broad pedagogical intro- 
duction to the BCS-Bose crossover problem. The finite temperature extension of 
these models naturally leads to a discussion of the Berezinskii-Kosterlitz-Thouless 
superconducting transition and the related phase diagram including the effects of 
OO quantum phase fluctuations and impurities. We stress the differences between sim- 
ple Bose-BCS crossover theories and the current approach where one can have a 
large pseudogap region even at high carrier density where the Fermi surface is well- 
<^> \ defined. The Green's function and its associated spectral function, which explicitly 
show non-Fermi liquid behaviour, is constructed in the presence of vortices. Finally 
"£3 ■ different mechanisms including quasi-particle- vortex and vortex- vortex interactions 
£h ! for the filling of the gap above T c are considered. 
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1 Introduction 



1.1 Pseudogap Phenomena in High-T c superconductors 



The discovery of high temperature superconductors (HTSC) [1] revealed new problems in solid state 
physics in general and in the theory of superconductivity in particular. A combination of factors, 
including unusual magnetic and electronic properties, a lowered dimensionality, closeness to the metal- 
insulator transition and relatively low carrier densities, makes the construction of an appropriate theory 
both difficult and far from resolved. There is no consensus as to the correct theoretical approach. It is 
even unclear as to which physical features of cuprate superconductors should be considered as the basis 
for the correct theory. 

One of the clearest differences between the BCS scenario of superconductivity and superconductivity 
in the cuprates is the existence of a pseudogap. This is simply a depletion of the single particle spectral 
weight around the Fermi level [2,3] (see also a recent excellent experimental review by Timusk and 
Statt [4]).0 The earliest experiments to reveal gap-like behaviour in the normal state were the NMR 
measurements of the Knight shift (see Refs. in [4]) which probes the uniform spin susceptibility. Based 
on the NMR data the phase diagram presented in Fig. 1 was suggested for the cuprate superconductors. 

3 Note that the word "pseudogap" for HTSC was originally suggested by Friedel [5], although in a different 
context. The notion of pseudogap has been used before in other areas such as ID conductors. 
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Fig. 1. Qualitative phase diagram of cuprates given in terms of the variables, temperature T and the delocalized 
carrier concentration x, taken from [4]. Pseudogap crossover temperatures were suggested from the NMR data. 
The crossover temperatures merge into T c slightly into the overdoped region of the phase diagram. 

In Fig. 1 there are two main phase boundaries. The first is the transition into a long-range antiferro- 
magnetic state at the Neel ordering temperature, T N , in the very low doping x regime. The second is 
the transition into a superconducting state, with the maximal critical temperature T c at approximately 
x = 0.2 holes/Cu0 2 unit cell. With respect to the maximal T c , under-, optimally and overdoped regimes 
are distinguished. 

In addition to the above-mentioned phase boundaries, two crossover lines, T* and T°, are evident in 
these materials, even at rather low doping levels. It must be noted, however, that the experimental 
determination of these lines is still somewhat controversial. We will return to the discussion of how the 
different theoretical approaches treat the phase diagram. 

At the upper crossover temperature T° the Knight shift changes behaviour. Above T° it is temperature 
independent, while below T° it decreases linearly with temperature. Below the lower crossover tem- 
perature T* the Knight shift decreases supra-linearly with temperature. It was suggested (see [4]) that 
these two crossover temperatures may originate from two different physical phenomena. In particular 
the lower crossover temperature T* appears to be the result of a gap appearing in the spectrum of 
elementary antiferromagnetic excitations. This manifestation of the pseudogap phenomenon was thus 
called a spin gap. 

Subsequently the optical conductivity data (which has been discussed in detail in the review [6]) showed 
in addition a gapping of the charge degrees of freedom. Specific heat data (see Refs. in [4]) also provided 
evidence that an electronic gap opens below T*. Considerable recent experimental progress in angle 
resolved photo-emission spectroscopy (ARPES) (see for example the reviews [7,8]) gives an experimental 
window on the single-particle spectral function and other fundamental quantities such as the electron 
self-energy [9] . In these experiments one can clearly see the presence of a strong reduction in the quasi- 



particle weight, or equivalently a pseudogap, below T*. Until recently most ARPES measurements were 
performed for the photon energy range between 19 and 25eV. The most recent data [10], which was 
taken with a higher photon energy of 33eV, is in contradiction with the previous data. This question 
is now under intensive investigation [11], particularly since ohmic losses may alter the ARPES spectra 
[12]. 

An explanation of the pseudogap phenomenon is regarded as one of the most important unresolved 
questions in the theory of superconductivity. The ARPES and the scanning tunnelling spectroscopy 
(STS) experiments (see Refs. in [4]) showed a smooth crossover from the pseudo- to superconducting 
gap, so that in these single-particle spectroscopies, the transition into the true superconducting state 
at T c is barely noticeable. Moreover the pseudo- and superconducting gaps also reveal experimentally 
the same ci-wave symmetry of the order parameter [13,14] (see also the reviews [15,4,8]). We note, 
however, that at this stage the question of the symmetry of the order parameter in HTSC is still under 
discussion [16]. For example, it is shown in [17-19] that the observed symmetry of the order parameter 
cannot be fitted by only the lowest harmonic of the rf-wave order parameter. Furthermore, a recent 
experiment with twisted Josephson junctions in the Bi-cuprates [20], is in favour of an extended s-wave 
order parameter, and has shown the absence of a d-wave part in the order parameter. It has also been 
suggested that a transition from the d x 2_ y 2 symmetry to the d x 2_ y 2 + id xy symmetry could be induced 
by an external magnetic field [21] and may be observed in the cuprates. Despite remaining uncertainties 
in the experimental data, the existence of a pseudogap and a corresponding crossover temperature is 
well established, and the underlying physics demands an explanation. 
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1.2 Current theoretical explanations for the origin of the pseudogap: magnetic, spin-density waves, 
stripes, etc. 

Due to the extremely complex structural and electronic nature of the cuprate systems [22], and the 
somewhat controversial nature of much of the current experimental data, there are many theoretical 



attempts to explain pseudogap behaviour. As an example, we show two theoretical phase diagrams [23]. 
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Fig. 2. Theoretical patterns for the phase diagrams, taken from [23]. 

The first is based on the theory developed by Pines and collaborators of the nearly antiferromagnetic 
Fermi liquid [24] (Fig. 2-A). The second is based on the arguments of Emery and Kivelson (Fig. 2- 
B) which are based on superconducting phase fluctuations [25,26] (see also [27]). Both these theories 
are attempting to interpret the experimental phase diagram presented in Fig. 1, and give in addition, 
predictions about the less well experimentally studied regions of the diagram. As an example, the 
continuation of the pseudogap lines below T c on Fig. 2 should be interpreted [23] as "what would be 
present if there was no superconducting phase" . This can be realized in practice by applying a magnetic 
field. 

At this stage most interpretations of the NMR experiments [28,29] support the magnetic theory of the 
pseudogap (or to be more exact the existence of a spin gap). On the other hand the ARPES [4,8,30] 
and STS [4,31,32] experiments can and have been interpreted using theories based solely on precursor 
superconducting fluctuations (see Sec. 1.3). This claim however depends on the precise way in which 
the experimental results have been interpreted, and how well developed the particular theory is for the 
given range of doping and/or the applied magnetic field if present. For instance, it is argued in [23] 
that the current NMR data still do not exclude the scenario based on the superconducting fluctuations 
and that more measurements in the overdoped cuprates are needed. An opposite point of view which 
supports the diagram in Fig. 2-A is presented in a very recent survey of the experimental data [33]. 
However, so evident disagreement between [33] and [23] (see also [3]) once more demonstrates vividly 



that there is no consensus about the true shape of the cuprate phase diagram. 

Another possible explanation relates the pseudogap to charge- and/or spin-density waves [16] (see 
also very recent paper [34] and Refs. therein). A theory based on the phase separation in the form 
of a strip structure of metallic and insulating domains with a spin gap in the AFM regions has been 
suggested by Emery, Kivelson and Zachar [35] (see also the more recent paper [36]). Such a separation is 
accompanied by the formation of a pairing amplitude through a proximity effect between non-magnetic 
metallic and AFM insulating stripes, but with large phase fluctuations. Below T* this theory has much 
in common with the more phenomenological approach based on the 2D attractive Hubbard model 
discussed here simply because both theories are based on phase fluctuations. A stripe-phase quantum- 
critical-point scenario for high-T c superconductors which has much in common with that of Ref. [35] has 
been developed in [37]. 5*0(5) symmetry [38] which unifies the antiferromagnetic and superconducting 
behaviour has also been used to explain pseudogap behaviour. Finally spin-charge separation [39] has 
been invoked. This last theory also predicts [40], in a quasi-2D system, a novel spin-metal phase with a 
nonzero spin diffusion constant at zero temperature. It is therefore entirely conceivable that the phase 
diagrams shown earlier are incomplete, and that new phases await discovery and/or confirmation, 
particularly if one applies a magnetic field (see for example [41]). 

1.3 Precursor superconducting fluctuations in particular phase fluctuations 

As mentioned above, one major thrust suggests that some form of precursor superconducting fluc- 
tuations or in particular phase fluctuations [25-27] are responsible for the pseudogap phenomenon. 
This assumption differs in principle from, for example, bipolaron theory of high-T c superconductivity 
(see the review [42] and textbook [43]), where a very strong electron-phonon coupling is assumed and 
where the bose-pairs are formed as stable units at temperatures higher than the temperature of the 
superconducting (in this case superfluid) transition. 

An incoherent pair tunnelling experiment [44] has been proposed which may be able to answer whether 
superconducting fluctuations are truly responsible for the pseudogap behaviour or whether another 
mechanism is involved. A second test of the nature of the condensation which comprises comparing 
two methods of determination of the gap has been proposed in [45]. The first method for finding 
the gap is related to the single-particle excitation energy A pair , measured by ARPES or STS. The 
second method is related to the energy coherence range A co h, defined by Andreev reflection. Finally, 
Andreev interferometry — the sensitivity of the tunneling current to spatial variations in the local 
superconducting order at an interface — has also been proposed recently [46] as a probe of the spatial 
structure of the phase correlations in the pseudogap state. 

One cannot however exclude the possibility that the pseudogap is the result of a combination of various 
factors, e.g. spin and superconducting (including phase) fluctuations^ (see, for example, models [47,48] 
where the interplay of antiferromagnetic and o?-wave pairing fluctuations has been recently studied). 

4 Indeed, the anisotropy of the spin subsystem in HTSC belongs to an easy-plane type which in turn implies 
that the magnetic vectors (averaged spins) can be described by a model similar to the XY model. The loss 
of long-range magnetic order under doping and the appearance of magnetic correlation length <^mag ^ x ' 
(see [22]) results in some sense in the formation of a gap ~ ^~^g i n the spin spectrum. Superconducting phase 
fluctuations, as will be seen in Sec. 5, can also be described by a XY-model and their phase transition behaviour 
is similar. 



Different authors argue that different types of superconducting fluctuations are responsible for the 
pseudogap. The simplest example of a theory which has the pseudogap/gap in the normal state is the 
so-called BCS-Bose crossover theory. 

Imagine starting from the weak coupling limit, which is well described by BCS theory [49], and gradually 
increasing the intensity of the attraction between the electrons. If the coupling is sufficiently strong, one 
no longer needs the presence of the Fermi surface which was essential in the BCS limit for the formation 
of Cooper pairs in metals [50]. The size of the Cooper pairs decreases until they can be regarded as 
almost structureless singlet bose-particles. The temperature of formation of these particles is given by 
the BCS pairing temperature, T pair which is determined by the strength of the coupling. However in 
contrast to the BCS case, these "pre-formed" bosons condense into a single superfluid quantum state 
not at T pair , but only at the temperature of Bose-Einstein condensation (BEC), T B <C T pair . It is obvious 
that even above Tb one will need to add extra energy A pair = A;BT pair to break up these pairs and this 
should lead to the gap-like features observed above Tb- This argument is strictly valid only in 3D. Note 
that in the space of lowered (2D or ID) dimensionality the formation of local pairs does not demand as 
strong an attraction as in 3D and the condensation process (see below) is not BEC. 

There is even some evidence [51,52] from time-resolved optical spectroscopy measurements of the quasi- 
particle relaxation dynamics that pre-formed pairs do exist in the underdoped YBa2Cu307_,5. We shall 
describe the history of the BCS-Bose crossover problem in Sec. 2 and consider a number of simple 2D 
models in Sec. 3 in the review since this problem appears to be very useful for a better understanding of 
superconducting fluctuations in general. Furthermore we believe that the BCS-Bose crossover problem 
is in its own right an exciting chapter in the history of research on superconductivity. 

It is, however, rather difficult to combine the arguments about the existence of the pre-formed local 
pairs with the fact that the one observes a gapped Fermi surface in HTSC [4,8]. In contrast the simple 
BCS-Bose crossover models (see for the reviews [53,54]) predict its absence in precisely the Bose limit 
of pre-formed pairs. It follows from these observations that the attraction between carriers in HTSC is 
not strong enough to destroy the Fermi surface, so that the pairs above T c appear to be short-lived (the 
same is true for bipolarons [42,43]). The presence of such resonant pairs above T c may still substantially 
affect the normal state properties of the system and this effect has been investigated in the papers of 
Levin et al. [55-57]. We also note that for spatially inhomogeneous systems it may be possible to have 
a well defined Fermi surface and yet regions of low carrier densities where pre-formed local pairs are 
present [58]. 

Up to this point our arguments were essentially independent of the dimensionality of the supercon- 
ductivity. Of course as mentioned above, if one considers the 3D case, stronger attraction is necessary 
to reach either the Bose limit or resonant pair regime than for the 2D system. It is known, however, 
that all HTSC are highly anisotropic quasi-2D systems with an almost 2D character of the conducting, 
superconducting and magnetic properties. Indeed because the coherence length is less than a lattice 
spacing in the direction perpendicular to the Cu02 planes (c-direction), the superconductivity in the 
copper oxides takes place mainly in the weakly connected interacting Cu0 2 layers (or their blocks). 
This is the reason why pure 2D models of HTSC are commonly accepted and studied. Undoubtedly the 
cuprate layers do exchange carriers even if these layers are situated in different unit cells (blocks). The 
mechanism for this transport (by coherent or incoherent (including pair) electronic tunnelling) is not 
yet established (see, for example, [6]). 

There is no doubt nonetheless that one must take into account the possibility of different (for instance, 
direct or indirect) interlayer hoppings to develop the full theory of HTSC. Strictly speaking therefore, one 



needs to consider quasi-2D models of these superconductors. In practice, however, predominantly either 
2D (or 3D) models have been considered with only a few attempts to consider quasi-2D models. This is, 
of course, related to the fact that the theoretical analysis even for the relatively simple attractive model 
is complicated. This is due to the necessity to go beyond the mean-field approximation, which proved 
good for BCS theory, and to take into account the well developed fluctuations of the superconducting 
order parameter. 

It can be seen from the anisotropy of the conductivity that the influence of the third dimension varies 
strongly from one family of cuprates to the other. For example, this anisotropy reaches 10 5 — 10 6 for the 
Bi- and Tl-based cuprates, while its value for YBaCu3C>6+5 compound is close to 10 2 . It seems therefore 
plausible that in the HTSC Bi(Tl)-compounds the transport in the c-direction is incoherent, while in 
the Y-ones this transport is coherent at least near the optimal doping. 

These brief remarks are intended to indicate how difficult it is to take into account the layered structure 
of HTSC. As a further example, we note that the interlayer tunnelling mechanism proposed by Anderson 
[59], and based on the layered structure of the cuprates, was considered as a leading candidate for the 
theory of HTSC superconductivity until a significant discrepancy with experiment was revealed [60]. 

One should also mention that the effect of interaction between the layers has been recently experimen- 
tally studied in [61] by intercalating an organic compound into bismuth-based cuprates. Even though 
the distance between layers increased dramatically the value of T c was almost unchanged from that 
for the pristine material. This provides unambiguous evidence that the superconductivity in layered 
cuprates is intrinsically of a 2D nature. 

The low-dimensionality of HTSC along with a relatively small (at least in comparison to ordinary 
metals) carrier density provide especially good conditions for the formation of different types of vortex 
excitations. Although the superconducting transition at T c itself belongs to d — 3 XY universality 
class, the large anisotropy of certain HTSC should lead to the Berezinskii-Kosterlitz-Thouless (BKT) 
[62] 2D XY regime over a very wide temperature range, until a crossover to 3D XY critical behaviour 
over a rather narrow temperature range (see, for example, the textbook [63]). One also expects that 
the transition temperature could well be practically unchanged [64], Tbkt < T c . Thus, outside the 
transition region the low-energy physics will be governed by the vortex fluctuations^ so one may 
expect the 2D models which we are going to discuss below to be especially relevant for the description 
of the pseudogap phase. Thus we arrive at the picture of Emery and Kivelson [25,26] (see also earlier 
paper by Doniach and Inui [66]) where the role of the phase fluctuations of the superconducting order 
parameter, particularly the vortex excitations is especially important. Indeed this picture has been 
recently supported experimentally [67,68] by the measurements of the screening and dissipation of a 
high-frequency electromagnetic field in Bi-cuprate films. These measurements provide evidence for a 
phase-fluctuation driven transition from the superconducting to normal state. 

Without pretending to present the whole picture (this is now practically impossible), we take the point 
of view that the effect of low dimensionality is leading to strong order-parameter phase fluctuations, and 
the dependence of the conducting properties on doping, are among the key ingredients of any proposed 
theory. We shall use general physical assumptions and relatively simple models (which permit much 



5 In a strictly 2D system the BKT transition is associated with proliferation of unbound vortex-antivortex pairs. 
Weak coupling between the planes leads to correlated motion between vortices in adjacent planes which form 
3D vortex loops close to the critical temperature. The transition to the disordered phase is then characterized 
by the appearance of vortex loops with arbitrarily large radii [65]. 



analytical calculation) to discuss the new properties of HTSC, including the pseudogap, based on the 
assumptions given above. As we have discussed above, there is no consensus about the origin of the 
pseudogap. We therefore believe that all currently existing theories of pseudogap behaviour should be 
developed until their predictions can both be experimentally tested and contrasted with the predictions 
from the competing theories. Thus the main goal of this review is to outline recent developments in the 
theory based on phase fluctuations. 

1.4 Outline 

This review focuses on the simplest pairing Hamiltonians but in the limit of two dimensions and low 
carrier density. As such we emphasize the separation into modulus and phase variables so essential in two 
dimensions. Not only are these models the simplest prototype Hamiltonians in which to discuss phase 
fluctuations and pseudogap behaviour, but as we shall demonstrate, illustrate many of the properties 
of the wider class of theories based on phase fluctuations. 

In Section 2 we consider a brief historical introduction to the BCS-Bose crossover as the simplest theory 
displaying a pseudogap. In particular we critically discuss the role of dimensionality and indicate why 
phase fluctuations are vital in the extreme two dimensional limit. In Section 3 we present the BCS-Bose 
crossover in a variety of models but at zero temperature where mean field theory still applies. In par- 
ticular we focus on how the inclusion of more realistic Hamiltonians influence the crossover. The effects 
discussed in this section include the role of <i-wave pairing, the role of the retardation of the interaction, 
and the effect of interlayer couplings. We stress that the theory of phase fluctuations is not in fact 
synonymous with the BCS-Bose crossover. In Section 4 we briefly sketch the corresponding results at 
finite temperature obtained using the T-matrix theory. These results have been derived predominantly 
in three dimensions and we indicate how this approach breaks down in the 2D limit. In Section 5 we 
present the BCS-Bose crossover in two dimensions at finite temperature. The superconducting transi- 
tion temperature in this case is the Berezinskii-Kosterlitz-Thouless transition temperature. Particular 
attention has been paid to the important effects of quantum phase fluctuations, Coulomb repulsions and 
the observed presence of non-magnetic impurities. In Sections 6 and 7 we construct the single particle 
Green's function and the corresponding spectral density which has been measured in recent ARPES 
experiments. We indicate the non-Fermi liquid behaviour induced by phase fluctuations and discuss 
how various mechanisms can "fill" the mean-field gap to give pseudogap behaviour. Section 8 presents 
our conclusions. 



2 A history of the BCS-Bose crossover problem 

2.1 Early history 

The idea that the composite bosons (or local pairs as they are sometimes called) exist and define the 
superconducting properties of metals is in fact more than 10 years older than the BCS theory. As early as 
1946, a sensational communication appeared saying that the chemist-experimentalist Ogg had observed 
superconductivity in the solution of Na in NH3 at 77K [69]. It is very interesting that the researcher 
made an attempt to interpret his own result in terms of the BEC of paired electrons. Unfortunately, the 



discovery was not confirmed and both it and his theoretical concept, were soon completely forgotten 
(for the details see e.g. [70,71]). 

It is instructive to note that the history of superconductivity has many similar examples. Probably, this 
explains why Bednorz and Miiller named their first paper "Possible high-T c superconductivity..." [1]. 
And even today there appear many unconfirmed communications about room-temperature supercon- 
ducting transitions. 

A new step in the development of the local pair concept, was taken in 1954 by Schafroth who in fact 
re-discovered the idea of the electronic quasi- molecules [72]. This idea was further developed in the 
Schafroth, Blatt and Butler theory of quasi-chemical equilibrium [73], where superconductivity was 
considered versus the BEC. Such a scenario, unfortunately, could not compete with the BCS one due 
to some mathematical difficulties which did not permit the authors to obtain the famous BCS results. 
Then the triumph of the BCS theory replaced the far more obvious concept of the local pairs and their 
BEC by the Cooper ones and their instability which takes place in the necessary presence of a Fermi 
surface (or more precisely a finite density of states at the Fermi level). In contrast for the local (i.e. 
separated) pairs the formation is not in principle connected to this density of electronic states. Also 
unlike the local pairs, the Cooper ones are highly overlapping in real space. More exactly there are no 
pairs in real space and the Cooper pairing should be understood as a momentum space pairing. 

Later experimental results reminded physicists of the existence of local pairs. Indeed, it was Frederikse 
et al. [74] who found that the superconducting compound SrTi0 3 (T c ~ 0.3K) has a relatively low 
density of carriers and, moreover, this density is controlled by Zr doping. The first deep discussion of 
the possibility of the BCS-Bose crossover and electron binding above the superconductivity transition 
temperature for a low density of carriers was carried out by Eagles [75] (see also his et al. relatively 
recent article [76]) in the context of SrTi0 3 . 

The subsequent history of the investigation of the BCS-Bose crossover has been often cited in the current 
literature. So, we note only that the features of 3D crossover at T = were considered in [77] (see also 
recent analytical investigation of the problem [78]), and its extension to finite temperatures was first 
given by Nozieres and Schmitt-Rink [79]. The 2D crossover in superfluid 3 He was studied in [80], where 
the very natural and convenient physical parameter the bound pair state energy, was first used. 

2.2 Relevance to HTSC 

The discovery of HTSC rekindled the interest in the problem of crossover and related phenomena. Let 
us discuss the reasons why the crossover problem appeared to be relevant to the general problem of 
the understanding of the HTSC. Indeed, these superconducting compounds have some peculiarities 
which place them much closer to the Bose or at least to the crossover region than the majority of low 
temperature superconductors. 

The ground state of the copper-oxide based materials forms due to the strong Hubbard repulsion 
resulting in antiferromagnetic spin fluctuations in the proximity of the metal-insulator transition (see 
e.g. reviews [22,81]). It seems plausible that the subsequent doping of these materials results in the 
appearance of weakly interacting (i.e. non-strongly correlated) itinerant carriers (holes). However the 
exact nature of the ground state including the strong electron-electron correlations in the presence of 
itinerant carriers is in fact not yet understood. We note that the experimental evidence for La2_ x Sr x Cu04 



indicates that the unusual "insulator-to-metal" crossover is even near optimal doping [82]. In addition 
the recent measurements of the of the transport properties in YE^CuaOz [83] also indicate that the 
metal-insulator transition takes place close to x — 6.42 which corresponds to the underdoped region. 
Thus it is an important challenge for theoreticians to include the metal-insulator transition into the 
theory of superconductivity. However currently most theoretical papers assume that superconductivity 
in HTSC develops on the metallic side of the metal-insulator transition and our review is not an 
exception. 

The density n/Q of these itinerant holes is not as large as in ordinary metals, so that the mean distance 
between them proves to be comparable with a pair size £ pa i r or a coherence length.Q This situation is 
significantly different from the conventional BCS theory where the parameter £ pair greatly exceeds the 

1 /2 

mean distance between carriers which is ~ . Experimentally the dimensionless value of /cp^pair (&f 
is the Fermi momentum) which describes the ratio of the pair size and the distance between carriers is 
about 5-20 for HTSC while for the low-temperature superconductors it is about 10 3 - 10 4 [22,53]. 

The quantitative differences between HTSC and low temperature superconductors were summarized on 
the well-known Uemura plot [84], where the superconductors are classified by the ratio of their critical 
temperature, T c and the superfluid density expressed in terms of the effective Fermi energy, ep. One 
can see from Fig. 3 that cuprates, organic superconductors and some other "exotic" superconductors 
have ksTc/eF as high as 0.01-0.1, much higher than those of conventional superconductors. This value 
of T c is however 4-5 times less than T-q. Another important issue which follows from Fig. 3 is the linear 
relationship T c ~ n s /m*. Uemura has interpreted this dependence as originating from BEC, but we 
would like to note that this dependence may also be understood within the phase fluctuation scenario. 
This linear dependence is absent in overdoped cuprates where the depression of T c is associated with 
a decrease of the superconducting condensate density n s in spite of the increasing normal-state carrier 
density [85]. 

In fact, as we will see later, the new materials are likely to be in an intermediate regime between the 
Cooper pairs and the composite bosons, at least when the doping is not large and the value of T c is far 
from the highest possible (optimal) one. 

Based on the phenomenological classification presented in Fig. 3 Uemura [84] suggested the following 
interpretation of the phase diagram from Fig. 1 in terms of evolution from BEC (in real space with non- 
retarded strong interaction) to BCS condensation (in momentum space with retarded weak interaction) 
shown in Fig. 4. 

There is an important outstanding question as to how the picture proposed by Uemura correlates with 
the ideas suggested by Emery and Kivelson [25,26] and we will discuss this point in Sec. 5.4. 



6 One should distinguish between the doping x which is usually connected to the chemical formula of the 
compound and nf which is defined as the carrier density in each separate layer and because of this can depend 
on the number of layers per unit cell, chemical interactions, etc. In particular, nj = n 2 ® = n^ D c, where n 3D is 
the 3D (bulk) carrier density and c is the lattice constant. 

7 It is good to bear in mind that the coherence length strictly speaking is distinguished from the pair size, 
especially at low carrier density, since the first one is related to the energy coherence range A co h and the second 
to the single-particle excitation energy A pa j r [45] (this question will be treated in Section 3.1.5). 



(D 

Z5 

O 

CD 
Cl 

£ 

(D 
I— 

C 

o 

'(f) 
c 
o 



n s 2/ V(mVm e ) or 0.66n s2d /(mVm e ) (cm 



10 



12 



10 



13 



10 



s2d/ 
14 



10 



15 



2223 



10 



10' 




10" 



10 



KB Q V 3 Si 



Chevrel 
'ET O c 
NbSe 2 
U R Fe# 



A Nb 



UPd 2 AI 3 
TMTSF^ * 

URu 2 Si 2 

UBe 13 
-•ft— - 
UPt, 



Sn A 



Al 
A 
Zn A 



10 10 10 10 10' 

Fermi Temperature T F = e F /k B 



10 



16 



K 



Fig. 3. A plot of T c versus the effective Fermi energy, ep ~ n&d/m*, where n S 2d is the 2D superfluid density 
and to* is the effective mass of the carriers in Cu02 planes. Tb denotes BEC temperature for non-interacting 
bosons with density n s /2 and mass 2m*. This plot was taken from [84]. 



2.3 BCS-Bose crossover and pseudogap: survey of the current literature 



There are so many published articles dealing with the BCS-Bose crossover study that we are only able 
to mention the review papers [53,54] where the interested reader may find more references. 

We also note the existence of the comprehensive review of Ginzburg [86] (see also [71]) which presents 
a historical overview of the parallel (and generally speaking independent) development of the macro- 
and microphysics of superfluidity and superconductivity. 

In the present review we focus on the more recent results related to the crossover problem, and with 
a stronger accent on the pseudogap phenomena. In most cases the attractive (negative-?/) 2D or 3D 
Hubbard model has been considered (see the review [87] for the results obtained prior to 1990). In 
particular this model has been studied within the "<3> derivable" in the sense of Baym [88] , conserving 
[89] and self-consistent T-matrix approximation both analytically [90-92] and numerically [93-101]. In 
all of the above papers the 2D model was considered with the exception of [90,96], where the 3D model 
was studied and [91] where the influence of the third direction was introduced. 

The non "$ derivable", but still conserving T-matrix approximation was considered in [55-57]. As we 
have already mentioned, the pseudogap in this approach has been related to the resonant pair scattering 
of correlated electrons above T c . This approach regards the pseudogap as a 3D phenomenon. 
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Fig. 4. Phase diagram describing BEC-BCS crossover with increasing carrier concentration n. The phase di- 
agram can be mapped to that of cuprates assuming that the pseudogap temperature T* corresponds to the 
formation of normal state pairs. This diagram was taken from [84]. 

A comprehensive critical analysis of the self-consistent T-matrix approximation which is, as we can 
see, widely utilized in the literature was performed for the 3D case in [102]. The authors came to the 
conclusion that a new class of the diagrams has to be included. In particular, they demonstrate the 
importance of the correct treatment of the residual interaction between the composite bosons in order 
to regularize the strong coupling limit. It is claimed [102] that the self-consistency of the fermionic 
Green's function in the T-matrix approximation maybe less important than due care about the residual 
bose-bose interaction. Their conclusion (in the 3D case) that the T-matrix approximation is in fact 
inadequate underlines our strong reservations about the applicability of the T-matrix approximation in 
the 2D limit (see Sec. 4.5). 

The pseudogap for the case of <i-wave pairing has also been studied in [103-106] (for the review see 
[107]) and quantum Monte-Carlo simulations for the 2D attractive Hubbard model were performed in 



There is in fact a big difference between the 2D and 3D theories mentioned here. In particular, the 
pseudogap obtained in the 3D approach [55,57] is the result of rather strong coupling, while to obtain 
the pseudogap in 2D theory one needs intermediate or even weak coupling [112,114]. Thus the 3D 
approach [55] has far stronger constraints on the mechanism of HTSC, than do the studies of the 2D 
models mentioned here (see also the papers [114-125] which we will discuss in the subsequent chapters 
of the review). These show, that almost independently of the mechanism of HTSC and the intensity of 
coupling, the fluctuations in the order parameter phase contribute to the pseudogap. Even if the origin 
of the pseudogap is related to some other mechanism e.g. those discussed in Sec. 1.2, the preconditions 
for the existence of a superconducting phase contribution to the pseudogap are so easily satisfied that 
such a contribution must always be present. 



[108-113]. 



It is assumed in 3D [55,57] and quasi-2D [56] models that T c is the temperature for the formation of long- 



range order which must go to zero in the 2D limit in accordance with the Coleman-Mermin-Wagner- 
Hohenberg (CMWH) theorem [126] .0 Indeed the theorem prohibits the formation of a superconducting 
phase with homogeneous order parameter (the latter related to breaking a continuous symmetry) in 
2D. This is due to the fact that the fluctuations of the order parameter phase destroy the long-range 
order [128] (see also [129]). Thus even this very general argument as to why long-range order is absent 
in 2D leads one to the conclusion that the phase fluctuations play a significant role in low dimensional 
models of metals like HTSC. 

The ideology of the 2D theories is completely different from those in 3D. As we have just explained, in 
2D there is no superconducting transition into a state with long-range order, and the only possible kind 
of transition is the BKT one [62]. (Since the physics of the BKT transition is rather subtle, and may be 
unfamiliar, we bring to the attention of the reader the following review articles [130,131] and textbooks 
[132-134]). Thus the temperature T c in these models is identified with the temperature of the BKT 
transition, Tbkt below which an algebraic order is established in the system. This transition does not 
require symmetry breaking and is therefore not forbidden by the CMWH [126,127] theorem. Although, 
due to the presence of even weak interlayer tunnelling, the transition is of d = 3 XY type with true 
long-range order, the critical temperature T c is expected to be close to the value of Tbkt calculated for 
the pure 2D model [64]. 

An approach to the 2D single-band repulsive and attractive Hubbard models which allows one to 
automatically satisfy the CMWH theorem, and enforce the conservations laws, the Pauli principle and a 
number of sum rules has been discussed in [135]. We refer also to another approach [136] developed for the 
Hubbard model with a strong repulsion which results in indirect attraction between carriers. It is worth 
to mention this approach here since it treats the emerging superconducting phase fluctuations in the way 
very similar to the method presented in the review. The phenomenology of superconductivity resulting 
from the Bose condensation of pre-formed pairs coexisting with unpaired fermions has been studied 
in [137]. It was shown that the phenomenology describes reasonably well the data in the underdoped 
Y-Ba-Cu-0 family but fails to describe these compounds in the optimally doped regime or underdoped 
La-Sr-Cu-O. The influence of the existence of the fermion bound states above T c was also considered 
in [138] and the theory of chemical equilibrium between these bound states and single fermions was 
applied in [139] to explain the normal state properties of HTSC (see also [140] which considers the BEC 
of such pairs). 

The large number of references cited here shows how popular this topic is currently and the authors 
apologize in advance to those whose work has been omitted. 

Certainly, it is necessary to point out immediately that the crossover phenomena in particular and the 
investigation of the attractive Hubbard model in general do not (and cannot) address the problem of the 
mechanism for HTSC. The problem of the mechanism for HTSC is very difficult and still controversial. 
Nevertheless the treatment of the crossover or, more generally, the low dimensional models with an a 
priori postulated attraction between carriers may shed light on some of the features of HTSC. 



We note that the CMWH theorem has been recently revisited for the special cases of ID and 2D Hubbard 
(both attractive and repulsive) models [127]. 



3 BCS-Bose crossover in 2D systems at T = 

This chapter is devoted to the description of the BCS-Bose crossover in a variety of 2D systems at 
T = 0. The reason why we decided to restrict our initial consideration to the zero temperature limit is 
to avoid the complexities related to the CMWH theorem [126] . Indeed, at T = the problem becomes, 
due to the integration over frequency (which replaces the summation over Matsubara frequencies for 
T 7^ 0), effectively a 3D one. Therefore, the 2D theorems [126] are not applicable to this case and one 
can safely consider long-range order in 2D systems at T = 0. Thus one might expect at T = the 
simplest BCS approximation with minor modifications to be able to adequately describe the physics of 
the BCS-Bose crossover. 



3.1 One band continuum model with s-wave pairing and non-retarded attraction 

Let us introduce a continuum field theoretical model of fermions with an attractive two-body interac- 
tion.^] Our goal is to consider how the energy gap A, the chemical potential /x, the coherence length £ co h 
and the pair size £ pa i r change as a function of the density of bare fermions, rif . We shall follow here the 
papers [54,141] (see also the original paper [80], the subsequent paper and review [142,53]) and discuss 
the difference between the coherence length and Cooper pair size using Refs. [143-146]. 

The simplest model is described by the Hamiltonian density, which is the field-theoretical analog of the 
negative-?/ Hubbard model (see Sec. 3.3), 



where x = r, r (r is a 2D vector and r is the imaginary time); ip a (x) is the Fermi field; m is the fermion 
effective mass; a =|, I is the fermion spin; U > is the attraction constant. The chemical potential fi 
fixes the average density rif of the free (bare) carriers. We choose units in which h = = 1 and the 
system occupies the volume v. 



9 It is worthwhile to emphasize once more that the application of a weak coupling model for the description 
of HTSC suggests weak coupling for the doped carriers only. Indeed, as stated above, HTSC compounds are 
antiferromagnetic insulators with a strong onsite (Hubbard) repulsion on Cu ions that enforces an almost local 
distribution of the spin density on these ions (localized fermions with spin S = 1/2, or copper 3d 9 configuration). 
This initial system should undoubtedly be described by strong coupling models near half-filling (n = 1). The 
doping of the system decreases the number of electrons on the copper ions (and also on the oxygen ions, if a 
pd- hybridization model is used), introducing the positively charged delocalized carriers (holes) which move in 
the medium with a short-range magnetic order. According to our assumption these current-carrying holes are 
weakly coupled and form correspondingly a small Fermi surface. These arguments justify to some extent why the 
simplest 2D models with a direct (negative-C/ Hubbard model) or indirect (for example, phonon-like) attraction 
are considered below to establish how the variable carrier density and the condensate phase fluctuations affect 
the properties of 2D superconductivity. One can hope that, under these circumstances, many of the physical 
properties described below persist in more realistic approaches and models. 




(3.1) 



3.1.1 Formalism 



The functional integral approach is appropriate for the problem studied here and will prove very useful in 
the subsequent chapters. Since we will need the results presented here in subsequent chapters (including 
the case of nonzero temperatures), we are going to use the Matsubara thermal technique from the 
beginning, but will take the limit T = in this chapter. Thus, let us consider the formalism used. 

Introducing the Nambu spinors for the fermion fields [147] (see also the textbooks [49,148]), 

= ] > &w=U\(x) mo), ( 3 - 2 ) 

one should rewrite (3.1) in the appropriate form: 

H = -^(x) (¥- + fj\ t 3 ^(x) - U&(x)t+V(x)V*(x)t-V(x), (3.3) 

\ Ziffh J 

where r 3 , t± = (ti ± ir 2 )/2 are Pauli matrices. 

Now the partition function can be expressed through the Hamiltonian (3.3) as: 

Z(v,fi,T) = J PtfP^exp j-J^dT J d 2 r[* t (a;)9 T *(a;)+W(r)]|, (3.4) 

where f3 = 1/T and D^T>^ denotes the measure of the integration over the Grassmann variables \I/ and 
satisfying the anti periodic boundary conditions: \I/(t, r) = — \P(r+/3, r) and ^(r, r) = — \Ld(r+/3, r). 

If it were possible to calculate the partition function (3.4) exactly one could obtain all thermodynamical 
functions from the thermodynamical potential 

Q(v, fi, T) = -T In Z(v, fi, T). (3.5) 

Introducing the auxiliary Hubbard-Stratonovich complex scalar field in the usual way one can represent 
(3.4) in an exactly equivalent form (see the review [149]): 



Z(v, fi,T) = J VW&V$V$* exp l-J^drjd 



2 r 



Mx)l 
u 



dj ~ r 3 + f?j - t + $(x) - r_$*(z) 



(3.6) 



The main virtue of this representation is the non-perturbative introduction of the composite fields 
$(x) = U^(x)t-V(x) = Uip^tp^x), $*(z) = U&(x)t+V(x) = Ui>\(x)i>\(x) and the possibility to 
develop a consistent approach. Specifically, the expression (3.6) turns out to be rather convenient for 
studying non-perturbative phenomena such as superconductivity. In this case the complex Hubbard- 
Stratonovich field naturally describes the order parameter arising due to the formation of the Cooper 
or the local pairs. The average value of |$|(= A) is proportional to the density of pairs, on the one 
hand, and determines the gap in the one-particle Fermi-spectrum, on the other. 



The integration over the fermion fields in (3.6) can be done formally even though $ and $* depend on 
the spatial and temporal coordinates. Thus, one obtains (formally exactly) 



Z(y, /i, T) = J exp[-pn(v, fi, T, $*(x))], 



(3.7) 



where 



pn(v,(i,T,Q(x),Q*(x)) = \-f dT \ d 2 r\$ 

U Jo J 



'x)\ z — TrLnG -1 + TrLnG, 



i-l 



(3.8) 



is the one-loop effective action. This action includes in itself a series of terms containing derivatives 
with respect to <&(#) and <&*(#). In the lowest orders it corresponds to the Ginzburg-Landau effective 
action. 

The operation Tr in (3.8) is taken with respect to the space r, the imaginary time r and the Nambu 
indices. The action (3.8) is expressed through the fermion Green's function which obeys the equation: 



-Id T + r 3 + fij + r + $(r, r) + r_$*(r, r) 



G{t,t) = 6{t)6{t) 



with boundary condition 



The free Green's function 



G(t + (3,t) = -G(t,v). 



G (t,t) = G(r,r)| $i $ v=0 



(3.9) 



(3.10) 



(3.11) 



in (3.8) is needed to provide the regularization in the calculation of 

Q(y, fi, T, $, $*). The representation (3.7), (3.8) is exact, although to perform the calculation in practice 
it is necessary to restrict ourselves to some approximation. Below we shall use the assumption, that one 
neglect the fluctuations of the fields $(x) and <&*(x) and thus replace them by their equilibrium values 
which will now also be denoted as $ and $*. This saddle-point approximation is in fact absolutely 
equivalent to the BCS mean-field approximation and works quite well for T = even in 2D. We note 
that the assumption that $ = const ^ implies that we are considering a state with long-range 
order which is of course invalid for T ^ 0. We will return to this problem in Chapter 5, where the 
generalization for finite T will be done consistently. 

The thermodynamical potential Q as well as the partition sum Z now depend on $ and $* which play 
the role of the order parameter. The order parameter appears due to the fact that one is not considering 
the exact potential Q(v,fi,T) which contains a functional integral over the auxiliary fields. Instead one 
is considering its saddle-point approximation Q(v, /i, T, <&*): 



&pot(v, /!, T, $, $*) = 

Qj J d 2 r\$(x)\ 2 -TTrLnCT 1 + TTrLnG 1 



Q=Q*=const 



(3.12) 



in which the integral over the auxiliary fields has been replaced by a single term in which the auxiliary 
fields take on their equilibrium values $ and $*, and about which one can subsequently expand to 
obtain better approximations. 



To avoid possible misunderstanding, we shall write down the formulae for the Fourier transformations 
which are used throughout the review. They connect the coordinate and momentum representations in 
the usual manner: 



F(iu n , k) = J dr J d 2 rF{r, r) exp(icu n r — ikr) , 



(3.13) 



F(t, r) = T J2 J j2nf F ^ iUn ' k ) ex PH w « r + ikr ). 



(3.14) 



where cu n = irT(2n + 1) are the fermion (odd) Matsubara frequencies. In the case of bosons the odd 
frequencies should be replaced by even ones: Q n = 2nnT. 



For example, the Green's function (3.9) has, in the momentum representation, the following form: 

iuoj + r 3 £(k) - r+$ - r_$* 



G(iu> n , k) = 



^ + e 2 (k) + |$| 



(3.15) 



where £(k) = e(k) - /i with e(k) = k 2 /2m, $ = ($(x)) and $* = are already taken to be 

constants which represent the complex order parameter (or more precisely, the complex ordering field, 
because in 2D at T ^ there is no order parameter in its usual sense). 

Substituting (3.15) into (3.12), one arrives at (see Appendix A and Refs. [141,54,115]) 



n pot (v,n,T, $,$*) = 
M$| 2 r d 2 k 



V J (2tt) 2 
d 2 k 



2Tln cosh 



2T In cosh 
e(k) 



^(k) + |$ 



2T 



-£(k) 



+ 



2T 



-e(k) 



(3.16) 



3.1.2 The Effective Action and Potential 

As has been noted above, it is impossible to calculate (3.8) for <3> dependent on x in general. However, 
if one assumes that the gradients of $ and <£>* are small, the action (3.8) can be naturally divided into 
kinetic and potential parts 



Q(v, /i, $*(a;)) = n kin (v, /I, $(x), $*(x)) + n pot (v, /I, $, $*), 



(3.17) 



where the effective potential has been defined by (3.12) (its final expression is given by (3.16)). The 
terms flki n (^(x),^*(x)) with derivatives in expansion (3.17) contain important physical information, 
which we shall consider in Sec. 3.1.5. 



Let us return to the effective potential. In the limit T — > (3.16) reduces to 

|$| 2 r d 2 k 



n pot (v,ii,$,$* 



U 



(3.18) 



where the terms which do not depend on $, $* and /i have been omitted. 



It is interesting that, by virtue of the invariance of the partition function (3.6) with respect to the phase 
transformation of the group U(l) 



$(x), 



^ f (x) 
&(x) 



$*(x)e 



-tars . 



-2ia 



(3.19) 



with real a the potential fl pot (v, /i, $, $*) (3.12) (see also (3.16) and (3.18)) can only be dependent on 
the invariant product $*$[ Ttr [ 

The analytic solution of the problem for the 2D case that we consider here is easier than the 3D one [78] 
which is obtained in terms of special functions (compare with Eqs. (3.25) and (3.27) below). Indeed, 
after performing the integration over k in (3.18) (which can be done straightforwardly due to the energy 
independence of the free fermion density of states) one obtains 



1 



rn 



U Att 



In 



W-n+ J(W -v) 2 + |$| 2 



/i 2 + |$| 



+ 



W-n + J(W-/j,) 2 + |$| 2 



fi 2 + |$| 2 - jJL 



(3.20) 



where the value W = k 2 3 /2m is the conduction bandwidth and k# is the maximal (in some sense 
Brillouin boundary) momentum. 



3.1.3 Main equations and analysis of solution 

If the quantity Af"""] is defined as the average value of |$|, then the equation for the extremum 

dn pot (v, 



<S>=$*=A 



yields, according to (3.20) 



9$ 



1 m W - fi- 
U~4^ 11 VAi" + A 2 - /i 



0. 



(3.21) 



(W - jj) 2 + A 2 



(3.22) 



while the condition 



1 dQ pot (v,fi,<5>,$* 



v 



dfx 



n f , 



(3.23) 



10 There is another transformation (when the sign of the phase a is defined by the fermion spin rather than 
the charge) under which the Hamiltonian (3.1) (or (3.3)) is also invariant. Such a transformation proves to be 
important for fermion- fermion repulsion (i.e. U < 0), or for the fermion-antifermion (electron- hole) channel 
of pairing. Apart from this difference the formalism for the case of a repulsive interaction is identical to that 
under consideration. The complete set of gauge transformations for the Hamiltonian under consideration were 
originally given by Nambu [147]. 

11 This is the parameter that is responsible for the appearance of a new (ordered, or with lowered symmetry) 
phase (see Sec. 3.1.1) which is permissible at T = 0. 



which sets the density of the particles in the system, takes the form 



W - ^{W - /i) 2 + A 2 + ^ 2 + A 2 = 2e F} (3.24) 

where we have used the simplest quadratic dispersion law e F = irnf/m appropriate to 2D metals. 

Equations (3.22) and (3.24), which were obtained in the mean field approximation (as is well adequate 
for T = 0)|33, f° rm a complete set for finding the quantities A and \i as functions of W and e F (or 
n/). It differs from the similar set in [142] by the explicit dependence on W that, in principle, can be 
important for the case of narrow or multi-band systems which we are going to discuss in Sec. 3.2. 

It should also be noted that the need to use the system of equations to find A and \x self-consistently has 
been known for a long time (see [49]). However, in ordinary 3D metals the number density is practically 
unchanged, so that, as a rule, the equation for /i is trivialised to the equation \i = e F and only the value 
of A is regarded as unknown. The importance of the second equation for small particle densities was 
first pointed out in the papers [75,77]. 

Equations (3.22) and (3.24) allow, along with the trivial solution (A = 0, fi = e F ) the nontrivial one: 



A 5 
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e F {W-e F ) 
sinh 2 (27r/m[/) 

e f coth 
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W ( , 2tt 

— coth 

2 V mU 



1 



(3.25) 



which is valid for any physically reasonable values of the relevant parameters. It is clear that this solution 
is meaningful physically when U > only. A solution of precisely this form was originally obtained by 
[116] (see Eq. (5.50) in Sec. 5.2). The functions A 2 (ei?) and n(e F ) are shown in Fig. 5, where the gap 





Fig. 5. The functions A 2 (e F ) and fi(e F ) that follow from the solution (3.25) as taken from [141]. 

value attains its maximum at e F = W/2, here A max = (W/2) sinh -1 (27r/[/m). At the same point the 
value fi = W/2 > 0. The symmetry of the functions A 2 (ep) and fi(e F ) about the line e F = W/2 reflects 
the particle-hole symmetry about half-filling. One should consider the hole or anti-particle picture in 
the region W/2 < e F < W and thus change the sign of the mass m — > — m. It is also very important 
(see below) that for small e F , there is a region where \i < 0, and that the sign changeover occurs at a 
definite point e F = W/2[l — tanh(27r/m£/)]. 
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Of course, there are also residual interactions related to quantum fluctuations, etc. 



The expressions that are found in [80,142] follow directly from (3.25) if, treating W as large and the 
attraction U as small, we introduce the 2D two-body binding energy 



Eb = — 2Wexp 



mil 



(3.26) 



which does not include any many-particle effects. The introduction of the expression (3.26) enables one 
to take the limit W ^> ef and thus justifies to a certain degree the use of the parabolic dispersion 
law. In addition the fitting parameter Eb is more physically relevant. For example it is well-defined 
even for potentials with repulsion. We stress that the introduction of Eb instead of U also permits the 
regularization of the ultraviolet divergence, which is in fact present in the gap equation (3.22). 

It should be mentioned that, in a dilute gas model, the existence of a two-body bound state in vacuum 
is a necessary (and sufficient) condition for a Cooper instability [53,80,142]. This statement becomes 
nontrivial if one considers two-body potentials U(r) with short-range repulsion (e.g., hard-core plus 
long-range attraction), so that one has to cross a finite (but really very weak [151]) threshold in the 
attraction before a bound state forms in vacuum. 



where the chemical potential now changes sign at €f = |£ft|/2. 

To understand the physical significance of these remarkably simple results we look at the two limits of 
this solution. For very weak attraction (or high density) the two-particle binding energy is extremely 
small, i.e. \eb\ <C €f, and it is seen that we recover the well-known BCS results with strongly overlapping 
(in r-space) Cooper pairs. The chemical potential /i ~ €f, and the gap function A <C ef- 

In the opposite limit of very low particle density (or a very strong attraction) we have a deep two-body 
bound state |£&| ef, and find that we are in a regime in which there is BEC of composite bosons, 
or "diatomic molecules". The chemical potential here fi ~ — |e&|/2, which is one half the energy of pair 
dissociation for tightly bound (local) pairs. 

It should also be kept in mind that in the local pair regime (// < 0) the gap E gap in the quasi-particle 
excitation spectrum equals not A (as in the case \x > 0) but rather V + A 2 (see [77], the review [53] 
and the discussion of the Bethe-Salpeter equation in Sec. 3.1.4). 

Leaving aside the analysis of Q pot for arbitrary values of the parameters (see [141,54]), let us consider 
the most interesting case of W — > oo, U ^ with finite e b . Then finding from (3.26) the expression for 
A-n/mU and substituting it into (3.20) one obtains 



n pot (v, n, $*) = 




(3.28) 



Using the potential (3.20) or (3.28) one can demonstrate that the superconducting state is indeed 
energetically favoured. 



Making use of (3.26), it is easy to simplify (3.25) [53,80,142] to 





(3.27) 



The equations (3.22) and (3.24) can also be written as 

yWA 2 -^ |e 6 |; y^ 2 + A 2 + n = 2e F , (3.29) 
respectively; their solution, (3.27), is quoted above. 

3.1-4 Pairs in the Fermi sea 

To study the possibility of the existence of the bound states in the presence of a collection of fermions 
(Fermi sea) the Bethe-Salpeter equation was analyzed in [141]. Here we shall only present the results 
of this analysis in the particular limit W — » 00, U —>■ with finite £&-H The analysis in [141] showed 
that if the charge symmetry is unbroken, i.e. <3> = $* = 0, the Bethe-Salpeter energies are given by 

^=-¥{^^1} < 3 - 3o > 

Thus the system contains non-decaying bound states, i.e. the energies are real, only if ep < e c p = 
I £6 1/8. When ep > ef the solutions acquire an imaginary part which, according to [49], suggests the 
development of a Cooper instability and vacuum rearrangement. 

For the state with the rearranged vacuum so that $ = $* = A the analysis of the Bethe-Salpeter 
equation firstly proved the existence of a gapless Goldstone mode. The second solution of the equation 
which is only real for /i < is in agreement with the quasi-particle excitation energy E gap described 
above. 



3.1.5 The gradient terms of the effective action and the correlation length versus doping 

Now we calculate the terms Qkin which contain the derivatives in the expansion (3.17). As before, we 
shall assume the inhomogeneities of $ and $* to be small, restricting ourselves to terms with only the 
lowest order derivatives. For simplicity we shall also consider the stationary case and calculate only the 
terms with second-order spatial derivatives. These terms make it possible to determine the coherence 
length £ co h and its doping dependence. 

With these restrictions, and taking into account the invariance of £l(v, /j,, <&*) (see equation (3.8)) 
with respect to the phase transformations (3.19), one can write the most general form for the kinetic 
part of the action as 



n Hn (v,fi,T,^(x),^(x)) =T j^dr j d 2 r T fcm ($, $*, V$, V$*) = 

(3.31) 



T ( dr I d 2 r 



ri(|$| 2 )|v$| 2 + ^r 2 (|$| 2 )(v|$| 2 ) 2 



13 It is important to stress the difference between the bound states in the vacuum which are always present in 
a 2D system with attraction and the possibility of existence of such states in the presence of other particles. 
To analyse the conditions for the existence of bound states in the latter case, the Bethe-Salpeter equation (see, 
for example, [152]) must be studied. 



Here there are no items with a total derivative since boundary effects are regarded as unessential, and 
the coefficients T 1)2 (|$| 2 ) are assumed to be unknown quantities. It follows from (3.31) that, if one 
includes derivatives up to second order, variations in both the direction (phase) of the field $ and its 
absolute value are included. 

The calculation for the coefficients T lj2 (|$| 2 ) is straightforward but rather lengthy. It has been carried 
out in [141] (see also [54]) using a general method for derivative expansion developed in [153]. The paper 
[141] also contain the expressions for the case of a finite band-width W. 

Knowing Ti(|$| 2 ) and T 2 (|$| 2 ) one can find the values for different observables. For practical purposes 
one can restrict oneself to considering the coefficients obtained at the point of minimal effective potential 
|$| = A. Instead of ( | <I> | 2 ) and T 2 (|$| 2 ) it proves to be more convenient to introduce the combination 
T 2 (|$| 2 ) = T\ ( | <5 1 2 ) + 2A 2 T 2 (|$| 2 ) which determines the change in the |$| value only and arises as the 
coefficient at (V|$|) 2 . Furthermore, using (3.26) and (3.27), one can derive 7\(A 2 ) and T 2 (A 2 ) in the 
very simple form: 



T\(A 2 



1 



lGirh 2 \e b \ 



T 2 (A 2 



1 (2e F -\e b \) 2 
2Aith 2 (2e F + \e b \) 3 ' 



(3.32) 



where for completeness we have restored Planck's constant. 



The explicit forms of T\(A 2 ) and a = v 1 d 2 Q pot ($, $*)/9$(9<l>*|| (I ,|2 =A 2 permit one to calculate the 
coherence length and the penetration depth. From Eq. (3.28) one has 



m 



2nh 2 2e F + \e b [ 

Then, according to general theory of fluctuation phenomena [154], one obtains 



(3.33) 
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(3.34) 



This formula shows the dependence of £ co h on (or rif). The zero temperature coherence length £ co h 
was also studied in [144] for 2D and 3D cases, where it is referred as the phase coherence length £ p h aS e- 

It is very interesting and useful to compare (3.34) with the definition of the pair size [145,146] (also 
sometimes incorrectly referred to as the coherence length (see [142,143])), namely 



3 pair 



Jd 2 rg(r)r 2 
J d 2 rg(r) ' 



(3.35) 



where 



9K T ) 



(^BCs|^}(r)Vl(0)|#Bcs> 



(3.36) 



is the pair-correlation-function for opposite spins and |^bcs) is the usual BCS trial function. For the 
2D model under consideration the general expression (3.35) gives [142,145] 
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(3.37) 



where A and /x are given by (3.27). So, we are ready now to compare (3.34) and (3.37). 



Again to understand the underlying physics it is worthwhile to look at the two extremes of (3.34) and 
(3.37). For high carrier densities, e F |£&|, one finds that £ cori ~ £ pa ir ~ Hv F /A, i.e. the well-known 
Pippard's result is reproduced correctly. Moreover, if we introduce the pair size = fo(m|£b|) -1 / 2 in 
vacuum [80], it is clear that in the high density limit < e F and /i m €f), £coh ~ £ pa ir < £&■ Thus both 
£ co h and £ pair prove to be of the order of the pair size, £&, which is much larger than the interparticle 
spacing. The latter statement follows from the value of the dimensionless parameter £ co h^F ~ £ pa ir^F ~ 
e F /A > 1. 

In the opposite limit of very low density (|£ 6 | ^> e F and /i < 0) one can see that <C £ cori ~ kp 1 while 
Cpair ~ £b- Consequently the correct interpretation of £ pa i r is the pair size (in presence of the Fermi sea 
(3.35)) rather than the coherence length. The former in the extreme Bose regime is much smaller than 
the mean interparticle spacing, i.e. £ pa ir^F 1) while £ C oh^F ~ 1- The meaning of £ cori is the coherence 
length because it remains finite and comparable with the mean interparticle spacing even in the limit 
of infinite binding i.e. when \ei,\ goes to infinity. This situation is consistent with the case of 4 He where 
the coherence length is nonzero and comparable with the mean inter-atomic distance although |e 6 | (or 
energy of nucleon-nucleon binding) is extremely large. 

Pistolesi and Strinati in [144] obtained the same results, namely the coincidence of the pair size, £ pa j r 
and the coherence length £ coh in the weak coupling limit £ pa i r /cF 1 and the inequality £ cori £ pair in 
the strong coupling or Bose limit £ pa i r A:F 1. They established that £ cori coincides with £ pair down to 

£p a ir^F ^10. 

The above discussion is a clear example of how the two different energetic scales discussed in [45] behave 
in the low- and high-density limits. Indeed, the single-particle excitation energy, A pair = A is obviously 
related to the individual pair size, £ pa i r , while the energy coherence range is defined by the coherence 
length, £ coh . It follows from the consideration of £ pair and £ cori above that the energy scales A pair and 
A cori are the same in the high-density limit. However they diverge in the low density limit where A pair 
is larger than A co h. 

It is interesting to note that the femtosecond time-domain spectroscopy [58] (see also the previous papers 
on this subject [155] and the papers discussed above [51,52]) shows the existence of two distinct gaps 
in the entire overdoped region of Y 1 _ 2 ,Ca x Ba2Cu307_ < 5, where the doping here refers to the chemical 
fraction of Ca included. It is claimed in [58] that one of them is a temperature-independent "pseudo- 
gap" A pa i r and the other is a T-dependent collective gap A coh (T). They suggest that the presence of 
two gaps is due to a spatially inhomogeneous picture, supported experimentally by the observation of 
an inhomogeneous charge distribution [156]. In low carrier density regions A pair is simply the energy 
scale for single pair formation. In high-density regions one has strong collective effects and a collective 
temperature dependent gap A coh (T) which shows BCS-type closure at T c . If such a two-phase picture is 
indeed correct, then one may have a Fermi surface in that the chemical potential remains positive yet 
at the same time in local low carrier density regions one can reach the Bose limit - limit of pre-formed 
pairs. 

Another important comment which we have to make here is that the analysis of the experimental 
data performed in [145] shows that the optimally doped cuprates are definitely on the BCS side of the 
BEC-BCS crossover in the sense that the chemical potential is close to e F if one assumes a spatially 
homogeneous system. To make such a conclusion the ratio fi/e F was extracted in [145] from the available 
experimental data for the coherence length. In particular, it was crudely estimated that for YBaCuO 



with T c = 93K the ratio /i/e F ~ 0.9998 - 0.9995; for BiSrCaCuO with T c = 100K — fj,/e F ~ 0.9991 - 
0.9978 and for TIBaCaCuO with T c = 125K — fi/e F ~ 0.9986 - 0.9965. It is clear that for a spatially 
homogeneous model that closeness to the Bose limit implies negative ratio n/e F or at least /i < e F 
which is apparently not the case for optimally doped HTSC. 

Note also that, since ~ £ P air&F and is directly related to the dimensionless ratio e F /\eb\, it 
can be inferred that £ pa ir&F is another physical parameter which can correctly determine the type of 
pairing. There is a very remarkable plot from [144], shown in Fig. 6, which examines the behaviour 
of the dimensionless chemical potential versus kpC, pair (\^b\ = e o i n the notations of [144]). This plot 
appears to be quite "universal", in the sense that it is remarkably independent of the specific model 
Hamiltonian and of the dimensionality (at least on the mean-field level) .\^\ The Fig. 6 also shows that 




Fig. 6. Chemical potential [i vs /cp^pair (at zero temperature) for "contact" potential and the dimensionality 
of space 2 < d < 3. Different curves are labeled by the values of d (in steps of 0.2). Positive values of [i are 
normalized by the Fermi energy e F = kp/2m, while negative values of \x are normalized by half the magnitude 
eo of the eigenvalue of the two-body problem in d dimensions. This plot was taken from [144]. 

the crossover between BCS and BEC regimes occurs in a rather narrow range of the parameter /cF^pair- 

Finally, note that the concentration dependence of the penetration depth was also studied in [141] 
and that the 2D crossover model (3.1) in the presence of a magnetic field was investigated in [157]. In 
particular the concentration dependence of the derivative dH&jcfT was studied and it was shown that 
[157] this derivative is substantially less in the Bose than in the BCS limit. 



Notwithstanding these similarities, one should be careful with the Bose limit for the discrete Hubbard model. 
This limit, as was firstly pointed out in [79], is quite different from that of continuum model (see Sec. 3.3). 



3.2 BCS-Bose crossover in the multi-band model: the coexistence of local and Cooper pairs 

Clearly the one band model as considered in Sec. 3.1 only describes the part of the phase diagram pre- 
sented in Fig. 1 corresponding to the superconducting state. Since we have considered zero-temperature 
we have in fact only described the single left-hand point with zero T c in Fig. 1. At this point T c = and 
it thus corresponds to the point A = e_p = in Fig. 5. We will return to possible theoretical explanations 
for the experimental phase diagram Fig. 1 when we consider the generalization of the one-band model 
to finite temperature. 

There is, however, an important question related to the peculiarities of the crystalline structure of 
the cuprates. It is known that the La 2 _ x Sr x Cu04 system discovered by Bednorz and Miiller [1] has 
T c PS 40K at optimal doping and contains a single Cu02 layer in a cell. On the other hand the later 
HTSC compound YBa2Cu306+«5 [158] has T c ps 90K and two cuprate layers per unit cell. The entire 
homologic family of cuprates has the composition A-M-Ca-Cu-O, where A = Bi, Tl, Hg; M = Ba, Sr, 
and the number of Cu0 2 layers per unit cell, 1 << N p[ < 6. Since these cuprate layers are situated 
in the same unit cell the arguments from Sec. 2.2 about the independence of the different Cu0 2 layers 
are not applicable in this case (see, for example, a recent paper [159]). Thus the possibility of coherent 
tunnelling and pairing between such layers must be taken into account. Taken together with the specific 
features of the BCS-Bose crossover problem discussed in Sec. 3.1, the effect of tunnelling and pairing 
between adjacent layers in the same unit cell produces rather interesting physics which we will briefly 
discuss in this section. 

Many-band models of superconductivity have been extensively studied in the context of HTSC [162] (see 
also the review [163]). Here we consider one such model [150] which is directly related to the geometrical 
structure of cuprates, and with particular emphasis placed on the BCS-Bose crossover. 

3.2.1 Model Hamiltonian 

Of course, a model which completely reflects the peculiarities of the crystalline and electronic structure of 
such HTSC materials would be too complicated for consistent calculations. For this reason the following 
simplifying assumptions were made in [150]: the HTSC material under investigation consists of identical 
metal blocks separated by spacers; superconductivity emerges in each iV^-layered metal independently 
so that all the actual interactions are confined to a single block; inside this block, the coupling differs 
from zero only for the nearest adjacent layers. 

Ultimately, we arrive at the Hamiltonian of a multilayered conductor possessing all the above mentioned 
characteristics (see also Fig. 7): 



H — H Q + Hy 



(3.38) 



where 
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Fig. 7. The schematic representation of the interactions included in the Hamiltonian (3.38). 



is the Hamiltonian of free charge carriers and 



Hint = "^EE // dh.dh^l^t^l^t) X 

xU jlj2 (r 1 - r 2 )ip j2 a{r 2 ,t)^ jl(7 {r 1 ,t) (3.40) 

is their interaction. 

In Eqs. (3.39), (3.40), the following notation is used: 4>} a (r, t), ipj a (r,t) are the Fermi operators of a 
particle with an effective mass m, spin <j(= —a) and r = (x,y) from the j-th (1 < j ' < N p i) layer 
satisfying the boundary condition ipo a {r, t) = il>N pl+ i t<7 (r, t); n is the chemical potential fixing the carrier 
density Nf(= JV p jn/); tjj+i(ri— r 2 ) is a j-independent one-particle interlayer tunnelling; [/^(ri— r 2 ) are 
the parameters of the intraplanar (ji = j 2 ) and interplanar (ji 7^ j 2 ) where the interaction depends only 
on \j 2 Here the sign of the interactions has been chosen such that Uj^ 2 {ri — r 2 ) > corresponds 

to attraction between carriers with opposite spins and we used real time t. 



The component of the Hamiltonian (3.38) corresponding to the free charge carriers permits diagonal- 
ization if we introduce new Fermi fields in accordance with the representation 



Wja{r, t) = 2^ u jv <p u<T {T, t) ; u jv = sin ■ 
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E U 3 V \ U 3 V 2 ~ J E U h vU h v 3jlj2- (3-41) 



As a result, the free Hamiltonian is transformed to the following diagonal form: 



pl r ( V 2 \ 

which describes an iV^-band metal in the effective mass approximation. In this case, a carrier belonging 
to the u-th band is characterized by the effective mass m u = (m _1 + Am -1 cos iris/ (N pl + 1)) and by 
the chemical potential 

\l v = fi + t cS cos nU . (3.43) 
IVpi + t 

The renormalization of mass and the position of each band on the energy scale are determined by the 
constant t, J+ i(r) of interlayer hopping since £ e ff = 2i,j + i(0), and Am -1 = a 2 V 2 ijj+i(r) | r=0 (a is a 
parameter having the dimensions of length and associated with the bandwidths in this approximation 
through the relation W v = 2n 2 / m v a 2 ) . The chemical potential (3.43) allows us to judge (see Sec. 3.1.3) 
what type of pairing (i.e. local or Cooper) occurs in the i^-th band simply by looking at the sign of the 
corresponding \i v . 

In the new variables (3.41) and for local interactions Uj 1 j 2 (r 1 — r 2 ) = Uj 1 j 2 (0)5(r 1 — r 2 ), the interaction 
(3.40) assumes the form 

1 Npl r 

Hrnt = —z E E / ^tfr.^lMl^Ml^M), (3.44) 

where the matrix element U\ v can be expressed directly in terms of the initial parameters Uj 1 j 2 (0). 
For the sake of simplicity we will retain only the following two values (see above): U intra = Ujj(0) and 
U inter = Ujj + i(0). As a result, expressions (3.42) and (3.44) can be regarded as a generalization of the 
well-known two-band model of superconductivity [160,161] to the case of an arbitrary number of bands. 
The symmetry of the chosen model (the equivalence of all planes) and the nature of intra- and interlayer 
interactions are evident. As an example Eq. (3.44) contains no terms corresponding to interband pairing 
which generally appear in phenomenological models [162,163]. Naturally, this does not mean that each 
band in the model (3.42) and (3.44) behaves independently because the rearrangement of the vacuum 
(the emergence of anomalous mean values or the condensate) for charge carriers for one of the bands 
immediately leads to the same rearrangement for charge carriers from the other bands. 

As noted previously many-band models of superconductivity have been studied in the context of HTSC 
[162] (see also the review [163]). The new idea that is brought to the many-band model by using the 
BCS-Bose crossover formalism is in establishing the correspondence between the sign of the chemical 
potential \i v of the v-Wi band and the nature of pairs. 

The further analysis of the many-band Hamiltonian (3.42), (3.44) is very similar to that discussed in 
Sec. 3.1 and we need not repeat it. Omitting the general equations derived for the many-band case [150] 
we show only the equations for the two-band model and their numerical solution. 

3.2.2 Bilayered cuprates 

It was mentioned above that bilayered materials include Y-Ba-Cu-0 compound. The chemical potentials 
for the two bands are given by /i u — [/,— (— l) v t e s/2. It proves convenient to introduce the dimensionless 
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giving the couplings of the superconducting order parameters from the same and different bands re- 
spectively. One may then obtain (see [150]) the following system of equations for the superconducting 
gaps Ai, A 2 and for the chemical potential /i: 
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(3.46) 



If one assumes that intra- and interplanar interactions are the same, i.e. Ui ntra = Uinter the system (3.46) 
can even be investigated analytically [150] leading to a better understanding of the underlying physics. 
Here however we present only its numerical solution shown in Fig. 8. To avoid confusion with the two 
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gaps A pair and A coh discussed above, we would like to stress that both the gaps Ai and A 2 in Fig. 8 are 
single-particle excitation gaps and have the same origin as A pair . The observation of different gaps in 
a similar context has been discussed in the literature [164,165]. It would be interesting to analyze the 
modern experimental data, for example, [51,52,58,155] bearing in mind that Y-Ba-Cu-0 cuprate may 
have in fact two single-particle excitation gaps. 

The behaviour of the chemical potentials Ati,2 shown in Fig. 8 demonstrates a very interesting feature of 
the many-band model. While A-fi > for almost the whole range of carrier densities which correspond 
to the Cooper pair regime in this band, one finds that fi2 < for a reasonably large range of carrier 
densities. This means that the pairs in this band are local and one has a coexistence of local and Cooper 



pairs in the same system. Of course, for higher carrier densities both [i\ and ^ become positive and 
the BCS regime described by the two-band model [160,161] is restored. 

In conclusion it is notable that the idea of mixed BCS vs BEC behaviour in a two gap model of 
superconductivity was recently readdressed in [166]. In this paper fermions from different parts of the 
Fermi surface with small and large Fermi velocities respectively were treated as belonging to two different 
bands. 



3.3 Peculiarities of the s-wave crossover on the lattice 



The physics of the BEC - BCS crossover can be well understood on the basis of continuum models. 
Nonetheless real superconductors are crystals and, if the pair size is not much larger than the lattice 
spacing (as in the case of HTSC), lattice effects are important and should be considered. 

Let us review the crossover at T = for the lattice version of the continuum model (3.1), given by 

H = -t ( c L c ma + h.c.) n n] n nl - /j, ^ n n(7 , (3.47) 

(nm),<r n n,cr 

where t is the transfer integral between neighbouring lattice sites n = (n x ,n y ), m = (m x ,m y ); c na ,c na 
are electronic creation, annihilation operators for the site n and spin a, respectively, U > is the strength 
of the on-site attractive interaction between two electrons occupying the same lattice site, n ncr = c na c na , 
and Ho is the chemical potential. This Hamiltonian is also called the attractive or negative-U Hubbard 
model. The word "negative" indicates that there is a negative sign before U in contrast to the positive 
sign in the original Hubbard model. 

Amongst others, this lattice model defined by (3.47) or (3.49) has been studied in [92,95,99,105,106,116] 
for T 7^ and in the many other papers mentioned in Sec. 2.3. 

The number of charge carriers per lattice site reads 

n = n(fi,T) = -^J^Jw) 

1 ' n,cr 

where n G [0,2]. In momentum representation the Hamiltonian (3.47) can be written as follows 

H = ]T(e k - Vo)cl a c ka - ^ S c kT c k+qt c p|C P - qi , (3.49) 

Icct k,p,q 

Here we have used the notation k as an index, = — 2t(cos k x + cos k y ) is the simplest band energy in 
the nearest neighbour approximation, we have set the lattice constant a — 1, TV is the number of lattice 
sites, and cj Kr (ck (T ) is the creation (annihilation) operator for momentum k and spin a. 

We shall predominantly follow the work in [167] and write down the standard coupled equations for the 
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where here (compare with (3.15)) = y + A 2 and £k = s^—^i, where the effective chemical potential 
is /i = /io + Un/2 i.e. it contains the Hartree-Fock shift —nU/2. This shift is essential to ensure that 
the true chemical potential /io coincides with — |e&|/2 in the Bose limit. 

The crossover point itself is simply the point where \x = 0. However, as motivated at the end of Sec. 3.1, 
it proves physically meaningful to consider the Bose and BCS regions as defined in terms of the pair-size 
£ P air (see Eq. (3.35)) for opposite spin-fermions: the condition k F £ pair < 1/ir identifies the Bose region 
and /cF^pair > 27r the BCS region. From the solution of the equations one can then construct the "phase 
diagram" [167] (the bound state energy there is denoted as e ) shown in Fig. 9 (a) which should be 
compared to that for the continuum contact potential shown in Fig. 9 (b) considered in Sec. 3.1. 

We note that the "phase diagram" for the lattice case is symmetric about half filling n — 1. From the 
diagram it is clear that to obtain a crossover from the BE to BCS regions as one changes the density 
(at fixed interaction strength U/t) one needs U/t < 2.4. In this case the BE and crossover regions 
occur only at extremely low densities. This is in contrast to the continuum case where one can obtain a 
density-induced crossover for all coupling strengths. The reason for this difference is the reentrant shape 
of the curves in the lattice case resulting from the van Hove singularity in the density of states as one 
approaches half-filling. Note that in contrast for the continuum case the density of states is constant. 

One should also note that, in contrast to the continuum case, one cannot always reach the "dilute" 
boson limit on the lattice at high densities even in the limit of very strong coupling. This effect is not 
the result of the finite size of the composite bosons which may be regarded as point-like in the infinite 
coupling limit. Instead it is due to the "overlap" of the centers of mass of the composite bosons since 
near half-filling the distance between any two bosons is a minimum of a single lattice spacing. Thus one 
expects the fermionic degrees of freedom to again predominate near half-filling on the lattice [79]. 

This argument may be quantified [167] by considering the commutator of the following boson-like 
operator 

^E^W (3.52) 

k 

where represents the pair wave function. The commutator [b, W] may be regarded as a c-number 
provided that the occupation number of each relevant state (n^) < 1. This can be achieved at high 
density (i.e. a large total number of particles) only if there is an infinite number of k states available. 
This condition is clearly not satisfied on the lattice where only a finite number of k states exist. Thus, 
in the lattice case, one must satisfy both the condition /^pair < 1/vr and the condition n « 1 for the 
system to be regarded as a composite Bose gas. 



Eqs. (3.22) and (3.24) give these equations for the case of a quadratic dispersion law e(k) (see (3.15)) 




|U|/t 

(a) 




(b) 

Fig. 9. (a) - Phase diagram (U/t,n) for the s-wave solution of the negative-lf Hubbard model; (b) - Phase 
diagram (eo = [£{,[, k 2 F /2m) for the contact potential. These diagrams were taken from [167]. 



3.4 Crossover in the models with d-wave pairing 



In this Section we consider the problem of d— wave pairing in 2D at zero temperature both in lattice 
[167-171] and in continuum models [172,173]. The motivation for considering this type of pairing is 
the experimental observation in the ARPES and other measurements [4,8,13-15] of a d x 2_ y 2 symmetry 
in HTSC. We shall analyse the crossover behaviour both as a function of density at fixed interaction 
strength and as a function of interaction strength at fixed density. The interesting question as to which 
type of pairing symmetry is actually present as one changes the interaction parameters and the density 
has been considered in [168,169,171] but will not be addressed here. 



Firstly one should note that pairs of o?-wave symmetry cannot contract to point like bosons. This auto- 
matically implies that the interaction between composite bosons which results from the Pauli exclusion 
principle is of finite range. This in turn leads to a larger range of correlations between these bosons 
which has a dramatic impact on the crossover. In particular at moderately high densities and large 
couplings these correlations suppress the Bose degrees of freedom and give rise to a larger (fermionic) 
BCS region [170]. Secondly for the simplest model which permits c/-wave pairing there exists a critical 
interaction strength for the formation of a two-body bound state on the empty lattice while this thresh- 
old is zero for s-wave pairing. This eliminates the possibility for a density induced BCS-BE crossover in 
this model. Thirdly the <i-wave symmetry has important implications for both the excitation spectrum 
and the momentum distribution since at least for the chemical potential greater than the critical (or 
crossover) value the excitation spectrum is gapless in certain directions. 

The lattice models that have been considered are usually 2D Hubbard models. To obtain a c?-wave 
solution the fermionic potential must contain an inter-site term of strength V in addition to the on-site 
term of strength U considered in the previous section. The Hamiltonian is then given by 

H = ~t ( C L C m CT + h.C.) ~UJ2 n ^n nl - /In n na ~ V C nT C m| C m|C nT , (3.53) 

{ma), a n n .°" <nm> 

where V > corresponds to attraction and the notation (nm) denotes as above nearest neighbour 
pairs. The case of V = and U > simply gives the Hubbard Hamiltonian (3.47) which only has an 
s-wave solution. When V > (inter-site attraction) and U < (on-site repulsion) one can also find a 
d— wave solution of the type 

A(k) = Arf(cos k x — cos k y ) . (3.54) 



We first review the results for this model, the so-called t — V model, with only nearest neighbour hopping 
and for which the dispersion relation takes the form 

£k = — 2t(cos k x + cos k y ) — /j (3.55) 

where /x is the effective chemical potential 

H = no + n(U + AV)/2 (3.56) 

which, as in the case of s-wave pairing, incorporates the Hartree shift. 

In this case there is a critical coupling strength V cr /4t = 1.83 below which no two-body ground state on 
the empty lattice exists and which gives the start of the BE region at zero density. This is in contrast to 
the s-wave case which has no threshold for pair formation. For a many-body s-wave Cooper instability 
to occur in 2D one needs the existence of a s-wave bound state for the two body problem on the empty 
lattice [53] J_] However no such condition applies to higher angular momenta and one can therefore 
obtain a BCS ground state for couplings less than V cr . There is however some controversy as to the 
minimal coupling strength actually required for such a paired ground state. In [167] it is claimed that 
pairing always occurs, or equivalently one always has a non-zero value of A^. As stated in this work 
[167] this is in direct contradiction to the earlier work of [170] which finds an unpaired region for very 
weak coupling. In our opinion, the result of [167] is more plausible since the presence of the Fermi 

16 We note that the two-body problem on the empty square lattice was solved many years ago without reference 
to superconductivity [174] and more fully investigated for the t — J model by Kagan and Rice [175]. 



surface reduces the dimensionality of the problem [50], and therefore favours the formation of a paired 
many-body ground state. 




(a) 




|V|/4t' 



(b) 

Fig. 10. Phase diagram (V/4t = \V\/4t,n) for the d-wave solution of the extended Hubbard model with 
attraction between nearest-neighbour sites in two dimensions, considering (a) - nearest-neighbour or (b) - 
second- and third-neighbour hopping. These diagrams were taken from [167]. 



For coupling strengths less than V CT no density induced BCS-Bose crossover is possible. For larger 
couplings strengths, V > V CT , it is clear from the phase diagram in [170] that one may achieve the 
crossover point for which /j, crosses the bottom of the band by increasing the density. One cannot 
however reach a true BCS phase for which /cF^pair > 27T by varying density as can be seen in Fig. 10 
(a). Thus independent of the coupling strength and in contrast to the s-wave case, a density induced 
crossover is not possible. This is a result of the presence of a pairing threshold [167] as evidenced by 
considering the t' — t" — V model, or extended attractive Hubbard model [103], with second and third 



neighbour hopping. For this model the dispersion relation is 

£k = At' cos k x cos k y + 2t" (cos 2k x + cos 2k y ) — \i (3.57) 

with t' > and t" > and where \x is again the effective chemical potential given by (3.56) .[^] In this 
case the critical value V cr vanishes for t" < 0.5t' and a density induced crossover is indeed possible as 
seen in Fig. 10 (b). 

For moderate densities, one becomes aware not only of the finite lattice spacing, but also of the intrinsic 
finite size of d-wave bosons. As stated above this finite size enhances the interactions between composite 
bosons leading to an enhanced fermionic region. For an s-wave system one can obtain a crossover from 
fermionic superconductivity to bosonic degrees of freedom with increasing interaction strength for any 
choice of the density. For a (i-wave system such a crossover is only possible in the dilute regime. 

One can describe the crossover in terms of the effective chemical potential which is plotted for the 
t — V model in [170] as a function of the coupling strength for fixed density. Here one confirms that 
at large doping the chemical potential differs little from its normal state value for a large range of the 
coupling strength in contrast to both the low-density and the s-wave results. The limiting case here 
is half-filling n — 1 for which the effective chemical potential is zero, and thus the system is always 
a fermionic superconductor. It is clear from this work [170] that one can access the crossover point 
at which the effective chemical potential reaches the bottom of the tight-binding band, albeit at ever 
increasing coupling strengths with increasing densities, for all values of the density below half-filling. 
One can never however reach the Bose-limit for higher densities (n > 0.016) for both this model and 
the t' — t" — V model [167] since for these densities /impair > V 71 " even for infinite coupling. This is 
apparent from the line labelled BE in the phase diagrams shown in Fig. 10. This should be contrasted 
to the corresponding line for the s-wave case. It is also confirmed for the t — V model by a plot of the 
pair-size in the infinite coupling limit as a function of density [170]. From such a plot one can see that 
the pair-size in the infinite coupling limit even diverges as one approaches half-filling. 

The chemical potential also highlights the interactions between the composite bosons in the BE region 
(n < 0.016). In the s-wave case the chemical potential approaches minus one half the two-body binding 
energy in the limit of infinite coupling corresponding to non-interacting bosons. This is only true for d- 
wave symmetry if the density is zero [170,167]. Defining the boson chemical potential by Hb — 2/Uo + | £ &| 
where is the two-body binding energy its behaviour in the BE region is given by 

Hb^-(U + 2V)n (3.58) 

One thus sees that one has /ie < for < —U < 2V corresponding to an average inter-boson attraction 
and fx b > for < 2V < —U corresponding to an inter-boson repulsion which increases with increasing 
interaction strength V (and —U). 

Returning to the pair-size it can be shown that it does not decrease monotonically with increasing V 
but converges asymptotically to a value larger then the lattice spacing for all n < 1. The asymptotic 
value increases with n and diverges as one approaches half-filling. This unusual behaviour is absent in 
the s-wave case and is related to quasi-long-range-order correlations among the bosons which reside 
individually on nearest neighbour sites. The pair-size in this case measures the correlation between 
opposite spin fermions belonging to different composite bosons. This can most clearly be seen in the 



This dispersion law captures the small-arc features of the Fermi surface detected in underdoped cuprates. 



lattice Fourier transform of the pair-wavefunction A^/E^ at half-filling. In the very large coupling limit 
V/t » 1 and at half-filling 



— = - sgn (cos k x - cos k y ) (3.59) 

with lattice Fourier transform 

0(Rn) = -J 1 2 J 2 ( 3 - 6 °) 

7T 2 n 2 - n 2 y 

where R n = (n x , n y ). It decays as a power law with increasing density when n x +n y is an odd integer and 
vanishes when n x + n y is an even integer (including the case n x = n y = for same-site correlation). In 
contrast in the s-wave case the pairing wavefunction tends to a constant with lattice Fourier transform 
Ar„,o- The divergence of £ pa i r is related to the strong interaction V/4t — > oo between fermions of opposite 
spin but on neighbouring lattice sites. It is not related to the divergence of the size of the composite 
bosons which is related to correlations between fermions of opposite spins but on the same site. For this 
reason the additional area in the top left of the phase diagrams (Fig. 10) refers to "correlated bosons" 
(CB) rather than a BCS state even though k F ^ pair > 2n. 

It is interesting to note that the phase diagram shown in Fig. 10 resembles the phase diagram obtained 
in [177] (see also the review [87]) for the same model, yet taking into account antiferromagnetic ordering 
competing with superconducting ordering. 

ci-wave pairing has also been considered in the context of a continuum model [173] where the system 
must be regarded as dilute for the approximations used. For low densities the results are similar to 
those for the t — V model [170] including the existence of a coupling threshold V CI . 

Lastly the ci-wave symmetry has an important effect on the excitation spectrum and the momentum 
distribution [173,170,176]. The gap parameter is zero at azimuthal angles = ±7r/4, ±37r/4. Thus for 
the continuum model with \i > 0, or the lattice model with \i above the bottom of the band, the quasi- 
particle spectrum is gapless at the four Dirac points given by = ±7r/4, ±37r/4 and k — k^ defined by 
£k M = 0. Near these points the excitation spectrum is linear. Now consider the direction = (or the 
equivalent directions = ±7r/2,7r) and let E g {k ll ) denote the energy gap at the point k = k^ in this 
direction. The E g {k ll ) is a non-monotonic function of fi reaching a maximum at intermediate values of 
/i above the crossover point. As fi reaches the crossover value, k^ becomes zero. The minimal value of 
the energy gap is thus at the single point k = and is given by -5^(0) = |A(0)| = implying that 
the spectrum disperses quadratically for small momenta for all 0. As soon as /i < for the continuum 
model, or /i (which is physically the same) is below the band edge for the lattice model, the spectrum 
becomes fully gapped but the minimal gap remains at k = and takes the value of \fj,\. This behaviour 
is strikingly different to that for the s-wave superconductor which is always gapped. 

Furthermore this behaviour of the excitation spectrum implies that the single particle distribution 
strength changes dramatically as one passes through the crossover value for \x for ci-wave pairing 
[173,170,176]. This discontinuity can be clearly seen in the momentum distribution at k = which drops 
to zero discontinuously as the chemical potential crosses the critical value [173]. No such discontinuity 
is seen in the s-wave case. It is also illustrated by the contour plots in [170]). For weak coupling, the 
chemical potential is essentially equal to its normal state value, and one has a slightly smeared Fermi- 
surface. The momentum distribution is thus maximal at the centre of the Brillouin zone (0, 0) and 
falls off least rapidly in the anti-nodal directions = 0, ±|, it. For strong coupling, when the chemical 



potential has crossed the bottom of the band, the momentum distribution is zero at the centre of the 
Brillouin zone and has a lobe-like structure centered around (±7r, 0) and (0, ±7r). These directions with 
high pair occupation probability may well be the preferred direction for the opening of the pseudogap 
[170]. 

This qualitative changes in the excitation spectrum and momentum distribution substantially affect the 
quasi-particle density of states [167]. For low frequencies this is always zero in the s-wave case because 
the system is gapped. For the <i-wave case it changes discontinuously with ji from linear in uj above 
the crossover value, to a constant at the crossover value, to zero below the crossover value where the 
spectrum is again gapped. 



3. 5 Peculiarities of the crossover in models with retarded attraction 

As seen in Sec. 2.2 Uemura suggested [84] (see Fig. 4) that the attractive interaction between carriers 
leading to superconductivity is non-retarded on the Bose side of the BCS-Bose crossover, while on the 
BCS side the interaction becomes retarded. We could not consider this feature of the Bose-BCS crossover 
in the previous sections since all the models discussed were based on postulated predominantly local (as 
in Sees. 3.1 and 3.2) direct non-retarded attraction. It should be admitted that this type of interaction is 
not realistic and must be regarded as a consequence of a more fundamental interaction whose nature has 
not yet been established for HTSC It cannot be ruled out that, as in traditional superconductors, this 
is the electron-phonon interaction [86] or the electron-magnon (or electron spin-fluctuation) interaction 
[24] in view of the peculiar magnetic properties of HTSC. In any case in any model describing HTSC 
correctly, the attraction between charge carriers must be derived from the exchange of excitations of a 
bosonic field. In such a model, we can not only obtain the attraction between charge carriers, but also 
estimate the role of such an important factor as the retardation of the interaction. The corresponding 
generalization introduces one more parameter into the system, viz., the characteristic frequency uj 
of elementary bosonic excitations. In low-temperature superconductors, uo is the Debye frequency for 
which the inequality ujq <C tp is always satisfied. However, in the case of HTSC with a low number 
density of charge carriers, the situation with u < ep is more likely if, for example, optical phonons are 
responsible for the interaction. Alternatively we can assume that u pa ep if the interaction is realized 
through magnons; it might even be the case that u > ep (see e.g. [22]). 

The inclusion of non-retarded interactions into the Eliashberg equations [182] (see also the textbooks 
[49,181,148] and the review [183]) has been extensively studied and we will give the references shortly. 
However here we focus on the simplest model with an electron-phonon interactionT 8 ] but with a variable 
carrier density which makes possible to obtain a BCS-Bose crossover as considered in [178] (see also the 
more formal treatment in [179]). This is to the best of our knowledge the only attempt to consider the 
BCS-Bose crossover itself using a model with indirect interaction and we will briefly present the model 
and discuss the main results. Another related and very interesting question which we discuss here is a 
possible violation of Migdal's theorem [180] in HTSC. 



It was noted above that elementary Bose excitations of different kinds may play the role of an intermediate 
boson, but we will simply refer to them as "phonons" . 



3.5.1 Model and basic equations 



It is known that (see [152,181]) the Hamiltonian density describing a system with electron-phonon 
interaction can be written in the form 

U = -^t( r ) [T- + iA ^( r ) + H ph {v{r)) + g el ^l{v)^{r)y{r), (3.61) 

where 7i p /i is the operator describing phonons with the dispersion u;(k); r is a 2D vector; ^v(r), (p (r) 
are the Fermi and Bose fields; a is the spin; g e i- P h is the constant of electron-phonon interaction. 

In order to clarify the influence of retardation effects on the electron spectrum as a function of rif, one 
must analyze a system of two equations. The first equation is the Eliashberg equation [182] (see also the 
textbooks [49,181,148] and the review [183]) for the self-energy of the electron, S(p) written in Nambu 
formalism [147] as 

= ig c \-ph J j!^T 3 G(p')T 3 T cl - ph (p,p-,p- p)D ph (p - p'), (3.62) 

where p = (po, p) (p being a two-dimensional vector) and one integrates over both the frequency and 
the momentum. r e i_ p h(p', P', p — p') is the vertex function and D p ^(p — p') is the "dressed" phonon 
propagator. 

The second equation is the well-known condition imposed on the Green's electron function G(p)[^ 
relating the chemical potential fi to nj = eprnj-n: 

/d^p 
exp(i5p T 3 )tT[T 3 G{po, p)]^— ^ = n f> 8 -> 0+ • ( 3 -63) 



Obviously, it is impossible to solve the system of equations (3.62), (3.63) in such a general form, and 
one must make certain simplifying assumptions. 

One such assumption, which is often made in an analysis of the Eliashberg equation, is to disregard 
effects that are not directly connected to the emergence of anomalous mean values. 

The second assumption follows from Migdal's theorem [180] (see also the textbooks [152,184,148]) stating 
that for ojq <C ep, the electron-phonon vertex function can be replaced by its "bare" value g e i- p h- Note 
that one cannot use directly the result obtained by Migdal since the standard inequality u ep can 
be violated in HTSC as mentioned above. However, Migdal's theorem was recently generalized to the 
case when u>o ^> ep [185]. Hence the vertex function T in (3.62) can also be replaced by g e i- P h in the 
case of sufficiently low densities (Bose limit) that c^o ^ Since the results will focus on the Bose and 
BCS limits, where Migdal's theorem appears to be valid, one assumes that the "bare" vertex can be 
used throughout. 

It is known [186,187], however, that the complete vertex should be taken into account for ujq ~ ep and 
this consideration leads to a significant increase in T c . The problem of solving the Eliashberg equation 
with a vertex correction [187] is very complicated and in fact deserves its own review. Nevertheless, 
for the limiting cases u -C ep and u ^> ep considered here, one does not need to include the vertex 



Note that Eq. (3.63) is exactly equivalent to Eq. (3.23) which was written in terms of the effective potential. 



correction. Nonetheless the final result obtained shows that, if one wants to explain why the BCS regime 
is not the most optimal for high T c , one may well need vertex corrections. The theory of the Eliashberg 
equation with vertex corrections is currently a very rapidly developing branch of superconductivity and 
we simply mention here some of the latest papers [188-191] where the relevance of this theory to HTSC 
and new theoretical developments were discussed. 

Furthermore, there is an interesting correlation between the need to include the vertex corrections into 
the Eliashberg theory and obtaining the correct limiting behaviour in 2D in the study of superconducting 
fluctuations. For example, it is known that the T-matrix approximation which is used extensively (see 
Sees. 2.3, 4.5) fails to describe BKT physics [95,135,92,102]. At precisely the same time vertex corrections 
are also crucial. 

Finally, the third assumption is associated with the choice of the form of the "dressed" boson propagator. 
As usual (see [49]), it is replaced by the "bare" propagator [152], 

^ w = a *^°' k ' ^ 0+ ' (3 ' 64) 



In the subsequent analysis, the Einstein model has been used for the phonon dispersion uj(k), i.e. 
cu(k) = ujq. Using this approximation and taking the function E(p) to be proportional to the order 
parameter one can transform Eq. (3.62) into an equation for the order parameter $ (see the 

details in [179]), which depends on the frequency p : 

S>( v ) = -iJ f d3p> (3 65) 

{P0) " 9d ~ ph J (2vr)3 (p' o y - £(p') 2 - |<&(p&)|2 + i6 (p - Po) 2 ~ ^ + i6 ' 

where again p = {p , p}- The problem is thus reduced to an analysis of the system of equations (3.65) and 
(3.63) as a function of the carrier concentration. It should be emphasized that the Eliashberg equation 
3.65 had been investigated and solved several times previously (see e.g. Refs. [49,152,181,182,192]) but 
only when the value of fi coincided with ep (i.e. one did not consider the BCS-Bose crossover problem). 
In this latter case the expansion £(p) Pf(p — Pf)/™, (pf is the Fermi momentum) was then used for 
the normal state spectrum. 



3.5.2 Electron spectrum 

The desired concentration dependence of the electron spectrum can be characterized by the quantity 
A(n/) = |$(po = 0, ri/) | which defines the energy gap. The system of equations (3.65), (3.63) was 
studied in [178] using the rather crude approximation that $(po) — const. It was also studied in [179] 
by applying a method (see Refs. in [179]) which allows one to reduce the integral Schwinger-Dyson 
equation (3.65) to a differential one, which can then be solved relatively easily. This more accurate 
treatment showed good agreement with [178]. Thus here we simply present the main classes of solutions 
of the system derived in [178]. We assume everywhere that the conduction bandwidth W 3> uj ,e F , A 
as is appropriate for HTSC. 



3.5.2.1 BCS or high density limit In the BCS limit, the inequality ep ^> u>o,A holds. In this 
case, it follows from (3.63) that ji ~ ep. 



a) If we assume that additionally A <C cj , we reproduce by solving (3.65) the following formula from 
BCS theory [152]: 



A = Abcs = 2u exp ( ^— ] , (3.66) 

Thus the above inequality is, as expected, satisfied if the coupling is weak (g^_ ph m/2-K <C 1). Note 
that Eq. (3.66) corresponds to the result of the classical Morel-Anderson [192] model (if we adopt the 
notation mg^_ ph / / 27r = A). The quantity A = $(0) was obtained in [192] by substitution of a suitably 
constructed trial function $(po) 7^ const into an equation similar to (3.65) and this represents the first 
known way of approximating the analytical solution of the Eliashberg equation. 



b) If one assumes the opposite limit, A ^> u , we obtain 

A = u 



This result is valid if we assume that the coupling is strong (g^ h m/16 3> 1). If A is much greater 
than ujq, it can be verified that the formula (3.67) corresponds to ultrastrong coupling [193] (see also 
[181]). Note that to be consistent in this limit the equation for the renormalization factor Z(p) must 
also be studied. 



3.5.2.2 Low-density limit Small values of rif correspond to the condition 6f, A <C loq. Since the 
value A in this case can be of the order of €f, Eq. (3.63) becomes nontrivial (see the discussion in 
Sec. 3.1.3) and does not lead to // = €f- 

a) If one also assumes that <C cj , one arrives at the following solution 

A = 2JuJtfiF~ exp [ ^ — ), n = -co exp( ^ — ) + e F - (3.68) 

V m 9li- ph J V m 9el- P hJ 

It can be seen that the condition ujq holds for g^_ ph m/4:7i <^ 1 i.e. weak coupling. Comparing 

(3.26) and (3.27) with (3.68), one can easily see that the latter equation can be written in the form 

, I ^.weak I 

A = \J2\e™ cak \eF, fi = - l -^ + e F , (3.69) 

where 

£ wcak = _ 2uJ(j exp ( (3>70) 



Since the expressions (3.70) and (3.26) are similar, it is natural to assume that er ak has the same 
meaning as defined by (3.26), i.e. is the energy of the bound state which emerges in this case due to the 
electron-phonon interaction. For this reason, expression (3.70) contains the quantity ooq characterizing 
the energy of phonons instead of the bandwidth W . 

At the same time, solution (3.69) differs significantly from (3.27). The latter expression is valid for all 
values of €f and all interaction strengths. The former applies only for weak coupling and low densities 



e F <C ojq- If the inequality e F ^> c^o holds, expression (3.66) must be used. For e F <C the charge 
carriers form an "adiabatically" slow subsystem, and the retardation of indirect interaction can be 
disregarded as is shown in Fig. 4. However, as the value ep increases (and especially in the region 
e F ^> u ), the role of retardation becomes decisive as is evident in formula (3.66). 

The numerical solution of (3.65) and (3.63) is shown in Fig. 11. It can be seen that the asymptotic 

0-5 i 1 — 1 




(3.69) for A and the value of Abcs determined from (3.66). 

expression (3.69) for A coincides with the numerical solution for e F <C u>o, while for e F ~ c^o and 
e F it lies above the corresponding numerical curve (Expression (3.68) can increase without bound, 

while if one incorporates the band structure then A achieves its maximal value at half-filling, see (3.25) 
and Fig. 5). The maximal value attained is however much lower in the system with an electron-phonon 
interaction. The gap cannot increase without limit as suggested by Eq. (3.68) but attains a saturation 
value equal to Abcs in the limit e F » u>o- 

The chemical potential fi (which is not shown on the figure) depends linearly on e F and is negative for 
very low concentrations of the charge carriers. It can also be seen from (3.69) that for e F < |^ eak |/2, 
the chemical potential \x < 0, and hence the local pair mode is realized (see Sec. 3.1.3). It is important 
to note that although the quantity was an independent parameter in Sec. 3.1 , in the case under 
investigation e^ cak is expressed in terms of electron-phonon interaction parameters. It can be seen from 
(3.66) and (3.70) that A BC s and ej eak are connected by the relation |e^ cak | = A BC s exp(-47r/m^_ ph ), 
i.e. ej eak is exponentially smaller than Abcs- Thus in the case of weak coupling, the region of existence 
of local pairs is very small and is only of theoretical interest. 

b) Let us now assume that although e F , A <C uj , the absolute value of ^> uj q (/i < 0). In this case 
one can obtain the following solution of the system (3.65) and (3.63): 



A = 2 



o; coth 3 \e F , 



(. 



47T 



/j, = —uj coth ( — ^ — ] + 2e F coth 2 



47T 



\mgl_ 



(3.71) 



It can be seen that the condition ^> u Q holds in the strong coupling limit (gl x _ ph m / An ^> 1). In this 
case, we can again introduce the energy 

£r " g s - 2 ^ C0 Mdfc) * (3J2) 

and write the expressions (3.71) in the form 

A = J2\ef ong \e F , /i = -^^ + 2e F coth 2 — ^— . (3.73) 

2 V m ^cl-ph/ 

Virtually all that has been said about the solution for (3.69) remains in force for (3.73). However, some 
differences do exist. In particular it can be seen from (3.72) and (3.73) that the region of existence of 
local pairs is substantially larger in the case of strong coupling. 

There is also some evidence that the strong coupling case may be relevant to HTSC due to the presence 
of an extended van Hove singularity. 



3.5.3 Discussion 

The results presented here correspond to an oversimplified model. Nonetheless they are a good illus- 
tration of the importance of the position of the Fermi level ep relative to the limiting boson frequency 
ujq, which is so transparently shown on Uemura's version of the HTSC phase diagram (see Fig. 4 where 
an arbitrary boson frequency oob replaces Uq). The relationship between carrier density and retardation 
is clear from the forms of the solutions both in the high density (see Eq. (3.66)) and low density (see 
Eqs. (3.69) and (3.73)) limits : 

!\j2\el h \ep, ep <^ u , (non-retarded), 
, 
\j2uj \el h \, €f > u , (retarded) 

where e\ is given by e™ cak from Eq. (3.70) or, for the non- retarded case only, by e s 6 trong from Eq. (3.72) 
respectively. 

However, the approximations above have raised certain questions. In particular the quantity A and 
hence T c attains its maximal value for e F 3> uj . The following question thus arises: why do HTSC then 
attain maximal values of T c at relatively small values of e F 7 Indeed if one assumes spatially homogeneous 
systems the presence of real local pairs as suggested by [51,52] even demands negative fi. This argument 
does not apply to spatially inhomogeneous systems. Moreover if one suggests that the optimally doped 
cuprates correspond to the BCS regime (see Sec. 3.1.5), the model above proposes an explanation as to 
why this limit gives the highest T c values. 

There also exists a possibility that one can exceed the value of Abcs by taking into account the vertex 
correction for the region e F ~ c^o [186-191]. 

We mention also other recent models [47,48] with an indirect (through spin fluctuations) interaction 
between carriers which have been studied recently and where the pseudogap behaviour is the result of a 
complex interplay between antiferromagnetic and rf-wave pairing fluctuations. It is very likely that the 



(3.74) 



ultimate explanation of the pseudogap behaviour is only possible within such models. The models with 
direct interactions between carriers, which we have predominantly discussed here, are only the very first 
step in understanding the unusual properties of HTSC. However, as we will see in the next chapter 
there are many unanswered questions even within the relatively simple attractive Hubbard model so 
often used for the study of crossover and pseudogap behaviour. 



4 Self-consistent T-matrix approximation and its limitations in the 2D limit 



In the previous chapter we considered a number of different models in the mean-field approximation 
which is of course justified at T = when continuous symmetry breaking is possible even in 2D 
[126]. Any generalization of these models to T ^ necessarily requires going beyond simple mean-field 
theory. Here we will describe the different forms of the T-matrix approximation widely used to study 
the Hamiltonian (3.47) or its continuum version (3.1) at nonzero temperatures following mostly the 
articles [55,93,95,99]. Unfortunately, as we will see in Sec. 4.5, even quite sophisticated approximations 
may not be sufficient for 2D systems. 



4-1 Self-consistent, conserving approximation 



We first discuss the T-matrix approach to the Hubbard model in the absence of a superconducting 
condensate i.e in the normal state. The T matrix consists of the sum of all particle-particle ladder 
diagrams (see Fig. 12) It is known that the T-matrix is the ladder approximation to the Bethe-Salpeter 




Fig. 12. Diagrammatic representation of Eq. (4.1) taken from [99]. The T-matrix contains the repeated scattering 
of two-particles. 



equation [184] 



T(ift m ,q) 



-U 



1 - Ux(itt m ,q) ' 

where x(^m, q) is the independent pairing susceptibility given by 

x(ito m , q) = G (ittni - iu n , q - k)G(iu n , k) 

' n, k 



(4.1) 



(4.2) 



Recall that u> n , and Q m are the fermionic (odd) and bosonic (even) Matsubara frequencies, respectively. 
This approximation is expected to be accurate in the low density limit when one can neglect particle- 
hole scattering and is valid when the product of the Fermi momentum, and the scattering length, 
a is small. Note that the function x(if2 m , K) is chosen to have different signs by different authors 



(compare, e.g. [95] and [99]), and the reader should take note of the consequence of this sign choice in 
the subsequent equations. 

The corresponding one-particle Green's function in (4.2) satisfies the Dyson equation (see the diagram 
in Fig. 13) 

G(iu n , k) = {G (iu n , k)- 1 - E(iu n , k))- 1 , (4.3) 

with 

G (icu n ,k) = — , (4.4) 

%UJ n - e k 

where T,(iu n , k) is the electron self-energy. 




Fig. 13. Diagram for the single-particle Green's function (solid line) in the self-consistent, conserving approxi- 
mation given by Eqs. (4.3) - (4.5). The thin solid lines represent the non-interacting Green's function Go- 

In the self-consistent T-matrix approach the self-energy is given by 

S(iw„, k) = T(i£l m , q)G(iVL m - iu n , q - k) (4.5) 

P-'* m,q 

which closes the system of self-consistent equations. 

In most cases the results for real frequencies are obtained from a numerical analytic continuation 
iuj n — > oj + iO [194] (see also the review [183]). Since this continuation is inaccurate (see an example in 
Sec. 7.1.3 and recent comprehensive analysis of the method in [195]) a real-time technique is also used 
to avoid this problem [99,100,104]. 

The above system of equations, valid in the normal state, is conserving, i.e. satisfying Kadanoff-Baym 
criteria [89] and "$ derivable" in the sense of Baym [88] and known to include Gaussian fluctuations of 
the pairing field in a self-consistent manner. "$ derivable" approximation means that the corresponding 
self-consistent approximation for the free energy (per lattice site) is given by the functional 

2 

Q[E, G] = -— exp(iuj n 5) {£(iu; n , k)G(iu n , k) 

' n,k 

+ ln[-G (^ n ,k)- 1 + S(^ n ,k)]} + <I>[G'], 5^0+, (4.6) 

with 

®\ G \ = \ E I 111 ! 1 " U X (in m , q)] + UxiiVm, q)} (4.7) 

P m,q 



and at the stationary point G, E and $ are related by Dyson's equation (4.3) and 



which reproduces Eq. (4.5). Using short-hand notation one may rewrite (4.6) as (compare with Eq. (3.8)) 

Q[T, /I, E, G] = -2Tr[SG + Ln(-G + £)] + $[G] , (4.9) 
where TrA = 1/{N(5) £ n>k 4(iw n ,k). 

The term "conserving approximation" implies that the approximation conserves the number of particles, 
the energy, the momentum, and the angular momentum. As proven by Baym [88] "$ derivable" theory 
is conserving [89], although this is not a necessary condition for the theory to be conserving. 

In fact the functional (4.6) (or its short-hand form (4.9)) is closely related to the Cornwall- Jackiw- 
Tomboulis effective potential [196] (see also the textbook [129]) which is widely used in quantum field 
theory for the description of systems of equations similar to those given here. For example, this method 
easily permits one to generalize the system (4.1) - (4.5) to a state with a broken continuous symmetry 
e.g the superconducting state in 3D. Note that such equations were obtained in [97] (see also [90,98,101] 
where these equations were also discussed and solved by different methods) for a 2D model using the 
functional derivative technique of [88,89] and solved approximately using a momentum approach. 

The electron density as a function of T and /i is given in terms of the thermodynamical potential by 

dQ 

n(jM,T) = — (4.10) 

If one calculates n(T, //) from fi[E,G] and exploits the stationary properties of this functional, one 
obtains immediately the result, 

n(n, T) = J- ]T G(iu n , k) exp(^ n <5), 5^0+, (4.11) 
which is equivalent to that given by (3.48) as expected. 

The superconducting instability of the normal state at T c is signalled by the divergence of the T-matrix 
at zero energy (at chemical potential) and zero momentum (the well-known Thouless criterion): 

1-1/X(q = 0,w = 0) = 0. (4.12) 

This equation, and the associated normal state properties, can be examined in a variety of approxima- 
tions. Approximations are often used to simplify the fully self-consistent, conserving theory described 
by Eqs. (4.1) - (4.5), (4.11) since they cannot be studied analytically without additional simplification. 
However they can be and have been studied numerically [96,99,101]. 

Many of the discussions related to the different predictions and scenarios of these theories, which 
are all based on superconducting precursor fluctuations, are related to the fact that these additional 
approximations lead to quite different results (see, for example, the different opinions about the most 
adequate approach expressed in [55,57] and [105,106]). It is thus very difficult to say, at this stage, 
whether a particular theory really predicts different physics or whether the differences are simply an 



artifact of the approximation used to solve the T-matrix equations.^] Our point of view (see Sec. 4.5) 
is, however, that even in its best form presented in this section, T-matrix approximation may not be 
sufficient for the description of vortex phase fluctuations which are important in the 2D limit at finite 
temperatures and are the main subject of this review. 



4-2 Non self- consistent, non-conserving approximation 



The self energy in this case is given by (compare with (4.5)) 

S (z^ n , k) = — T (iQ m , q)G (iQ m - iu n , q - k) . (4.13) 

The index "0" indicates the use of free Green's functions so that the T-matrix (4.1) is also calculated 
in terms of the free Green's functions 

T Q \iQ m , q) = -— + Xo{i£lm, q) 

= -— + G o(^m ~ iu n , q - k)G (iu n , k) . (4.14) 

U n,k 

The summation over Matsubara frequencies for this lowest order approximation can be easily performed 
and the T-matrix (4.14) is given in terms of the corresponding pair susceptibility (4.2) which acquires 
a rather simple form 

q) = ^E 1 - /( W) - > ,4.15) 

^ k 4k+q/2 + 4-k+q/2 ~ " 

with £k — ^k — A* an d f( x ) = ( e + i s the Fermi function. 

The full Green's function in this non self-consistent and non-conserving approximation is the expansion 
to the first order in the self-energy of (4.3): 

G(iw n , k) = G (iuj n , k) + G (iu n , k)S (^„, k)G (iu n , k) , (4.16) 

which means that G in (4.16) and E (4.13) are not calculated self-consistently. 

This system corresponds to the approximation used in the earlier work of Nozieres and Schmitt-Rink 
[79] (see also the reviews [53,54]) to analyze the BCS-Bose crossover in 3D theory, by solving the system 
of equation comprising the Thouless criterion (4.12) 

1-L/Xo(0,0) = (4.17) 



20 There are probably two criteria one can use as a guide to decide which of the existing approaches is the 
best. Firstly, a good approach should violate as few exact results as possible. Secondly, a good approach should 
reproduce the exact or numerically accurate results whenever these are known. Based on these criteria the best 
approach of the T-matrix type for the repulsive Hubbard model in the weak and intermediate coupling regimes 
was probably presented in [135]. We refer the reader to this paper where the generalization of these results to 
the negative-?/ Hubbard model is also discussed in detail. 



in terms of xo given by Eq. (4.15), and the number equation (4.11) (see below its final form given by 
Eq. (4.18)). In particular, to calculate the fluctuation correction in the number equation (4.11) they 
replaced Q[E, G] by J1[0, G ] before differentiating. In practice, Nozieres and Schmitt-Rink replace the 
frequency sum in Q by the standard contour integral representation [49,148,184], and deform the contour 
to lie along the real axis before differentiating with respect to fi. An alternative (but formally equivalent) 
procedure is to retain the explicit frequency sum, and to evaluate the derivative — <9$[Go]/<9/i by first 
differentiating $ with respect to Go using Eq. (4.8), and then differentiating Go with respect to /i. The 
result is 

2 

n(n, T) = n + — G (iu n , k)£ (^ n , k)G (iu n , k) , (4.18) 

P n,k 

which exactly corresponds to the use of the Green's function (4.16) defined above. In (4.18) the density 
n is expressed via non-interacting Green's function G . 

Note also that the thermodynamical potential 

Q L = $[G ] = W (M 1 - Uxo(M m , q)] + Uxo^m, q)} (4.19) 

P m,q 

which includes the contribution from the sum of all particle-particle ladder diagrams evaluated with 
non-interacting Green's functions was firstly studied by Thouless for separable electron-electron inter- 
action [197]. The equivalent representation of this contribution to the free energy is given by Gaussian 
fluctuations of the pair field in a functional integral representation of the partition function, as was 
pointed out by Langer [198] (see also the more modern consideration in [149]). 

An extension of the Nozieres and Schmitt-Rink approach [199] to the case of 2D continuum gas in 
the same approximation as described here suggested the breakdown of the Fermi-liquid picture for any 
attraction. Moreover, for any electron density the chemical potential goes to the energy of two particle 
bound states as T — > in evident contradiction with the results of the mean field theory presented in 
Sec. 3.1. This specific behaviour of \x has been associated in [199] with the existence of large q bound 
states in the 2D system mediated by s-wave attraction. This was, possibly, one of the earliest papers 
where the composite hard core bosons were replaced by less stable objects in an attempt to explain the 
observed properties of cuprate superconductors. 

However, the approximation scheme applied in [199] is not sufficient for 2D systems because it is not 
conserving [93,116], and neglects the interaction between fluctuations (see also Sec. 4.5). In particular, 
it was shown in [116] that the values of entropy, S and specific heat, G become negative in some 
temperature interval if the thermodynamical potential Q[0, Go] which leads to the non self-consistent 
Green's function (4.16) is used to calculate them. The Green's function (4.16) even leads to a negative 
spectral weight, A(u,Vl) = — lmG(u + iO, k). Indeed, taking G (iu n ,k) given by (4.4) and £ (iu; n , k) = 
const or simply a function that does not have a zero at iu n = s(k), one can see that one of the 
contributions to A{u, k) coming from the last term of Eq. (4.16) is 

-Imf -—] = Im^-( —) = n-^-6(u - e(k)) . (4.20) 

\uj + i0 -e{k)j du \u + iO - e(k) J du v K " K J 

The derivative of the delta function is negative in certain regions, as one can see by taking one of the 
many representations of the delta function as the limit of a regular function. In other words, all double 
poles in (Go(iu n , k)) 2 result in a negative contribution to the spectral weight. That is one of the reasons 



why Dyson's equation (4.3) is necessary: to sum high-order poles to infinite order and transform them 
into an eigenvalue shift. 

All of this represents a serious gaffe of non-conserving, non self-consistent approximation. 
4-3 Non self- consistent, conserving approximation 

The conserving character of the approximation described in Sec. 4.2 can be, however, restored if we do 
not expand the Green's function (4.3) in terms of £ as we did with Eq. (4.16). At the same time one 
retains the same approximations for the T matrix (4.14) and self-energy (4.13) so that the Thouless 
criterion is unchanged and given by Eq. (4.17). This conserving but still non self-consistent approxima- 
tion was proposed by Serene [93] and the Green's function in this case is shown diagrammatically in 
Fig. 14. As discussed in [93], such an approach entails taking the sum of all repeated scatterings of an 



+ 



Fig. 14. Diagram for the single-particle Green's function (solid line) in the non self-consistent, conserving 
approximation given by Eqs. (4.3) and (4.13). The thin solid lines represent the non-interacting Green's function 
Gq. This diagram has been taken from [99]. 

electron by independent pair fluctuations, but omitting the interaction between fluctuations and vertex 
corrections. There is a debate (see, e.g. [93,116,99]) as to whether the restoration of the conserving 
character of approximation will preserve the physics suggested in [199]. Finally we write the number 
equation (4.11) for the non self-consistent, conserving approximation 

n{n, T) = ]T exp(iu n 5) [G^\iuj n , k) - S (^ n , k)]- 1 , 5^0+. (4.21) 

It is clear that Eqs. (4.18) and (4.21) differ because the first of them is obtained in non-conserving 
approximation. The two approximations are equivalent when the self-energy corrections are small, but 
if the corrections are not small, Eq. (4.21) might yield a sensible result when the non-conserving ap- 
proximation does not. 

4-4 Pairing approximation 

An alternative approximation to that discussed in the previous section has been used in the set of 
papers by Levin et al. [55-57]. 

The feedback effect of the self-energy £ on the T-matrix has been included into the theory by using the 
following equation for the T matrix (compare with Eqs. (4.1) and (4.14)): 

T-\i{l m , q) = -I + _L £ G (iQ m - iu n , q - k)G(iu n , k) , (4.22) 

• n,k 




where the full Green's function G is given by the Dyson equation (4.3) and the self-energy is calculated 
using T-matrix (4.22) and the bare Green's function Go 

S(zcj n , k) = — T{^m, q)G (iQ m - iu n , q - k) . (4.23) 

Thus one can see that the feedback effect is introduced into the theory via the presence of only one full 
Green's function in Eq. (4.22). And as usually this system should be completed by the equation for the 
chemical potential (4.11). 

It is stated in [55] that the system (4.22), (4.3), (4.23) and (4.11) is conserving because it satisfies Baym- 
Kadanoff [89] criteria. It is clear, however, that the theory is not "^-derivable" and non self-consistent. 
Based on the early experimental data which showed the absence of the pseudogap in the optimally 
doped cuprates it was assumed in [55] that the correct theory should be consistent with BCS theory. 
Claiming that the fully self-consistent theory discussed Sec. 4.1 does not reproduce the BCS limit in the 
weak coupling regime (see, e.g. [90,96]) it was argued in [55] that the mode coupling approach is more 
appropriate, since it more closely reproduces the results of BCS theory when the coupling is small. More 
recently however it has become clear that pseudogap phenomena are also observable in the optimally 
doped region [31], but over a far smaller range of temperatures than for the underdoped cuprates. This 
implies that one needs reproduce the BCS limit only asymptotically (see Sec. 5.1.2), so that the original 
objection against the fully self-consistent theory [55] no longer seems valid. 

It was also shown in [55] that the mode coupling approximation leads to the desired gap-like self-energy 
E R (u;,k) ~ 1/(c<j + ek + iT ), where r is pair-breaking scattering (see Sec. 7.2). However it is stated 
in [106] that this result is in contradiction with the fully self-consistent theory. This clearly shows that 
more studies are necessary to determine which approximations are sufficiently good that the results 
obtained are not simply an artifact of the approximation. 

4-5 Limitations of the T-matrix approximation 

In the previous Section we have discussed the discrepancies between different forms of the T-matrix 
approximation when it is applied to predominantly 3D models. We note also that this approximation 
neglects the particle-hole channel important for high densities. In 2D there is, however, a more funda- 
mental problem with the T-matrix approximation used in any, including the fully self-consistent, form. 
The essence of the problem is rather well-known: since the superconducting fluctuations are treated in 
an RPA-like or Gaussian approximation the BKT transition is not recovered by a T-matrix approach 
[95,92,135] (see also [200,201], where the effect of interaction between fluctuations was considered). 
One may hope however to derive BKT physics, or to be more exact the properties of the system for 
T < TbktPH only going beyond the T-matrix approximation, e.g. including the equation for the vertex. 
The T-matrix approximation in fact represents only the lowest order approximation to the infinite chain 
of Schwinger-Dyson equations for the fermion and boson Green's functions where these have been trun- 
cated by neglecting the equation for the vertex function and using the bare vertex. It is known that in 
the T-matrix approximation in the form discussed here there is no BKT solution below the mean-field 
critical temperature T C MF , i.e. there is no solution for the boson propagator behaving like q~ 2a with 



The region T > Tbkt is evidently more complicated due to the presence of vortices which are singular 
objects. 



a < 2 (see Sec. 6.3).[^] As was shown in [202] for the 4D gauged Nambu-Jona-Lasinio model solving the 
equation for vertex leads to the desired behaviour of the boson propagator, but the solution obtained is 
valid only on the critical line and cannot be analytically continued into a broken phase (which exists in 
this model). This clearly shows that a more elaborate scheme for the T-matrix approximation should 
be investigated. Note also, in this connection, that the importance of the consistency between the one- 
and two-particle properties was stressed in [135]. 

We mention also the recent paper [101] where it is stated that T-matrix approximation does not yield 
the pseudogap behaviour found in quantum Monte-Carlo simulations [108-113] and always leads to a 
Fermi-liquid normal state. 

At the same time there is some evidence for pseudogap behaviour even for the repulsive 2D Hubbard 
model, but within a more refined form of the T-matrix approximation [135]. If the coupling constant 
is not too large, one can convince oneself that the spectral weight reveals a two-peak structure near 
the Fermi level below a crossover temperature, Tx < T 6 MF . According to this paper (see also the more 
recent Ref. [113]) the pseudogap starts to appear (in d = 2, but not in d — 3) when the two-particle 
correlation length becomes larger than the single-particle de Broglie wave length and the system enters 
into a so-called renormalized classical regime, as confirmed both analytically and numerically. A similar 
picture applies to pairing fluctuations slightly away from half-filling in the attractive Hubbard model 
where the pseudogap opens, initially in the region of fc-space where the Fermi velocity is smallest. 

The Gaussian fluctuations destroy long-range order in 2D [128], and if solving equation (4.12) one 
searches for the T c at which long-range order sets in, one obtains zero in accordance with the CMWH 
theorem [126,127]. In some sense this is an excellent check whether one has a good solution to the 
T-matrix equations: the correct solution should give T c = in 2D. 

The most unambiguous acknowledgement of the problem with T c is evident in [106] , where the authors 
phenomenologically introduce three-dimensionality and define the superconducting critical temperature 
to correspond to a factor of 100 enhancement of the superconducting susceptibility, i.e. the value of T c 
follows from the equation (compare with Eq. (4.12)) 

1 _ Ux{q = 0, u = 0) = 0.01 . (4.24) 

One can see that this choice of the definition of T c is entirely arbitrary and with the same success one 
could propose an enhancement of x by a factor of 1000. In fact to avoid problems with the correct 
definition of T c most T-matrix 2D theories are restricted to the discussion of the phase above T c and 
one cannot then investigate the evolution of the pseudogap into the superconducting gap. 

There is no doubt that the quasi-2D character of cuprates must be taken into account and many 
attempts have been made to consider quasi-2D models (see, for example, [27,56,65,203-206]). However 
for a pure 2D theory one cannot use the T-matrix approach to find the superconducting temperature 
and we outline an alternative approach in the next chapter. 



Note that studying the Schwinger-Dyson equations one can speak about the BKT phase transition only in 
terms of the propagators. 



5 The superconducting transition as a BKT transition and the temperature scale for the 
opening of the pseudogap 



5. 1 Phase diagram based on the classical phase fluctuations in the absence of Coulomb repulsion 



As discussed in the previous chapter most previous analyses of the behaviour of 2D systems at T ^ 
have been based on different extensions of the Nozieres and Schmitt-Rink [79] approach (see, e.g. 
[199,93,200,201]). This approach, and its T-matrix extensions as discussed in Sec. 4.5, do not go beyond 
the Gaussian approximation for the pairing fluctuations which may explain the abovementioned diffi- 
culties faced in these calculations [53]. It is an interesting problem in the theory of superconductivity 
to develop a finite temperature crossover approach in 2D but based on the BKT transition. Such an 
analytical treatment did not however prove easy to develop [53]. 

Nonetheless, there has been some progress using approaches which are specifically designed to deal 
with the phase degree of freedom. The initial development of the crossover from superconductivity 
to superfluidity in 2D as a function of the carrier density rif was carried out as long ago as 1992 by 
Drechsler and Zwerger [208] (see also Ref. [172], where the ci-wave case was studied). By means of a 
Hubbard-Stratonovich transformation a statistical Ginzburg-Landau theory was derived. By applying 
the correct normalization to make the order parameter a Bose field the following free energy functional 
results 
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where h = V x A is the magnetic field, c is the speed of light and the effective chemical potential //*, 
effective mass m* and effective coupling strength g* are all determined by the parameters of the original 
Ginzburg-Landau theory. The critical temperature must still however be determined by the introduction 
of a third dimension to avoid divergences in the number equation [208,157,172]. It can then be shown 
that one obtains the correct limiting behaviour for the effective parameters e.g. the chemical potential 
at the critical temperature ji c = fJ>*(T c ) attains the value — |e&|/2 in the Bose limit and ep in the Fermi 
limit. Remarkably one also recovers certain features of BKT physics e.g. in the Bose limit the critical 
temperature is close to the BKT critical temperature for a dilute hard core Bose gas of boson mass 2m 
and density nf/2 on a lattice 

T c ~ Tbkt = 0.89^ (5.2) 
4m 

We defer the discussion of the Nelson-Kosterlitz jump in the full (renormalised) superfluid density 
to Sec. 5.4. As noticed in [117], the method of [208] is fully justified only close to T C MF where the 
linearized Ginzburg-Landau theory is valid. Nor does this method satisfy the important requirement of 
self-consistency [116]. 

A crossover theory based directly on the BKT transition which avoids these problems has however been 
studied in relativistic 2 + 1-theory [207]. We now present this approach which has recently been applied 
to the theory of superconductivity in a continuum model [115] (see also the generalization to the phonon 



h 2 



model (3.61) in [209]). Analogous results, which will also be discussed, have previously been obtained 
for the discrete Hubbard model [117,120,121]. 

For a 2D system one should rewrite the complex ordering field $(x)[^] in terms of its modulus p(x) 
and phase 6{x) i.e. $(x) = p(x) exp[i9(x)], as was discussed by Witten in the context of 2D quantum 
field theory [210] (see also the review [211]), and used in the theory of the superconducting properties 
of quasi-lD systems [212] and thin films [213]. It is impossible to obtain $ = (3>(x)) ^ at finite 
T since this would correspond to the formation of symmetry breaking homogeneous long-range order 
which is forbidden by the CMWH theorem [126]. However, it is possible to obtain p = (p(x)) ^ 
but at the same time $ = p(exp[i9(x)]) = due to random fluctuations of the phase 9(x) (i.e. due 
to transverse fluctuations of the order field proceeding from the modulus conservation principle [214]). 
We stress that p ^ does not imply any long-range superconducting order (which is destroyed by the 
phase fluctuations) and, therefore, is not in contradiction with the abovementioned theorem. At the 
same time one must take into account that p is a collective variable which corresponds to some sort of 
collective behaviour in the 2D system. 

The basic concepts that originate from the modulus-phase representation can be summarized as follows. 
Firstly, as in the paper of [208], one can develop a low-energy theory for the phase fluctuations in terms 
of an effective XY model with a superconducting temperature that is given by the BKT temperature 
[117,120,121,115] (see also [215]). Secondly the amplitude (or modulus) of pairing fluctuations can be 
non-zero [64,207,115,120,121,215] even above the BKT transition temperature and first becomes non- 
zero at roughly the pairing temperature. We now illustrate these concepts using the simple continuum 
model of [115] for which one can do much of the calculation analytically and then discuss the effect of 
the lattice [117] and of quantum phase fluctuations and Coulomb repulsion [120,121,215]. 

We stress that the two concepts above are not specific to the pairing Hamiltonian and/or the BCS-Bose 
crossover but are related to the dimensionality of the system. Any 2D theory which is characterized 
by a single complex order parameter has a corresponding effective low-energy theory in terms of the 
phase of this parameter. The resulting Hamiltonian is often of XY type with a corresponding BKT 
superconducting transition temperature. Furthermore the modulus or amplitude of the order parameter 
can be non-zero above the superconducting transition temperature since this non-zero value does not 
correspond to long-range order. Thus these ideas were first proposed in a phenomenological model by 
Emery and Kivelson [25,26] in which there are two temperature scales as discussed above. 

In the continuum model (and in fact in general) three regions can be identified in the 2D phase diagram. 
The first one is the superconducting (here BKT) phase with p ^ at T < T BKT , where T B kt is the 
temperature of the BKT transition, which plays the role of T c in pure 2D superconducting systems. In 
this region there is algebraic order, or a power law decay of the (<&*<&) correlations. The second region 
corresponds to a phase with non-zero p (Tbkt < T < T p ), where T p is the temperature at which the 
homogeneous (by definition) p becomes zero. In this phase although p is still non-zero the correlations 
mentioned decay exponentially. The third is a normal (Fermi-liquid) phase at T > T p where p = O-PH 
Note that $ = everywhere as is the case for all other correlators which violate the symmetry e.g. 
($(x)$(0)) = 0. Note that while this phase diagram was derived for the idealized 2D model, there 
are indications that even in layered systems as complicated as HTSC the critical temperature may 
nonetheless be well estimated using Tbkt [64,204], even though the transition undoubtedly belongs to 



Because an order parameter is absent in 2D systems we will use this notation in what follows. 
It might be the case that p(x) never becomes strictly zero. In this case one would have crossover behaviour 
instead of a phase transition, see, for example, P. Curty, H. Beck, Phys. Rev. Lett. 85 (2000) 796. 



the 3D XY class (see also [216]). The identification of the intermediate phase with the pseudogap phase 
is somewhat controversial and we will return to this point in some detail. We note however that it 
has been pointed out by Abrikosov [64] that a nonzero gap in the one-particle excitation spectrum can 
persist even without long-range order, being a local characteristic (see the discussion in Sec. 6.1). 

The proposed description of the phase fluctuations and the BKT transition is of course very similar 
to that proposed earlier by Emery and Kivelson [25,26]. However in their phenomenological approach 
the field p(x) does not appear explicitly while in the present microscopic approach it occurs naturally. 
We mention here also the application of similar ideas to the 3D case [216], where instead of the 2D 
temperature T B kt one has the temperature of the phase transition in the 3D XF-model, T c xy . 

There is no need to write down the model Hamiltonian which is studied here, since it is identical with 
that described in Sec. 3.1. 

The desired phase diagram consisting of normal, anomalous and superconducting phases was calculated 
in [115] employing the Hubbard- Stratonovich method (see Sec. 3.1.1, the equations (3.6) - (3.8)). In 
the 2D case at nonzero T, however, instead of using the accepted method for the calculation of the 
partition function Z(v,p,T, ) (see (3.6)), one must perform the calculation in modulus-phase variables. 
This allows one to avoid any subsequent treatment of the phase fluctuations at Gaussian level only. 
Thus, one is able to take into account the phase degree of freedom with the needed accuracy. 

The modulus-phase variables were introduced in accordance with [217], where the parameterization 

$(x) = p{x)e ie{x \ $*(x) = p{x)e- iei - x \ (5.3) 

was used. At the same time as this replacement (5.3) is implemented, one makes the gauge transforma- 
tion 

iM*) = xA*Y e(x)l \ itf (*) = xl{x)e~ W{x)l2 . (5.4) 

[210] (see also Refs. [211,213,217]). Physically, this amounts to replacing the charged fermion ip a (x) with 
a neutral fermion x<r( x ) an d a spinless charged boson e ie ^ x ^ 2 . These parameterizations, (5.3) and (5.4), 
prove more appropriate for presenting the corresponding functional integral in two dimensions [210]. In 
Nambu variables (3.2) the transformation (5.4) takes the following form 

y( x ) = e iTs6 ^ 2 T(x), V\x) = Y^e"^^ 2 . (5.5) 

Note that the neutral fermion variables are still Grassmann fields since they satisfy the appropriate 
anti-commutation algebra. 

Making the corresponding substitutions (5.5) into the representation (3.6) and integrating over the 
fermi-fields T and we arrive at the expression (compare with (3.7) and (3.8)) 

Z(v,n,T) = J pVpV9exp[-pn(v,n,T,p(x),de(x))], (5.6) 

where 

PSI(v, p, T, p(x), Mx)) = ^ ! dr [ d 2 rp 2 (x) - Tr LnG" 1 + Tr LnG 1 (5.7) 

U Jo J 

is, as (3.8), the one-loop effective action, which, however, depends on the modulus-phase variables. The 
action (5.7) is expressed through the Green's function of the initial (charged) fermions that has in the 



new variables the following operator formp 5 ] 
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(5.9) 



The free fermion Green's function Go = G\^^ p fi = o provides a convenient regularization in the process of 
calculation. It is important that neither the smallness nor slowness of the variation of the phase of the 
order parameter is assumed in obtaining expression (5.7). In other words, it is formally exact. 

Since the low-energy dynamics of phases for which p ^ is governed mainly by long- wavelength fluc- 
tuations of 9(x), only the lowest-order derivatives of the phase need be retained in the expansion of 
n(v,fi,T,p(x),dO(x)): 



where 



Q(v, p, p(x), dO{x)) ~ n kin (w, p, T, p, dO(x)) + n po t(v, p, T, p), 



(5.10) 



tt kin (v,p,T,p,d9(x)) = TTr£ -(££) 



n=l 



n 



p=const 



(5.11; 



and 



Q pot (v,p,T,p)= Qj J d 2 rp 2 -TTrlng- 1 + TTrlnG 1 



p=const 



(5.12) 



The kinetic i\m and potential Q pot parts can be expressed in terms of the Green's function of the 
neutral fermions, which satisfies the equation 



~ IdT + Ts 2m" + ^ +TlP 



G(r, t) = 6(t)6(t) 



(5.13) 



and the operator (5.9). 



The representation (5.10) enables one to obtain the full set of equations necessary to find Tbkt, p(7bkt), 
and /x(Tbkt) at given ep (or, for example, p{T) and p{T) at given T and ep). While the equation for 
Tbkt will be written using the kinetic part (5.11) of the effective action, the equations for p(Tbkt) 
and /x(Tbkt) ( or p(T) and p(Tj) can be obtained using the mean field potential (5.12). It turns out 
that in a phase for which p ^ 0, the mean-field approximation for the modulus variable describes the 
system quite well. This is mainly related to the non-perturbative character of the Hubbard-Stratonovich 
method, i.e., most of the pairing effects are included in the nonzero value of p. 

It is clear that the CMWH theorem [126] does not preclude nonzero (p) and, as a consequence, an 
energy gap for fermion x, since no continuous symmetry is broken when such a gap appears. Despite 



25 It may be obtained as a solution of a differential equation with anti-periodic boundary conditions (see (3.£ 
and (3.10)). 



strong phase fluctuations in the two-dimensional case, the energy gap in the spectrum of the neutral 
fermion x still persist in the spectrum of the charged fermion if) (see Chap. 7 and [125]), even well above 
the critical temperature.p 6 " 



Thus, within the proposed scenario, the pseudogap properties might conceivably be attributable to the 
energy gap of a neutral fermion introduced in the way described above, so that the pseudogap itself 
can be considered as a remnant of the superconducting gap. The condensate of neutral fermions is 
completely unrelated to the superconducting transition; the latter is only possible when the superfluid 
density of bosons becomes large enough to stiffen the phase 6(x). The temperature T p at which nonzero 
(p(x)) develops should be identified in this approach with the pseudogap onset temperature T* (see 
Sec. 5.1.2 item a and the discussion about heavily underdoped region in Sec. 5.1.3). 

The strategy of treating charge and spin degrees of freedom as independent seems to be a useful, and 
at the same time general feature of two-dimensional systems [210,213] and low-dimensional systems in 
general (see also the review [218]). 

5.1.1 Derivation of self- consistent equations for Tbkt, neutral order parameter, and chemical potential 

When the model under consideration is reduced to some known model describing the BKT phase 
transition, one can easily write down an equation for Tbkt, which in the present approach can be 
identified with the superconducting transition temperature T c . Indeed, in the lowest orders the kinetic 
term (5.11) coincides with the classical spin XV-model [130-134], which has the continuum Hamiltonian 



Here J is the some coefficient (in the original classical discrete XV-model it is the stiffness of the 
relatively small spin rotations) and 9 is the angle (phase) of the two-component vector in the plane. 

The temperature of the BKT transition is, in fact, known for this model: 



Despite the very simple form of Eq. (5.15), it was derived (see, e.g., Refs. [130-134]) using the renor- 
malization group technique, which takes into account the non-single-valuedness of the phase 9. Thus, 
fluctuations of the phase are taken into account in a higher approximation than Gaussian and this 
approximation may well be better than the T-matrix even in its fully self-consistent conserving form. 
However, the disadvantage of such an effective theory of the phase fluctuations is that one loses contact 
with the underlying fermions (they are essentially integrated over to obtain the effective theory^). 
From this point of view, these two approaches are essentially complementary. 

The XY~-model was assumed to be adequate for a qualitative phenomenological description of the 
underdoped cuprates [25-27] (see also Refs. [66,120,221]), and the relevance of the BKT transition 

26 We note that the specific heat experiments [4,107] demonstrated a loss of entropy at temperatures much 
higher than T c . This can be considered indicative of a degenerate normal state, consistent with the existence 
of a nonzero order parameter {p{x)). 

27 This contact, however, would not be lost if one retained the source terms for the fermi-fields in the corre- 
sponding functional integral. See also the approach of [124] where the low energy fermions are not integrated 
out. 




(5.14) 



Tbkt — — J- 



(5.15) 



to Bose- and BCS-like superconductors was recently discussed in Ref. [172] (see also [222]). The most 
important difference between superfluidity in 4 He films, which is adequately described by the XF-model 
[130], and superconductivity is that the superconducting liquid is charged. This results in coupling to the 
magnetic field, which inevitably accompanies any current flow in the superconductor, e.g. in the presence 
of vortices. However, as was shown in [222,172] the effective penetration depth for the magnetic field in 
thin films is of the order of 1 cm, and thus for sample sizes smaller than this, magnetic screening becomes 
irrelevant. In such a situation the difference between charged and neutral superfluids becomes irrelevant 
and superconducting films may undergo the BKT transition. Note also that this difference leads also to 
electric coupling between Cooper pairs, in particular those situated in different superconducting layers. 
The effect of this coupling is considered below in Sec. 5.3. 

To expand f2k in up to ~ (V6>) 2 , it is sufficient to restrict ourselves to terms with n — 1,2 in the 
expansion (5.11). The calculation is similar to that employed in [217], where only high densities rif were 
considered at T = 0. (We mention also the derivation of the effective action at T = in [223] where the 
fluctuations in the density (conjugated to the phase 9) were also treated.) Thus, to obtain the kinetic 
part, one should directly calculate the first two terms of the series (5.11), which can be formally written 
= TTr(£E) and = ±TTr(£E£E). We note that E has the structure E = r 3 Oi + I0 2 , where 
0\ and 2 are differential operators (see (5.9)). One can see, however, that the part of E proportional 
to the unit matrix / does not contribute to il k ; n . Hence, 



= T l "rj^j^Zj^^M + i^-) , (5.16) 
where (compare with (3.15)) 

r/;,, Jr\ W + r 3 £(k) - np 

^> k ) = - ^+ak)+P 2 (5 - 17) 

is the Green's function of neutral fermions in the frequency-momentum representation, with £(k) and 
^(k) defined after (3.15). The summation over the fermionic Matsubara frequencies uo n = 7i(2n + 1)T 
and integration over k in (5.16) can be easily performed using the sum (A. 7); one thus obtains 

id T e (yef 



^ = r/>/,W,,r,p) ^ + ^L , (5 , 8) 



where 
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n F (p, T, p(/i, T)) = — vV + p 2 + fi + IT In 



1 + exp 



T 



(5.19) 



This has the form of a Fermi quasi-particle density (for p = the expression (5.19) is simply the density 
of free fermions given by the first term uq in (4.18)). 

For the case T = [217], in which the real time t replaces the imaginary time r, one can argue from 
Galilean invariance that the coefficient of d t 9 is rigorously related to the coefficient of (V6>) 2 . It therefore 
does not appear in fij^. We wish, however, to stress that these arguments cannot be used to eliminate 
the term (V6 1 ) 2 from when T ^ 0, so that we must calculate it explicitly. 



The 0\ term in E yields 



rp ,p , rp oo „ 

U$L{0{) = --] Q drjd\— J d 2 kTr[g(iu n ,k)r 3 g(iu n ,k)r 3 ) 
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Using (A. 11) to calculate the sum over the Matsubara frequencies, we find that 



(5.20) 



n^(Oi) = ~jyrJd 2 rK( t ,,T,p)^d T e + ^) , (5.21) 



where 



Kfr T, p{p., 71) = £ (l + ^ tanh . (5 , 2) 



Obviously, the Oi term does not affect the coefficient of (V#) 2 . Further, it is easy to make sure that the 

cross tei 
yieldF] 



(2) (2) 

cross term involving 0\ and O2 in fl^j, is absent. Finally, calculations of the O2 contribution to 



nffiW =T f dr f d 2 r-f- £ / rf 2 A; k 2 Tr[^(^„, k)/^(^ n , k)/]^£. (5.23) 
Thus, summing over the Matsubara frequencies (see Eq. (A. 12)), one obtains 

^gi(0 2 ) = ~ f'dr f d\ — I— [ d*k ^ (W) 2 . (5.24) 



7 128vr 2 m 2 
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e 2 (k) + p 2 



2T 

As expected, this term vanishes when T — > 0, but at finite T it is comparable with (5.18). 
Combining (5.18), (5.24), and (5.21) we finally obtain 

fikm = \\l dr \ rfV [^0, T ; P)W + J(ji, T, p) (W) 2 + K(ji, T, p) (<9 T fl) 2 ] , (5.25) 

where 

J (p, T, p(p, T)) = -^-n F (n, T > P) - 7- / rfx x + /j,/2T 

Am 47T J-n/2T / „2 

cosh \l x2 + ^j7i 

characterizes the phase stiffness and governs the spatial variation of the phase 0(r). One can see that 
the value of the phase stiffness J(T = 0) coincides with the nonrenormalized stiffness used in Ref. [25]. 

During the derivation of the effective action for the phase field we have neglected the so-called Landau 
terms (see [224] and Refs. therein) which have a non-local form and thus cannot be presented in terms 
of the derivative expansion. These terms are important since they make the #-mode decaying for any 
T^O. However, as recently shown in [224], for s-wave superconductors their contribution is negligible 



28 Derivatives higher than (V#) 2 were not found here. 



at least for T < 0.6T P . This is not the case for ci-wave superconductors where these terms are expected 
to be important even in the low temperature region. 



The quantity J(p,T,p) vanishes at p = 0, which means that above T p the modulus-phase variables are 
meaningless; to study the model in this region one must return to the old variables $ and $*. Near T p 
one can obtain from (5.26) in the high-density limit (see below) 

J(p ^e F ,T^ T;,p^ 0) = ^jf^ F ~ 0.016^e F , (5.27) 

where ( is the zeta function. 

Direct comparison of (5.25) with the Hamiltonian of the XF-model (5.14) makes it possible to write 
Eq. (5.15) for T BK t directly: 

Tbkt = | J(n, Tbkt, p(/i, T BKT )). (5.28) 

Although mathematically this reduces to a well-known problem, the analogy is incomplete. Indeed, 
in the standard XF-model (as well as the nonlinear cr-model) the vector (spin) subject to ordering 
is assumed to be a unit vector with no dependence on T. In this case this is definitely not the case, 
and a self-consistent calculation of T B kt as a function of nj requires additional equations for p and p, 
which together with (5.28) form a complete set. We note, however, that the dependence of J on T was 
neglected in [25], where the nonrenormalized phase stiffness J(T = 0) was used to write Eq. (5.28). 
Using this nonrenormalized phase stiffness one would obtain that Tbkt ~ £f which is, as we will see, 
apparently wrong for high carrier densities since one does not recover the BCS limit. 

Using the definition (5.12), one can derive the effective potential Q pot (v : p, T, p) (see Appendix A). Then 
the desired missing equations are the condition 

dQ pot (p)/dp = that the potential (A. 10) be minimized, and the equality v _1 <9f2 po t/<9/i = —rif, which 
fixes rif. These are, respectively (compare with Eqs. (3.21) - (3.24)) 



' = ■ 1 tanhV e(k)+ ^, (5.29) 

U J (2tt) 2 2 / e(k)+p2 2T y ' 



n F (p,T,p) = n f , (5.30) 

where uf(p,T,p) is defined by (5.19). 

Equations (5.29) and (5.30) comprise a self-consistent system for determining the modulus p of the 
order parameter and the chemical potential p in the mean-field approximation for fixed T and rif. 

It is very interesting to note that to resolve the problems which are present in the non-conserving, non 
self-consistent approach (see Sec. 4.2), the following form of the phenomenological pair susceptibility 
(compare with Eqs. (4.2) and (4.15)) was suggested in [116] 



(5.31) 



where (n^) is a function with a free parameter tp, and has the form (compare with (5.19)) 



<n k > = g 



: tanh 



IT 



(5.32) 



It was stressed in [116] that tp is not an off diagonal long-range order parameter but is related to 
pair binding. The phenomenological susceptibility (5.31) was then used in [116] to write down the 
thermodynamical potential fi[0, G] (see 4.6) where $[G] (see 4.7) is given in terms of (5.31). This 
thermodynamical potential was then used to derive equations for the chemical potential (which coincides 
with (5.30)), entropy and specific heat. An additional equation for tp was derived in [116] by minimizing 
this potential with respect to tp subject to the number constraint 2][] k (n k ) = nj. It was shown in [116] 
that the presence of tp avoids problems with negative entropy and specific heat which are present in the 
non- conserving, non self-consistent approximation discussed in Sec. 4.2. 

Eqs. (5.29) and (5.30) seem to yield a reasonable approximation at high densities rif, since they include 
the condensed boson pairs in a non-perturbative way via nonzero p. Nonetheless they must certainly 
be corrected in the strong coupling regime (low densities rif) to take into account the contribution 
of noncondensed bosons (this also appears to be important for Eq. (5.28), which determines Tbkt)- 
The extent to which this alters the present results is not completely clear. Previously, the best way to 
incorporate noncondensed pairs appears to have been the self-consistent T-matrix approximation (see 
e.g. [55,91,93,95] and other references in Sees. 2.3 and 4), which allows one to account for the feedback 
of pairs on the self-energy of fermions. However as we discussed above, the T-matrix approach, at least 
in its standard form presented in Sec. 4, fails to describe the BKT phase transition, for which one must 
consider the equation for the vertex. On the other hand, in this approach the BKT phase transition is 
realized by the condition (5.15), while an analog of the T-matrix approximation in terms of propagators 
of the p-particle and the neutral fermion x has yet to be developed. 

Again the energy of two-particle bound states in a vacuum defined by (3.26) is more convenient to use 
than the four-fermion constant U . For example, one can easily take the limits W — > oo and U — > in 
Eq. (5.29), which after this renormalization becomes 



In 



du 



lu 2 + 
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(5.33) 



Thus, in practice, we solve Eqs. (5.28), (5.30), and (5.33) analytically and numerically to study Tbkt 
as function of n/ (or equivalently, of the Fermi energy ep = imfjm). It is easy to see that at T = 0, the 
system (5.30), (5.33) transforms into a previously studied system (3.29). Setting p = in Eqs. (5.29) 
and (5.30), we obtain (in the same approximation) the equations for the critical temperature T p : 
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v/ 2T p tanh it 
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(5.34) 



where 7 ~ 1.781 (C = In 7 ~ 0.577 is the Euler's - Mascheroni's constant) and the corresponding value 
of /i: 



T p ln 




(5.35) 



Note that these equations coincide with the system that determines the mean-field temperature T^ 2D ^ MF (= 
T p ) and p(T^ 2D ^ MF ), evidently as a result of the mean-field approximation for the variable p used here. 
Thus in this approximation the average amplitude p of the pairing fluctuations is simply given by 
the mean-field gap A MF which is non-zero below y c ( 2£) ) MF (= T p ). This non-zero amplitude then plays 
the role of the "pseudogap". Note that the paramagnetic susceptibility (see [115]), specific heat in the 
mean-field approximation exhibit kink-like behaviour at T p while at the superconducting transition 
temperature Tbkt they do not show any singularities. The picture presented here is thus to all practical 
purposes identical to that given in [120,121]. 

There is an important difference between the temperatures T 2D and T p . Specifically, if one takes fluc- 
tuations into account, T 2D goes to zero, while the value of T p remains finite. The crucial point is that 
the perturbation theory in the variables p and 9 does not contain any infrared singularities [210,219], 
in contrast to the perturbation theory in $, $*; thus the fluctuations do not reduce T p to zero. This is 
why the temperature T p can have a physical meaning: incoherent (local or Cooper) pairs begin to form 
(at least at high enough rif, see Sees. 5.1.2 item c and 5.1.3) just below T p . At higher temperatures, 
only pair fluctuations exist (see, e.g. [220]). 



5.1.2 Main results 



The interested reader may find a more detailed discussion of the results for the continuum model in 
[115]. Here we present only the issues directly related to the phase diagram of the model (3.1). A 
numerical and analytical investigation of the systems (5.28), (5.30), (5.33), and (5.34), (5.35) yields the 
following results, which are displayed graphically as the phase diagram of the system. 

a) For rather low carrier densities, the pseudogap phase area (see Fig. 15) is comparable with the BKT 
area. 



T c /\e b \ 




Fig. 15. Tbkt and T p versus the noninteracting fermion density (taken from [115]). Dots represent the function 
p{sf) at T = Tbkt- The regions of normal phase (NP), pseudogap phase (PP), and BKT phase are indicated. 



For high carrier densities (e^ > 10 2 — 10 3 |efe|), inserting into (5.27) the well-known dependence of p(T) 



(see, for example [49]) 



^ T - T ^ - m T > t 1 - S (5 - 36) 

and then substituting (5.27) into (5.28) one obtains the following asymptotic expression for the BKT 
temperature (also given in [222,107,216] and [209], where the corresponding asymptotical expression for 
the boson-exchange model (3.61) was obtained) 

/ 4T \ 

Tbkt = T p (l - — p -j , Tbkt < T p , (5.37) 

which shows that the pseudogap region shrinks asymptotically. This behaviour qualitatively restores 
the BCS limit. This is a good place to consider the difference between the two high-density limits shown 
in Fig. 2-A and -B respectively. As one can easily see, Eq. (5.37) implies that the present results are 
in favour of the picture presented in Fig. 2-B.[^] The merger of T* and Tbkt in the approach based 
on the superconducting fluctuations is probably not that surprising since both these temperatures have 
the same superconducting origin. However, as one can see from Fig. 2-A this merger would be rather 
unusual if one accepts the antiferromagnetic pseudogap scenario and considers T* to have a magnetic 
nature, different from that of T c . Thus, as was suggested in [23], additional studies on the overdoped 
cuprates in a magnetic field may also help to clarify the pseudogap origin. 

It has been shown that for optimal doping, the dimensionless ratio can be estimated as ep/\sb\ ~ 
3 • 10 2 — 10 3 [146]. This implies that the size of the pseudogap phase predicted by (5.37) to be of order 
Tg KT /e_F [222,107,216] is (even in the underdoped region) much smaller than the difference T* — T c 
observed in cuprates. One might however be able to explain this discrepancy either by the presence of 
impurities as done in Sec. 5.4 or by the formation of spatial inhomogeneities e.g. stripes. 

b) For ep < (10 — 15)|£b|, the function T B kt(cf) is linear, as is confirmed by the analytic solution 
of the system (5.28), (5.30), and (5.33), which yields Tbkt = €p/8. As we discussed in Sec. 2.2 (see 
also the Uemura plot in Fig. 3) such behaviour of T c (ep) is observed for all families of HTSC cuprates 
in their underdoped region though with a smaller coefficient of proportionality (0.01 — 0.1). Again 
this discrepancy can be in principle explained by the presence of non-magnetic impurities (see Sec. 5.4). 
Other reasons such as, for example, a contribution due to noncondensed pairs in Eq. (5.28) which defines 
Tbkt, could also be important. 

c) Finally the calculations showed (see Fig. 16) that the ratio 2A/Tbkt for the carrier densities shown 
in the figure is greater than 4.7 and grows as the doping decreases. This result is also confirmed by 
the more general consideration of Abrikosov [64] which is also based on phase fluctuations of the order 
parameter but takes into account interlayer tunnelling. This concentration behaviour is consistent with 
numerous measurements of the ratio 2A/T C in HTSC [4]. The value 2A/T P (= 2A/T C MF ) is, however, 
somewhat lower and reaches the BCS theory limit of 3.52 only for ep > \e\\. It is very interesting 
and important that the measurements of [32] show that for all doping levels 2A/T* ~ 4 — 5 which 
is very close to the BCS value, 2A^/T C MF ~ 4.3 but for a <i-wave superconductor. This presents a 
strong argument in a favour of the theories which relate the pseudogap formation at T* to the precursor 
superconducting fluctuations. In this case the scale for T* is the BCS temperature T C MF , or as in [115] T p 

29 Of course, Tbkt does not reach zero at high doping due to the oversimplified character of the model. For 
example, we saw in Sec. 3.1 that Tbkt may decrease due to band filling or at the very least stabilizes as in the 
model with phonon interaction from Sec. 3.5. Nonetheless one clearly sees that T p and Tbkt are merging. 
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Fig. 16. 2A/Tbkt and 2A/T p versus the non-interacting fermion density (taken from [115]). 

which, as we showed, has a quite rigorous definition in the modulus-phase approach. Another important 
conclusion which one can draw from the experimental values of 2A/T* [32] is that the simple Bose-BCS 
crossover theories with pre-formed local pairs existing in the region with \xjep < 1 (or even \x < 0) 
with e F /\e b \ < 1 are unable to give 2A/T p < 4 since Fig. 16 clearly shows that the ratio grows for 
very low carrier densities. As we discussed in Sec. 3.1.5 one may still have local Bose pairs in a more 
sophisticated picture where the system is inhomogeneous (see also the discussion in Sec. 5.4 about the 
differences between the approaches of Uemura [84] and Emery-Kivelson [25,26]). 

Note also that the divergence of 2A/T BK t and 2A/T p at e F — > is directly related to the definition of 
A at \i < 0. 

5.1.3 Pairing temperature T pair versus carrier density. More about the meaning ofT p 

There is no disagreement concerning the asymptotic behaviour of Tbkt (or T c shown in Fig. 1) propor- 
tional to ep in the region of very low carrier densities. In contrast, the behaviour of the temperature 
T*, below which we assumed pairs are formed, cannot be considered to be generally accepted. It is an 
open question as to how it behaves in the heavily underdoped cuprates: namely whether it increases 
being related to the Neel temperature Tn or decreases. The opinion of theoreticians on this subject is 
also divided. 

For example, on the phase diagram suggested by Uemura in Fig. 4, the temperature T* is taken to be 
the temperature T pair of local uncorrelated pairing which increases with decreasing rif. Randeria (see 
Ref. [53] and references therein), to define the pairing temperature T pair , uses the system of equations 
for the mean-field transition temperature and the corresponding chemical potential, which is essentially 
identical to the system (5.34), (5.35). Thus his T pair — > as nf — > 0. This is accordance to the behaviour 
of T* shown in the theoretical phase diagram from Fig. 1 suggested by Zachar [23]. 

It is also well known [53,90,199] that in the low-density limit, it is vital to include quantum fluctuations, 
at least in the number equation [79], in the calculation of the critical temperature at which long-range 
order forms in 3D. 

Certainly quantum fluctuations are also important in the calculation of T p in the limit nf — > and, 
in particular, in the number equation. However, as already stressed in Sec. 5.1.1, these corrections are 




quite different from that obtained using the variables $, $*, since perturbation theory in the variables 
p and 9 does not contain any infrared singularities [210,219], and the fluctuations do not yield T p = 0. 



In our opinion, the temperature T p has its own physical interpretation: this is the temperature of a 
smooth transition to the state in which the neutral order parameter p ^ 0, and below which one 
can observe pseudogap manifestations. In this respect the temperature T p can be compared with the 
crossover temperature Tx introduced for the repulsive and attractive Hubbard models in [135] (see also 
the explanation in Sec. 4.5). At the same time note that both these temperatures should not be strictly 
identical to the mean-field temperature T C MF (see the end of Sec. 5.1.1). 

There is also a very interesting and important question about the character of the transition (see, for 
example, [215]). Certainly in the simplest Landau theory one appears to have a second-order phase 
transition, since p takes a nonzero value only below T p [207] and since p(x) has only been treated in the 
mean-field approximation i.e. one has neglected the fluctuations in both p(x) and 9(x), a second-order 
phase transition was obtained at T p . However, as it was stressed in Introduction, experimentally the 
formation of the pseudogap phase does not display any sharp transition and the temperature T* observed 
in various experiments is to be considered as a characteristic energy scale and not as a temperature 
where the pseudogap is reduced to zero [4,31]. We believe that taking into account the fluctuations of 
p(x) may resolve the discrepancy between the experimental behaviour of T* and the temperature T p 
introduced in the theory. 

To define the temperature T pa i r properly, one should study the spectrum of bound states either by 
solving the Bethe-Salpeter equation [141] (see also Sec. 3.1.4) or by analyzing the corresponding Green's 
functions as is done here. It turns out that there is no difference between T pair and T p in the Cooper pair 
regime (p > 0), while in the local pair region (p < 0) these temperatures exhibit different behaviour. 

Indeed, let us study the spectrum of bound states in both the normal (p = 0) and pseudogap (p ^ 0) 
phases. We are especially interested in determining the conditions under which real bound states (with 
zero total momentum K = 0) become unstable. For this purpose one can look at the propagator of the 
p-particle in the pseudogap phase: 



x l(35 2 n(v,p,T,p(r,r),d9(r,r)) 
T p (r,r) = - — 



Sp(r,r)5p(0,0) 



(5.38) 

P=Pmin=COnst 



where p min is defined by the minimum condition (5.29) (or (5.33)) of the potential part (A. 10) of 
the effective action (5.7). In the momentum representation, the spectrum of bound states is usually 
determined by the condition 

0, (5.39) 



where T^u, K) is the retarded Green's function obtained directly from the finite temperature Green's 
function (5.38) in frequency-momentum representation, T p (iQ n , K) using analytical continuation iQ n — > 
uj + iO. Recall that such analytical continuation must be performed after evaluating the sum over the 
Matsubara frequencies. In the case of vanishing total momentum K = 0, one arrives at the energy 
spectrum equation 



1 _ 1 r d 2 k g(k) tanhy^W+p 2 ^ 
T p ^,0)-f/ + J (27T) 2 / e(k)+p2 ^-4[e(k)+p2] 



= 0. (5.40) 



From the explicit expression (5.40) for T^{u, 0), this function obviously has a branch cut at frequencies 



M > 2min v /£ 2 (k) + p 2 = 

Thus, bound states can exist only below this cut. 

Real bound states decay into two-fermion states when the energy of the former reaches the branch 
point 2 min^£ 2 (k) + p 2 . Since 1/Tf (cj,0) is a monotonically decreasing function of u 2 , the solution for 
the bound-state energy must be unique and is given by uo = — £&(T) = 2p(T). At this point Eq. (5.40) 
coincides exactly with the mean-field equation (5.29) for p(T) so that p{T) takes on the mean-field 
value. It is also clear that for p < we have real bound states with energy £&(T) = — 2p below the 
two-particle scattering continuum at uj — 2v ' p 2 + p 2 , while at p > there are no stable bound states. 
The line p(T,e F ) = in the T-e F plane at p ^ separates the negative p region, where local pairs 
exist, from that in which only Cooper pairs exist (positive p). This line (see Fig. 17) begins at the point 
T = (e 7 /7r)|efe| ~ 0.6|£fc|, e F ~ 0.39|eb| and ends at T = 0, e F = |£&|/2. (The latter follows directly from 
the solution at T = given by (3.27).) 



2p, p>0 
2^p 2 + p 2 , p<0. 



(5.41) 
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Fig. 17. Phase diagram of the 2D metal at low concentrations (taken from [115]). The dotted line corresponds 
to p = 0, and the pairing temperature T pa i r (squares) separates the abnormal normal phase (ANP) from the 
normal phase (NP). The region below the solid line is the pseudogap phase (PP). The critical temperature 
Trkt is not shown. 



To find the equivalent line in the normal phase with p = i.e. above T p , we consider the corresponding 
equation for the bound states. The propagator of these states (in imaginary time formalism) is defined 
to be 



T 



_ x f36 2 Q(v,p,T,^(T,r),^(r,r)) 



5$*(r, r)5$(0,0) 



(5.42) 



(In the normal phase, where p = 0, we must again use the initial auxiliary fields $ and $*.) Then in 
the momentum representation (after summing over the Matsubara frequencies) we have 



T-\iQ n ,K) 



1 1 r d 2 k tanh£ + (k,K)/2T + tanh£_(k,K)/2T 
U~2J (2^)2 f + (k,K)+f_(k,K)-if2 n 



(5.43) 



where k is the relative momentum of the pair. It is seen that Green's function (5.43) coincides with the 
T-matrix T (iQ n ,~K) (compare with Eqs. (4.14) and (4.15)) in the non self-consistent, non-conserving 
approximation, up to an overall minus sign. The spectrum of bound states is given again by Eq. (5.39). 
Using the energy e h (see Eq. (3.26)) of the bound state at T = 0, for K = we obtain the following 
equation for the energies of these states in the normal phase: 



Such states can exist provided —2/j, — \sb\ < uj < — 2/i. The left-hand side of Eq. (5.44) is positive at 
uj = —2fi — \ej,\ and tends to +oo (p > 0) or — oo (p < 0) when uj — > —2/i. This equation thus always 
has a solution between —2/i — and —2/i for /i < 0, i.e. bound states with zero total momentum exist 
for negative p. 

For p > 0, the analytic analysis becomes more complicated, and requires numerical study. One can 
easily find from (5.44) that for the special case T = 0, which is not in fact in the normal phase but for 
which one can use the variables $ and $*, that such stable bound states exist up to p < |£&|/8 (see 
also proof in Sec. 3.1.4). In fact, numerical study for T > T p shows that the trajectory p(T, e_p) = 
(or T = ei?/ln2, see (5.35)) approximately divides the normal phase into two qualitatively different 
regions - with (/i < 0) and without (/i > 0) stable (long-lived) pairs. This result is the same as that 
in the phases with p ^ 0, so that the line /i(T, ep) = (Fig. 17) separates the regions with stable and 
unstable pairs for all non-zero temperatures. However this line (squares) does not correspond to any 
phase transition. 

Given the two-particle binding energy as a function of temperature, it is natural to define the pairing 
temperature T pair by T pair ps |£ft(T pair , /i(T pair , e^))|. This equation can be easily analyzed in the region 
ep <C \eb\, for which we directly obtain T pair ps \eb\, which clearly coincides with the standard estimate 
[84,95]. This means in turn that the curve T pair (eir) starting at T pair (0) ~ will decrease until the point 
2p a i r (0.39ei?) ~ 0.6|e&|, which lies on the line T p {ep) (see Fig. 17). It is important that this line is not 
a phase transition curve; it merely divides the fermion system diagram into temperature regions where 
one has either a prevailing mean number of local pairs (T < T pair ) or unbound carriers (T > T pair ). This 
is the region of the abnormal normal phase where one has some density of pre-formed boson pairs. It is 
widely accepted, however, that this case is only of theoretical interest, since there is no Fermi surface 
(/i < 0) in this phase (see nevertheless the discussion in Sees. 1.3, 3.5.3 and Refs. [51,52,42,43]). The 
phase area or the difference Tp(e^) — T p (ep) is an increasing function as ep — > 0, which corresponds to 
the behaviour usually assumed [84,25]. 

When /i > there are no stable bound states (£&(T) = 2p(T) = 0) in the normal phase. Formally, 
using p(T) = in Eq. (5.40), we immediately obtain (5.34) or, in other words, here T pa j r = T p . Such a 
conclusion is in accordance with the generally accepted definition of T pair in the BCS case [95]. 

Thus the phase diagram of a 2D metal above T c acquires the form shown in Fig. 17. It is interesting 
that if there are no stable pairs i.e. T pa ,i r (ep) is not well defined, then the temperature T p (€f) is the line 
below which pairs reveal some signs of collective behaviour. Moreover, for T < T p one can speak of a 




(5.44) 



real pseudogap in the one-particle spectrum, while in the region T p < T < T pair only strongly developed 
pair fluctuations (a finite number of pairs) exist, although they probably suffice to reduce the spectral 
quasi-particle weight, and to produce other observed manifestations that mask pseudogap (spin gap) 
formation. 



5.2 The peculiarities of the phase diagram for the lattice model 

The discrete version of the attractive Hubbard Hamiltonian (3.47) was also mapped [117] into the 
classical XF-model. Although we disagree on some of the fine points in the derivation we agree entirely 
with the results. It is instructive to consider these fine points since they yield a better understanding 
of the modulus-phase formalism described in Sec. 5.1. The statistical sum (3.7) is written in [117] for 
(3.47) in the momentum representation. Then the saddle-point approximation described by (3.12) is 
considered. Of course, this approximation leads to nothing but the finite temperature generalization 
of the zero temperature equations for A (3.21) and \i (3.23) considered in Sec. 3.1.3. The value A in 
these equations is called in [117] "order parameter" in accordance with the saddle-point character of 
the derived approximation for the complex Hubbard-Stratonovich fields $ and $*. These equations are, 
of course, formally identical to Eqs. (5.29) and (5.30) derived in Sec. 5.1.1, but as we have discussed 
the meaning of p is completely different allowing nonzero p but keeping ($) = p(exp(id(x))) = 0. 

To go beyond the saddle-point approximation, it is proposed in [117] to consider a variation of the 
Hubbard-Stratonovich fields around their saddle-point values. In doing this, one looks for the low 
energy excitations. It is claimed in [117] that choosing the order parameter to be real (or taking it 
any specific value), one breaks its continuous symmetry, provided A / 0. According to Goldstone's 
theorem, this symmetry breaking implies that a soft mode exists as q — > 0. These modes are associated 
in [117] with the phase fluctuations of the order parameter. It is clear, however, that this breaking of 
the continuous symmetry and appearance of the Goldstone mode in 2D at T ^ are strictly forbidden 
by the CMWH theorem [126]. In fact the above conclusions of [117] are obviously related to the use of 
the saddle-point approximation in the form described in Sec. 3.1.1. To be consistent with the CMWH 
theorem one should (following by Witten [210]) use the set of transformations (5.3) and (5.4) for any 
2D theory. As we have explained in Sec. 5.1 these transformations make the perturbation theory in the 
modulus-phase variables free from symmetry violating terms. Furthermore, as we will see in Sec. 6.3 in 
2D the propagator associated with the phase fluctuations does not have the canonical behaviour ~ 1/q 2 
and becomes softer. 

All these comments do not however change the results of the calculations in [117] but rather reflect 
subtle points in their interpretation. Thus we present the effective action for the phase mode derived in 
[117H 

^eff[0] =E Jeff(q)0q0-q (5-45) 

q 

with 



Note that the phase was assumed to be time-independent in [117]. 



j eS (ci) =—h-uY: 



(u + u_ + v + v_) E+ + E — 

f , - f ~n 

(5.46) 

where A is the solution of the mean-field gap equation (see e.g. Eq. (5.29)), the subscript ± indicates 
k ± q/2, u\ = |(1 + £±/E±), v± = |(1 — £±/E±), and f± = f(E±) is the usual Fermi functions with 




Ek = + A 2 - A s i n Sec. 5.1.1 one may check that J e fj is nonzero only below the mean-field (BCS) 
critical temperature where A^O. 

To provide a link with the XF-model, one expands J e s up to second order in q to obtain the low energy 
behaviour 

H = J£q 2 ^- q — J I dr[Ve{r)} 2 . (5.47) 

q J 

This expansion of J c g- in momentum space is an equivalent way to understand the approximations used 
in the derivative expansion in Sec. 5.1.1. Note that the definition of the phase stiffness J used here and 
in Sec. 5.1 differ by a factor of 2. The effective model is then discretized in [117] by introducing a lattice 
parameter a: 

H = Ja d ~ 2 ^2(6i - 6 3 f , % - 6j\ < 1 , (5.48) 
where d — 2 is the dimensionality of the system. 

One can also derive an effective XY- model starting from the Ginzburg-Landau theory [203,208,172]. 
As noticed in [117], this procedure which allows one to write parameters of the effective Xy-model 
in terms of the Ginzburg-Landau parameters, is only fully justified close to T C MF where the linearized 
Ginzburg-Landau theory is valid. However, as we saw in Sec. 5.1.2 the most interesting region i.e. that 
with a relatively large pseudogap (T* — Tbkt)/T* ~ 1 corresponds to p/T ~ 1 (or A/T ~ 1) which is 
clearly beyond the limits of validity of the Ginzburg-Landau theory. 

Arguing that the phase stiffness J is only weakly temperature dependent up to a temperature of the 
order of the BCS temperature T C MF it was suggested in [117] to use 

^.^^^AMiM^ (549) 

2 k 

where ^ = ^d^/dk x and £^ — \d 2 i^/ d 2 k x . The low temperature value of J as a function of 1/U is 
shown in Fig. 18. One can see that for large U, J goes as 1/U. 

Assuming a rectangular density of states of width 2W (i.e. a square band extending from — W to W) and 
an electron density n consistent with its definition given by Eq. (3.48), one obtains [116] the following 
zero temperature solution for the system of Eqs. (3.50) and (3.51) (compare with Eq. (3.25) taking into 
account the difference in the definitions of the bandwidth in Sec. 3.1.3 and here) 
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Fig. 18. The low temperature coupling constant J vs the inverse of U, calculated for a square lattice with 
electron density n = 0.2. This figure has been taken from [117]. 

where n is the carrier density per lattice site (see Eq. (3.48)). Note that the chemical potential /i includes 
the normal Hartree-Fock energy shift. This solution is in fact equivalent to (3.25), but we present it 
here for convenience, because it is expressed in terms of often used density n per lattice site for the 
lattice model (3.47). Then, if in contrast to Eq. (3.26), one considers the strong coupling limit W <C U, 
the solution (5.50) can be written as fi = —(1 — n)U/2 and A = \Jn{2 — n)U /2. In the same limit J 
has been obtained in [117] as 



3n(2 -n) 
2 U 



J = 

which implies that J <C A for strong coupling. 

The equation for T B kt was written in [117] for the discrete version of XF-model (5.48): 

Tbkt = AJ , 



(5.51) 



(5.52) 



where A is a dimensionless number of the order of 1 which depends on the details of the short distance 
physics [25]. For the model on the square lattice considered here A ~ 0.89. Note that using the definition 
of J [117] in this section the value of A for the continuum model is n/A ~ 0.79. 

The authors of [117] classified weak and strong coupling regimes by the ratio between J and A. In the 
weak coupling limit at low temperatures, J 3> A and pair breaking, BCS-type excitations dominate 
the low temperature thermodynamical properties. In this regime, the phases exhibit quasi-long-range 
(algebraic) order up to the temperatures of the order of the BCS critical temperature T C MF . It is only very 
close to this temperature, when J rapidly decreases to zero, that the phase becomes really disordered. In 
this regime, the BKT temperature T B kt lies near T C MF , and the transition to the superfluid phase occurs 
concurrently with the appearance of Cooper pairs. In the strong coupling limit, the low temperature 
value of the phase stiffness is such that J< A and the thermodynamic properties at low temperatures 
are dominated by phase fluctuations with Tbkt *C T c mf . In this regime, when the temperature is 
decreased, first Cooper pairs without off-diagonal long-range order are formed and then, below Tbkt, 



the system becomes superfluid. In Fig. 19, Tbkt is shown as a function of U, and for comparison T c 
is shown on the same graph. 
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Fig. 19. The critical temperature, T c of the mean-field (dashed line) and BKT critical temperature, Tbkt 
(full line) vs U, calculated for a square lattice with electron density n = 0.2. This figure has been taken from 



Thus, apart from some details related to the formalism and the behaviour of the system at very low 
carrier densities /i < the phase diagram built in Sec. 5.1 and its physical interpretation for the 
continuum model (3.1) in fact coincides with that for the discrete model [117]. 



5. 3 The effects of Coulomb repulsion and quantum phase fluctuations 



In this section we will discuss the influence of Coulomb repulsion and quantum phase fluctuations on 
the phase diagram built in the previous section. We will follow the papers of Ariosa et al. [118-121]. 
Since the discussion of the pairing mechanism is beyond the scope of the current review, we will not 
consider any effects directly related to the pairing mechanism suggested in [121]. 



5.3.1 Derivation of the quantum XY -model 



[117]. 



Our starting point is the Hamiltonian [120] which describes charge carriers on a lattice, subject to an 
on-site attraction and a long-range Coulomb repulsion, e 2 Vc, acting between particles on different sites: 



H = H + Hu + He , 

H = ~ tc li* C ™* - H £ n " > 

(n,m), cr n 

TT _U_ \ ^ f f 
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^ n,o- 

= 2 Z ( Un ~ n )e 2 ^c(n, m)(n m - n) , (5.53) 

where [/ > 0. Here, as before, (see Hamiltonian (3.47)) cj^ (c ncr ) are creation (annihilation) operators, 
(n, m) denotes denotes pairs of nearest-neighbour sites, n is the neutralizing background charge which 
is of course equal to the carrier density per lattice site, and n n = c Lr c no-- 

The partition function can be written as a functional integral by means of two successive Stratonovich- 
Hubbard transformations [225] (see also [149,226] and compare with Sec. 3.1.1), decoupling the two 
interaction terms in H with the help of a complex field $(n, r) and a real field V(n, r) both of which 
depend on the site n and the imaginary time r: 



Z = Tre-P H 

= J V$V<5>* J w|Tr e -m> Tre -ij; lfl dT(H(*V,T) + e(*y,T)) | ? (5>54) 

with a part linear in the Hubbard-Stratonovich fields 

if V,t) = J2 r)c nT c ni + $(n, r)^ + iV(n, r) (n n (r) - n)} , (5.55) 

n 

and a part containing quadratic terms 

e(<t>, V,r) = ^Yl M™, r)\ 2 + ^ £ V>, r) ^(n, m)V(m, r) . (5.56) 

n n^m 

As in the earlier sections $(n, r) denotes a complex pairing field which couples to the local pair operator 
c n| c nT while the second real field represents a local electric potential which acts on the local charge 
density fluctuations n n (r) — n. One evaluates the trace over the electronic degrees of freedom, with the 
definition : 

Tr |e-^ T r e~*io~ </3 dr[H(*,V,r)+ £ (*,V,r)] | = g -i <M1(*, V,t) (557) 

Then expanding the "free energy" f2 to fourth order in $, to second order in V and to leading terms in 
the space and time gradients of the two fields one arrives at: 

n = ft + ft($,V,T) , (5.58) 



where fig gives the result for the free electron contribution and 
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,EV(n, r) V s - C '(n, m)K(m, r) + 6 £ |#(n, r)| 4 . (5.59) 



Here the Ginzburg-Landau coefficients a, b, c and d are related to the free electron particle-particle 
propagator [53,149,208,172,220,225] and are thus temperature dependent in general; Vgc (n, m) = 
V ( J 1 (n, m)+x (n, m) is the screened Coulomb potential (approximated by its static limit) with Xo( n > m ) 
being the electronic polarizability. 

The terms containing the Ginzburg-Landau coefficients a, b, c and d give the usual time-dependent 
Ginzburg-Landau free energy^, which is usually supposed to be valid only very close to the transition 
temperature equal to T C MF or T p in our case. Thus these coefficients are derived in the same approxima- 
tion as the phase stiffness (5.27) for the BCS limit. This, of course, does not allow one to investigate the 
equation equivalent to the self-consistent equation (5.28) in the regime p/T ~ 1 (see the discussion in 
Sec. 5.2). However, the advantage of the representation (5.59) is its relative simplicity which allows one 
to integrate out the dependence on the electric potential V. Commenting further on Eq. (5.59) note that 
its third term describes the coupling between the two space- and time-fluctuating Hubbard-Stratonovich 
fields; the fifth term in the above equation takes into account the coupling between the local electric 
potential and the number density fluctuations of the charge carriers, while the last double-sum is the 
contribution of the screened Coulomb potential to the thermodynamical potential. 

Integrating over the electric potential V yields 



fi($,r) =^[a|$(n,r)| 2 + 6|$(n,r) 



(n, r) „ ' - 2$(n, r) 
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P( n > r )^sc(n, m)p(m, r) 



(5.60) 



where p(n, r) = 2d|$(n, t)| 2 + (n n ) — n represents the number density fluctuations in the presence of the 
fluctuating pairing field at site n. Roughly speaking, if one interprets \/d&(n, r) as the wave-function 
of a local pair at site n and, thus, its squared modulus as a "pair-density", we note that p(n, r) is 
composed of two separate contributions, namely one from the "unpaired" charge carriers ("normal" 
channel) and one from the "paired" fermions ("superfluid" channel), respectively. Obviously, the charge 
neutrality constraint imposes Yin P( n ; r ) = 0- 

Next the amplitude and phase of $ was introduced in [120] in a way similar to that discussed in 
the previous section. In the following discussion, relatively strong coupling was assumed [120] so that 



Note that the coefficient c is related to the kinetic term which is not obvious from the discrete form (5.59). 



the fermions bind into on-site singlet pairs at a temperature of the order of the mean field transition 
temperature T C MF , well above the superconducting phase transition, the latter being finally triggered by 
the onset of phase order. Below T C MF , a < 0, so that the average amplitude has a non-vanishing mean 
value A = ($(n, r)) given by a + 2bA 2 = 0. Charge neutrality implies 2dA 2 + (n n ) — n = and, for 
t < U [53,149,208,172,220,225]: 

2t 2 1 1 1 

c=— , rf = — , A 2 = -n(2 - n)U 2 « -nf/ 2 for n < 1. (5.61) 
(7 d (7^ 4 2 

Note that, as explained above, the expressions for c, d and the relationship between a, 6 and A are, 
strictly speaking, valid in the vicinity of T C MF while one needs to use them well below this temperature. 
This is the reason why in [120] the proposal was made to use these expressions together with the zero 
temperature expression for A taken from [117] (see Eq. (5.50) and the strong-coupling limit below). 
Splitting the number of pairs at a given site into |$(n, r)| 2 = A 2 + ^U 2 Sn p (n, r), the integration over 
Sn p (n, t) was performed in [120] which (neglecting gradient terms) yields a free energy functional for 
phase fluctuations only: 
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(5.62) 



with the phase stiffness [225,117] J = 2cA 2 = 2nt 2 / £7.0 Note that since the value of J was expressed via 
the zero temperature gap A and neglects the gradient terms for 8n p (l, r), the treatment of the classical 
part of the phase fluctuations becomes equivalent to that in the previous section. The interaction matrix 
W(n, m) in (5.62) is given by 



Vso(n,m), n^m 

Kin, m = < , 5.63 

(2e) 2 ^, n = m. 

The expansion of fi(<&, V, r) in powers of $(n, r) has yielded a spurious on-site repulsion (2e) 2 b/d 2 
between pairs. However, due to the exclusion principle, two pairs cannot sit on the same lattice site. 
Thus, in the following, n = m is excluded in the last term of (5.62); this is equivalent to the usual 
introduction of hardcore bosons. 



The final step in the derivation in [120,226] is a standard Legendre transformation from the "phase 
velocities" dO/dr to the conjugate momenta p(n, r), that are interpreted as the number density fluctu- 
ations of charge carriers, at site n and imaginary time r. The partition function that emerges from the 
latter transformation is equal to the one that would have been obtained from the following Hamiltonian: 



It may also be called Josephson coupling due to the close analogy to the quantum analysis of the Josephson 
junction, see the review [227] and Refs. therein. This analogy is widely used in [118-120,226] for physical 
interpretation. 



H = \Y, (?( n ) - I) [(2e) V(n, m)] (p(m) - | 



+J ]T [1 -cos(0(n) -0(m))]. 

(n,m) 



(5.64) 



The Eq. (5.64) is the central result of the present section: the partition function of the Hubbard model 
with local attraction and long-range Coulomb repulsion has been written as a functional integral with 
a generalized action (or thermodynamical potential) involving a fluctuating pairing field $(n, r) and a 
local electric potential V(n, r). After functional integration over the latter Hubbard-Stratonovich field 
and over fluctuations in |$(n, r)| 2 , the final form of the generalized action involves: 

(i) an effective phase coupling term between the local spatial variations of the phases of $(n, r), which 
presents the same analytical structure as the Josephson coupling term appearing in the Hamiltonian of 
the classical XF-model; 

(ii) a "kinetic energy" term, representing the screened Coulomb interaction between local fluctuations of 
the number density of charge carriers. The Hamiltonian (5.64) represents a non-local (extended) version 
of so-called 2D quantum XY -model which, in addition to the usual Josephson coupling term, involves 
a non-local kinetic energy term; in the case of underdamped Josephson junction arrays, the latter 
describes capacitive (quantum) effects which may arise between the superconducting islands themselves 
but not between the superconducting islands and the normal substrate. 



5.3.2 Mapping to the quantum XY -model with local kinetic energy. Phase diagram 

Here, following [120,121,226] we discuss the application of the Hamiltonian (5.64) for calculating the 
transition temperature of strongly anisotropic superconductors, such as superlattices and bulk systems 
in the underdoped regime. To do this the following approximations are made: 

(i) H is restricted to one superconducting layer; 

(ii) the screened Coulomb interaction, which takes into account the electric coupling between layers, is 
modelled by an analytical form, which is similar to the usual Yukawa potential: 

t/ n e2 ex P(- r /^TF) R _v 

Vsc{r) = —j , 5.65 

ea 2 r 

where e is the dielectric constant of the interlayer material, a is the lattice constant and Atf is Thomas- 
Fermi screening length depending on the density n of charge carriers according to the well-known 
formula 

^-ir^" 1 !- (5 - 66) 

In the last expression ae is the Bohr radius and r s = (3/Ann) 1 ^ 3 is a sphere diameter, related to the 
density of charge carriers n. 

(in) considering only n C 1, the "background shift", —n/2, in the first term of (5.64) is neglected. 



At this point one can make a connection to another, simpler model [118,119] (which is rather similar to 
that proposed by Doniach and Inui [66]), based on the quantum phase fluctuation of the XF-model, in 
which the BKT transition temperature, T B kt is identified with the superconducting critical temperature 
in the absence of charging effects and in which Tbkt is depressed from its classical value by charging 
(quantum) effects [228] (see also [221]). The basic Hamiltonian for this system is the quantum XF-model 
with local kinetic energy: 

#ioc = J E [1 - cos (0(n) - 6>(m))] + ^ £ N 2 n , (5.67) 

(n,m> ^cap n 

where N n is equivalent to the number of "2e-charges" which are present at site n, C cap is the self- 
capacitance of the site. A fundamental property of the quantum XF-model with a local capacitive 
term, is provided by the fact that n n and 9(n) are canonically conjugate variables and therefore satisfy 
the following commutation relation (we restored temporarily the Planck constant h in this relation): 

[N n ,e(m)] = -iH5 n , m . (5.68) 

Thus, the fluctuations of these two are intimately related. In the approach that has been used in 
[118,119,66], the relation between T B kt and real superconducting temperature, which we denote as T c , 
is given by (again we restore the Boltzmann constant fce): 

where (C cap ) is the average effective self-capacitance and a is the ratio between charging and Josephson 
energy. Thus, within the approach [118,119], one only has to evaluate the appropriate (C cap ). It is obvious 
that the quantum phase fluctuations decrease T c with respect to Tbkt evaluated in the presence of the 
classical phase fluctuations only : T B kt is an upper bound to the actual transition temperature [25]. 

It is also important to mention the results of [205] where for the quasi-2D version of the Hamiltonian 
(5.67) it was shown that phase fluctuations produce a linear increase in the penetration depth A(T) with 
increasing T at low temperatures. This in turn implies a linear decrease of the full superfluid density 
(see the discussion in Sec. 5.4.3). 

In order to apply the framework sketched above one must map the Hamiltonian (5.64) onto a "local 
capacitance" Hamiltonian (5.67) using the following approximation: 



2 

where 



Wp(l) [(2e)V(n,m)U(m) ^E?(n) 2 , (5.70) 



2 Ccap i! 



' 4e%0) = ?! T^. (5.71) 



2C cap e 2 „ ea 2 



In Ref. [118] the critical temperature T c for the Hamiltonian density (5.67) has been evaluated in the 
"self-consistent harmonic approximation" (SCHA) which gives results for a <C 1 in good agreement 
with Monte Carlo simulations. A good overall fit of the numerical SCHA result is (see Refs. in [120]) 



T c (a)t*T BKT Jl 

V "c 



(5.72) 



When a approaches a c =6.2, T c goes to zero. 



One has now derived an effective Hamiltonian of the form (5.67), where the capacitance C cap is given 
by Eq. (5.71), and J is the phase stiffness, which is in fact temperature dependent. As seen earlier, this 
temperature dependence is particularly important in the high-density limit where Tbkt — T C MF and 
thus J{T BKT ) ~ J(T C MF ) = J(T = 0). In [121] the temperature dependence is approximated by a 
simple linear form chosen precisely such that J(T C MF ) = i.e. 

J(T) = J(T = 0) (l - , (5.73) 

where J(T = 0) is the zero-temperature phase stiffness. Given this fit to the phase stiffness, one can 
then calculate the critical temperature for the classical phase fluctuations as a function of doping. 
This is nothing but the BKT transition temperature, Tbkt, given by Eq. (5.52) with the value of 
A = 0.9 taken from Monte-Carlo simulations. The critical temperature for the quantum fluctuations 
(i.e. including charging effects) can then also be calculated as a function of doping using Eq. (5.72). 
Note that, it follows from the form of the equations, that both T c and Tbkt are bounded above by T C MF . 
The phase diagram based on only the classical phase fluctuations is shown in Fig. 20. This version of 




Fig. 20. The phase diagram based on the contribution of the classical phase fluctuations. The x-axis shows the 
number of carriers per Cu atom, and here Tbkt has been denoted by Tpn- The non-monotonic behaviour of 
the BCS temperature Tmf is related to the peculiarities of the pairing mechanism. This diagram is taken from 
[121]. 

the phase diagram is in fact very similar to the phase diagram presented in Fig. 15 except for details 
related to the pairing mechanism in Ref. [121] which results in a non-monotonic behaviour for T c MF (n/). 
The idea that, in the above context, the so-called pseudogap be identified with the standard mean field 
superconducting gap is also proposed [121]. Nonetheless, these gaps should not be identical due to the 
quantum corrections, as is the case for the temperatures T C MF and T p or T x from [135] (see the discussion 
in Sec. 5.1.3). 

The phase diagram where the effect of the quantum phase fluctuations is taken into account is shown 
in Fig. 21. The most important effect of these fluctuations is that they push the critical temperature 
to zero in the underdoped region. Thus one can see that the quantum fluctuations lead to a greater 
overlap between the experimental phase diagram in Fig. 1 and its theoretical version in Fig. 21. 
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Fig. 21. The phase diagram from [121] based on both classical and quantum phase fluctuations. Again the 
x-axis shows the number of carriers per Cu atom. 

Furthermore, a possible mechanism which also pushes T c to zero for the overdoped region was suggested 
in [121]. Its discussion is however beyond the phase fluctuation scenario of this review. Finally, to make 
a connection with the next section, we note that another reason why T c is zero for overdoped systems 
is the presence of impurities which can also destroy superconductivity. 



5.4 The effect of non-magnetic impurities 

Up till this point we have considered models for the crossover and the pseudogap behaviour which do 
not take into account the effect of impurities. However, a discussion of the effect of impurities seems to 
be crucial for any realistic model of the HTSC. Indeed, it is known, that the itinerant holes in HTSC 
are created by doping which in turn introduces a considerable disorder into the system, for instance, 
from the random fields of chaotically distributed charged impurities [229]. Thus the purpose of this 
section is to study the influence of non-magnetic impurities on superconducting properties including 
the superfluid density. 

In the theory of "common metals" the Fermi energy €f and the mean transport quasi-particle time 
r tr are independent quantities which are always assumed to satisfy the criterion ei?r tr ^> 1. In HTSC 
which are "bad metals" [26] both ep and r tr are dependent on the doping and the abovementioned 
criterion may fail [229]. In particular it has been shown [229], that for strongly disordered metallic 
systems, superconductivity is absent if the scattering-to-pairing ratio exceeds a critical value. Secondly, 
for such systems [229], even if this ratio is not exceeded, superconductivity only exists for a finite range 
of doping. 

We shall study two cases here. The first is the case of weak disorder originally studied in the papers 
of Anderson [230] and Abrikosov-Gor'kov (AG) [231] (see also [152]). when the superconducting order 
is preexisting and the criterion epT tT ^> 1 is satisfied. In this case superconductivity in disordered 
thin films was first considered by Ramakrishnan [213] using an approach very similar to the method 
presented in the review. The second case is that of strong disorder where this criterion is not satisfied, 
and superconductivity is destroyed by sufficient disorder. 




We concentrate on several important aspects. The first is the influence of impurities on the phase diagram 
particularly at high-densities where the Fermi surface is well-developed. The second is the behaviour of 
superfluid density. Its temperature dependence at low (T <C T c ) and high (T c < T < T*) temperatures 
is not yet understood and presents a challenge for competing theories. One of the consequences of the 
presence of impurities is that the value of the zero temperature superfluid density can be less than 
the total density of carriers, so that their presence may help to explain the experimental results for 
these quantities, [25,26]. The third is the creation of inhomogeneities and the ultimate destruction of 
superconductivity by strong disorder. 

It is interesting to note an alternative explanation for the lowered zero temperature superfluid density as 
a result of quantum phase fluctuations has been recently presented [232] (see also [68]). In the previous 
section the importance of the contribution from the quantum phase fluctuations is clearly evident. It 
follows that the lowering of the superfluid density is the result of at least two factors - namely impurities 
and quantum phase fluctuations. 

5.4-1 Weak disorder limit 

Anderson's theorem [230] (see also [148]) states that in 3D the BCS critical temperature is unchanged in 
the presence of weak disorder in the form of non-magnetic impurities. As discussed in Sec. 5.1, the BCS 
critical temperature (in 2D) is the temperature T* at which the pseudogap opens and is unchanged 
by the presence of impurities. On the other hand the superconducting transition temperature is the 
temperature Tbkt of the BKT transition. In contrast to the former, the latter (as seen in Sec. 5.1) 
is defined in terms of the bare superfluid density or the phase stiffness J (given by the delocalized 
carriers) which decreases with an increasing concentration of impurities. Thus in 2D the temperature 
of the superconducting transition, T B kt decreases with increasing impurity concentration. 

Thus, in the model with impurities, the relative size of the pseudogap phase, (T* — T BK -[)/T* could be 
larger in the presence of impurities [114] than in the clean limit [115] and might therefore be observed 
over a wider range of densities. 



5.4.1.1 Model and main equations Our starting point is, as in Sees. 3.1 and 5.1, a continuum 
version of the 2D attractive Hubbard model defined by the simplest Hamiltonian (3.1) in the presence 
of impurities [152] 



where as before x = r, r denotes the space and imaginary time variables, ip a {z) is a fermion field with 
spin a =|, |, m is the effective fermion mass, fi is the chemical potential, U is an effective local attraction 
constant, and V imp {r) is the static potential of randomly distributed impurities. 

From the Hamiltonian (5.74), following the derivation given in Sec. 5.1.1, one can again derive the 
effective Hamiltonian of the classical XF-model (5.14) with the bare (i.e. unrenormalized by the phase 
fluctuations, but including pair breaking thermal fluctuations) superfluid stiffness: 



H 



+ V imp (r)ipl(x)ip a (x) , 



(5.74) 



J(p,T,p) = E J d 2 ktr[r 3 (g(tuj n ,k))] 
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J <Pkk?tr[(g(iu> n ,k))(G(iUn,k))]- (5-75) 



n=— 00 



Here 



with 



(^(lWn,k)> = 7 2 , 2 s 2 , c2n \ , 5.76 



1 + (5.77) 



is the AG [152] Green's function of neutral fermions averaged over a random distribution of impurities 
and written in the Nambu representation [233,148]. In writing (5.75) we assumed that (Q(iu n , h)Q(iu n , k)) ~ 
(Q(icu n ,\s.))(Q(iu! n ,\s.)). This approximation, as shown by AG [152], does not change the final result for 
J. Note also that the Green's function (5.76) is valid only when e F T tT 3> 1 which demands the presence 
of a well developed Fermi surface, which in turn implies that p ~ ep. Thus one cannot use the expression 
(5.76) in the Bose limit with p < 0. 

Substituting (5.76) into (5.75), and using the inequalities p ^> T,p to extend the limits of integration 
to infinity, one arrives at 

^ 47T J? J -00 l* 2 + (a* + p'Wn I* 2 + lf»l + P 2 )V 2 } 2 ) ' ^ ' ' 

Eq. (5.78) is formally divergent and demands special care due to the fact that one, in contrast to 
Sec. 5.1.1, has to perform the integration over x before the summation [152]. Finally one can formally 
cancel the divergence [152] to obtain 

uo 2 T 00 1 



-°° \ul + P 2 



The temperature of the BKT transition for the Hamiltonian (5.14) is determined by the self-consistent 
equation (5.28). 

The self-consistent calculation of T BKT as a function of the carrier density rif — mep/ir as before requires 
additional equations for p and p, which together with (5.28) form a complete set (see Sec. 5.1.1). 

When the modulus of the order parameter p(x) is treated in the mean field approximation, the equation 
for p (compare with (5.29)) takes in the presence of impurities the following form 

9n 00 r rl 2 h 

77= E j -^ 2 tr[T,mu n ,m (5.80) 

which formally coincides with the gap equation of the BCS theory. This coincidence confirms Anderson's 
theorem [230] which states that in this limit the dependence of p{T) is the same as that for the clean 
superconductor and is not affected by the presence of non-magnetic impurities. 



We recall, however, that there is both physical and mathematical differences between the gap in the 
BCS theory and p (see Sec. 5.1 above). The main point, which we would like to stress here, is that due 
to the Anderson theorem [230] the value of T p does not depend on the presence of impurities at least 
when epT tr ^> 1, while the temperature T B ktj as we will see below, is really lowered. 

The chemical potential \x is defined by the number equation written for the Green's function (5.76) in 
the presence of impurities: 

00 f d 2 k 

/ tr [ r 3 (£?(««>„, k))]exp(i<W 3 ) = n f . (5.81) 

n=-oo J \ Z7T ) 

Since we are interested in the high carrier density region the solution of (3.23) is p, ~ ep, so that in 
Eqs. (5.79) - (5.80) one can replace p by e F . 



5.4.1.2 The effect of weak disorder on the pseudogap region In the clean limit the transport 
time r tr is infinite, so that 

, „2rp OO 1 

j{cMf ' t)) = ^^WT7W- (5 ' 82) 

Near T p one can obtain from (5.82) its asymptotical form (5.27) which has been re-derived here using 
the reverse order for the summation and integration. 

As has been explained in Sec. 5.1.2, if one inserts the dependence of p(T) given by (5.36) into (5.27) and 
then substitutes (5.27) into (5.28) one obtains the asymptotic expression (5.37) for the BKT temperature 
in the clean limit for high carrier densities. 

In the high density limit one can also use the following solution of Eq. (5.34) for T p 

T p = ^J2\e b \e F . (5.83) 

71 v 

(Compare Eq. (5.83) with (3.27) which reproduce together the standard BCS ratio 2A/T p = 2n /j ~ 
3.5.) 

It is obvious from (5.37) and (5.83) that for a clean system the pseudogap region shrinks rapidly for 
high carrier densities. It is therefore natural to ask (see, for example, [107,222] and the discussion in 
Sec. 5.1.2) whether precursor superconducting phase fluctuations can explain the pseudogap anomalies 
which are observed over a wide range of temperatures and carrier densities, since in the clean limit the 
relative size of the pseudogap region (T p —T B kt)/T p is, for instance, less than 1/2 when the dimensionless 
ratio €F/\£b\ ^ 1287 2 /7r 2 ~ 41. If one makes a crude estimate for the dimensionless ratio for optimally 
doped cuprates one obtains ep/\sb\ ~ 3- 10 2 - 10 3 [146] which indicates that in the clean superconductor 
the pseudogap region produced by the phase fluctuations is too small. Of course, all these estimations 
are rather qualitative due to the simplicity of the model. 

We now turn to the dirty limit where the transport time r tr is small (r tr <C p _1 (T = 0)). This inequality is 
not in contradiction with the weak disorder limit but merely implies high density i.e. that ^/Tp » \J\tb\- 
In this limit one can neglect the radical in the bracket of (5.79) [152]. The remaining series is easily 



summed and one obtains the following expression for the bare superfluid stiffness, 



J(e F ,T,p{e F ,T),r tr ) 



4 




(5.84) 



As in the clean limit, due to Anderson's theorem, the expressions (5.36) for p and (5.83) for T p remain 
unchanged in the presence of impurities. Again substituting (5.36) into (5.84) one obtains 



One can see that the size of the pseudogap region is now controlled by the new parameter r tr which is 
an unknown function of e F for HTSC If one assumes that one indeed has weak disorder i.e. e F r tr » 1 
clearly the pseudogap region remains small. If one however moves away from the weak disorder lint and 
uses values for r tr extracted from experiment [4,6], for which r^ 1 ^> p(T = 0) ~ T p then Eq. (5.85) one 
does obtain an increase in the pseudogap region with impurities. Such an estimate is of necessity only 
qualitative. 

Certainly impurities which are always present in HTSC do increase the size of the pseudogap region 
in the weak disorder limit. The effect however appears to be small in this limit and it seems likely 
that the effect of quantum phase fluctuations [232] is more important in allowing one to reproduce the 
experimentally observed extent of the pseudogap region. It would however be interesting to readdress 
this problem in the case of strong disorder. 

The difference in the temperature scales for pairing and phase coherence for higher densities where 
p > is vital to the debate, previously mentioned in Sec. 2.2, as to how the approaches of Uemura [84] 
and Emery-Kivelson [25,26] (see in particular [234]) are related. While Uemura [84] considers that the 
approach based on the phase fluctuations as "essentially identical to the BEC-BCS crossover picture" , 
Emery and Kivelson [234] insist that the separation of the temperature scales for pairing and phase 
coherence is a consequence of the fact that HTSC are doped insulators. Here we would like to present a 
point of view which is closer to that expressed in [234]. Indeed, as we saw in the previous chapters the 
indication that one is on the BEC side of the BEC-BCS is a negative sign for the chemical potential p. In 
Sec. 5.1.2 we saw that in the region with a large pseudogap p was indeed mostly negative phase, so that 
the statement of Uemura [84] applies. However, as we discussed optimally doped cuprates correspond 
to higher carrier densities where p is already positive and p ~ e F . The presence of quantum phase 
fluctuations, and to a lesser extent impurities, however permits one to have a relatively big pseudogap 
region even on the BCS side of the crossover. Another argument in favour of the picture presented by 
Emery and Kivelson is that the value 2A pair /T* ~ 4 — 5 [32] is almost independent of doping while 
it increases on the BEC side of crossover, see Sec. 5.1.2. Thus the BEC-BCS crossover and the present 
approach based on the phase fluctuations need not be related. 

We stress however that this argument about the sign of the chemical potential is not valid in the presence 
of spatial inhomogeneities (see discussion in Sec. 3.1.5) as for example seen in the stripe picture of Emery, 
Kivelson and Zachar [23]. This implies that BEC-BCS crossover and phase fluctuation picture are not 
mutually exclusive since for a spatially inhomogeneous system it is possible to have local regions with 
so small a carrier density that pre-formed Bose pairs do exist there, while in average the whole system is 
still on the BCS side of the crossover. This is, in some sense, similar to the many-band model considered 
in Sec. 3.2 where coexisting local and Cooper pairs were separated not in real, but in momentum space. 





(5.85) 



5.4.1.3 The superfluid density In this section we contrast the experimental results for the bare 
and full superfluid densities with those for the clean superconductor and indicate the likely role played 
by impurities. 

The bare superfluid density is straightforwardly expressed via the bare phase stiffness, n s (T) = 4mJ(T). 
In particular in the clean limit, it follows from (5.82) (or (5.26) with (5.30)) that n s (T = 0) = n/. This 
is not surprising since n s {T = 0) must be equal to the full density n/ for any superfluid ground state 
in a translationally invariant system [235] and the clean system is translationally invariant. We note, 
however, as stated above that in HTSC n s (T = 0) Cn/ [26]. Substituting (5.36) into (5.27) one obtains 
for T close to T p the bare superfluid density 

n.(T-T-) = 2n/(l-r/T p ). (5.86) 

This behaviour of the bare superfluid density is formally the same as the behaviour of the full superfluid 
density in the BCS theory. Nevertheless it is important to remember that the full superfluid density in 
the present model undergoes the Nelson-Kosterlitz jump at T B kt and is zero for T > T B kt- 

We note that one can experimentally probe both the bare superfluid density in high-frequency measure- 
ments [67] and the full superfluid density in low- frequency measurements [236]. The former describes 
mainly the local condensate density distribution which is, of course, nonzero above T c in the precursor 
superconductivity theories, while the latter corresponds to an average condensate density and is zero 
above T c . 

A direct comparison with experiment reveals further discrepancies between experiment and the predic- 
tions given above for the clean limit. The experimental shape of the bare superfluid density (see Fig. 22 
taken from [67]) disagrees with that given by Eq. (5.86). Moreover, one of the most striking features 




T [K] 

Fig. 22. The phase stiffness, G(T) of the Bi2Sr2CaCu20g +( 5 sample with T c = 71K plotted at 140 GHz, open 
squares, and 600 GHz, open circles. The long-dashed line is the bare phase-stiffness temperature. The dotted 
line is the prediction for G in the low- frequency limit, illustrating the universal superfluid jump. This figure is 
taken from [67]. 

is that although n s (T) remains nonzero well above T c , it does not persist throughout the pseudogap 



temperature regime. Instead it falls linearly to zero at a temperature T„ s well below T c . Even though 
the presence of impurities leads one closer to the experimentally observed behaviour of n s (T) in the 
region T c < T < T p (see [114]), the superfluid density still becomes zero only at T = T p . Thus the 
observed behaviour of the bare superfluid density at finite temperature has not yet been fully explained 
in this approach and presents an open question to the current theories of phase fluctuations (for another 
problem related to the superfluid density in the theories based on phase fluctuations see also the end 
of Sec. 5.4.3). 

On the other hand, the zero value of n s (T) above T Hs < T* may be well understood within other 
scenarios which claim that pseudogap is not the result of the superconducting fluctuations but is related, 
for example, to charge- and/or spin-density waves as discussed in [16,34]. In such a case there is still a 
superconducting contribution to the pseudogap even in the region T c < T < T ns . Since T ns is less than 
T*, one can speculate that it could be much easier to explain the observed size of the superconducting 
pseudogap region even within a simple model from Sec. 5.1. Further experimental studies of the doping 
dependence of T Hg and its clear resolution from T c and T* are necessary. 

The experimental data of [236] also shows that the full superfluid density, A/" S (T), does not display the 
Nelson-Kosterlitz jump. This is probably related to the influence of the interlayer coupling as discussed 
in [27]. 

One can however directly obtain a reduced value of the zero temperature bare superfluid density in the 
dirty limit. In this limit the value of the zero temperature superfluid density is given by n s (T = 0) = 
nnfT tr p <C rif since r tr p <C I. This does not contradict the results of [235] because the system is not 
translationally invariant in the presence of impurities. 

We conclude that the low value of the superfluid density in HTSC [26] may at least in part (for other 
mechanisms, for example, based on the quantum phase fluctuations, see [232,68]) be straightforwardly 
related to the impurities which are inevitably present in HTSC. On the other hand its temperature 
dependence remains an important unsolved problem. 



5.4-2 Strong disorder 

As seen in Sec. (5.4.1) impurities are vital for an understanding of the superconducting properties of 
HTSC. We now turn to a discussion of these properties in the complementary case of strong disorder. 
In the presence of strong disorder Anderson's theorem no longer holds, and the pairing amplitude p 
need no longer be constant in space but is in fact spatially dependent i.e. p = p(r). An analysis of this 
effect on the superconducting properties of s-wave superconductors has been presented in [237] and we 
simply represent the main results here. 



5.4.2.1 Model and main equations The model is the 2D attractive Hubbard model (see Eq. (3.47)) 
but with an additional site dependent disorder potential \/„ mp with values chosen randomly from the 
uniform distribution [— V, V] where V controls the strength of the disorder 

H = -t £ (ci a c ma + h.c.) - f/]Tn nT n ni + ]T(K mp - n)n^ , (5.87) 

(nm),<r n n,cr 



Firstly one treats the spatial fluctuations of the modulus of the order parameter in the mean field 
approach, using the standard Bogolyubov-de Gennes (BdG) equations [238] 
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where the site-dependent modulus is given by 

P(r„) =f/^«i(r„K(r n ), 




(5.88) 



(5.89) 



£«i(r n ) = — t J2s u i( T n + <5) + {V^ mp — //)«i(r n ) and similarly for t>i(r n ). Here 5 = ±x, ±y as appropriate 
to a square lattice and p n = p + Un n /2 incorporates a site-dependent Hartree shift dependent on the 
local number density, 



n 



n = J2 l V i( r «) 



The chemical potential, p, is simply given by fixing the average number density, n. 



n 



(5.90) 



where N is the number of sites. These equations are non-linear and have to be solved self-consistently 
for the BdG eigenvalues E\ and the eigenvectors [iii(r n ), i>i(r n )]. In practice they are solved numerically 
on a finite lattice with periodic boundary conditions. 



5.4.2.2 Single-particle spectral gap and bare superfluid stiffness Several random realizations 
for the disorder are considered. For each realization one solves the BdG equations to find p(r n ) and then 
averages over the disorder to obtain a probability distribution P(p) for the local pairing amplitudes. For 
weak disorder, V < 0.25t, P(p) is strongly peaked at its mean-field (BCS) value and one thus recovers 
a spatially uniform p. The distribution broadens significantly for V ~ t and is rather skewed for V ~ 2t 
with a large weight near p = 0. 

One can then also construct the disorder averaged one-particle density of states from the BdG eigen- 
values 

iV M = ^E^-^) ( 5 - 91 ) 

i 

and this has been plotted in [237] for three different disorder strengths V. The most important feature 
of these plots is the fact that there is a finite spectral gap even at high disorder. For each individual 
representation of the disorder, the lowest eigenvalue E h remains non-zero and of the order of the zero 
disorder BCS gap and this is the origin of the spectral gap. This occurs despite the fact that, for high 
disorder, a large fraction of the sites have near zero pairing amplitudes. 

There is however a particularly simple explanation for this behaviour in the high disorder limit as follows. 
Plotting p(r n ) one observes superconducting "islands" , (i.e. clusters of sites where p(r n ) is large) which 
correspond to regions where |V^ mp | is small. These are surrounded by a larger "sea" of sites with p ~ 
and where |V^ mp | is large. The lower lying BdG eigenvalues correspond to the superconducting islands 



and so have finite gap. The BdG eigenvalues corresponding to the "sea" have even higher energies 
because the corresponding disorder potential |V^ mp | is large. 

One notes that in this procedure [237] one calculates the spectrum E\ and N{u) for each individual 
disorder representation and then disorder averages. This will clearly give different results from the 
procedure [239] of constructing a disorder averaged p{vi) and then looking at the energy gap. In the 
latter case one smoothes out the amplitude fluctuations and, in contrast to the picture presented here, 
one obtains a gap which closes with increasing disorder [239]. Note that it is not possible to prove that 
the single-particle strength at u = is exactly zero (particularly for weak coupling U/t « 1) since the 
numerical calculations are by necessity in a finite system. Furthermore, note that the gap will be filled 
by the presence of phase fluctuations, as in the case of zero disorder which will be discussed in Chap. 7. 

One can also calculate the bare superfluid stiffness, D° s) given by [240] 

D° 

_JL = (-k x ) - A xx {q x = 0, g w -> 0, w = 0). (5.92) 

7T 

The diamagnetic term is one-half (in 2D) the kinetic energy (— /C) and the paramagnetic term is the 
(disorder averaged) transverse current-current correlation function. From Fig. (23) one can see a reduc- 
tion in Dg of two orders of magnitude with increasing disorder. It never however becomes zero showing 
that superconductivity is not destroyed by amplitude fluctuations. 

N=12xl2 (n) = 0.875 N=12xl2 (n) = 0.875 
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Fig. 23. The T = BdG bare superfluid stiffness D® and the spectral gap E gap plotted as a function of disorder 
strength for two different values U(> 0) of attraction. Note that the gap persists with increasing disorder while 
stiffness decreases. This figure has been taken from [237]. 

In the strong disorder limit « E gap regardless of the dimensionless coupling U/t. It is thus vital 
to consider phase fluctuations and preliminary calculations in [237] show that it is these fluctuations 
which are responsible for the destruction of superconductivity at high disorder. 




5.4-3 Discussion 



Since in HTSC the pairing scale T* is different from the superconducting transition temperature the role 
of non-magnetic impurities is not traditional and they in fact define the superconducting properties of 
a "bad metal" . In particular even for weak disorder we have seen that these impurities strongly depress 
the superfluid density at zero temperature bringing it closer to that experimentally observed. In the 
case of strong disorder, impurities lead to large inhomogeneities in the system with the existence of 
isolated superconducting islands. This effect then dramatically changes the superconducting properties, 
as evidenced by a greatly reduced bare superfluid stiffness, and ultimately leads to the destruction of 
sup er conduct i vity. 

Of course, the results presented are only qualitative since only the simplest model with non-retarded s- 
wave attraction and an isotropic fermion spectrum was considered. However, there are some indications 
that for d-wave pairing, the properties obtained will persist [2411. P 3 ] 

Further studies of the superfluid density are necessary. For example, it is important to explain the 
concentration dependence of the full superfluid slope, dn s {T)/dT at T = [236,17] As already 
mentioned in Sec. 5.3.2 the classical phase fluctuations-^ result in the linear slope of the full superfluid 
density. Nevertheless, as shown in [243] (see also [39]) the linear dependence can be also produced 
by the nodal (Dirac) quasiparticles with a relativistic dispersion law which are present in a d-wave 
superconductor. It is now widely discussed whether the main contribution to the observed slope of the 
low temperature superfluid density is due to these quasiparticles (see e.g., [244,245]) or phase fluctuations 
[27,246]. It is likely that the dominant contribution to the slope varies with doping as suggested in 
[246]. While in the optimally- and over-doped samples it is produced by the nodal quasiparticles, in the 
underdoped samples it can be compatible with the model of classical phase fluctuations. We stress only 
at the end that the problem of the superfluid slope in HTSC should be considered also in the presence 
of defects of different origin (including dopants). 



6 The Green's function in modulus-phase representation and non-Fermi liquid behaviour 

The main quantity of interest in the present chapter is the one-fermion Green's function G(u;,k). The 
associated spectral function A(u>, k) = — (l/7r)ImG(u; + iO, k), being proportional to the intensity of the 
ARPES [7,8], encodes information about the pseudogap and quasi-particles. Following [125] the Green's 
function for the charged (physical) fermions is given by the convolution (in momentum space) of the 

33 Note that <i-wave superconductivity appears to be completely destroyed by relatively weak disorder if one 
describes the system within the AG theory used in Sec. 5.4.1. However, it is shown in [241] using the finite 
temperature version of the BdG equations (5.88) for d-pairing, that d-wave superconductivity may well survive 
even for strong disorder. This difference results from the more adequate treatment of the spatial fluctuations of 
the complex ordering field modulus p[r) in [241]. One also notes that d-wave superconductivity can be stable 
in the strong coupling limit, or on the BEC side of the BCS-BEC crossover [242]. 

It appears to be more difficult to understand why n s — > as nj —* in the weak coupling model. However, as 
mentioned in the beginning of Sec. 3.1 it is the doped carriers that are weakly coupled and in fact a half- filled 
band in the strong coupling model corresponds here to the empty conduction band. 

35 While the quantum Xy-model was considered in [205] for more exact quantitative estimates, it was shown 
in this paper (see also more recent publication [27]) that the linear slope of the low temperature full superfluid 
density is produced by the classical phase fluctuations. 



propagator for neutral fermions which has a gap p 7^ and the Fourier transform of the phase correlator 
{exp(tT 3 9(x)/2)exp(-ir 3 6(0)/2)). 



Thus the approximation employed here assumes the absence of coupling between the spin and charge 
degrees of freedom. This can however be taken into account at the next stage of approximation. Nonethe- 
less, as we shall see, the phase fluctuations result in non-Fermi liquid behaviour of the system both 
below and above the superconducting transition temperature, T B kt- 



6. 1 The modulus-phase representation for the fermion Green function 



The starting point is a continuum version of the two-dimensional attractive Hubbard model (3.1). As 
discussed in Sec. 3.1.1 the model with the Hamiltonian density (3.1) is equivalent to a model with an 
auxiliary BCS-like pairing field. This equivalent model has a Hamiltonian density, which is given in terms 
of the Nambu variables ty(x) and ^(x) (see Eq. (3.2)) and the auxiliary complex Hubbard-Stratonovich 
fields $(x) = U&(x)t-V(x) = C/V^T and ® ( x ) = U&(x)t+V(x) = Uip\ip\, as 



n = ^ f ( 



X 




/! - T + ®(x) - T^*[x) 



ty(x) + 



\$(x) 
U 



(6.1) 



Let us consider the full fermion Green's function in the Matsubara finite temperature formalism 

G(x) = (V(x)&(0)), x = t,t. (6.2) 

The frequency-momentum representation for (6.2) in the case of the BCS theory in 3D (i.e. in the mean 
field approximation) is given by (3.15). 

However, a problem arises when one tries to apply Eq. (3.15) directly to 2D systems at T ^ 0, since 
a non-zero value for $ = (&{x)) is forbidden by the CMWH theorem [126] (see also the discussion 
in Sees. 2.3). Nevertheless as was discussed in Sec. 5.1 one can assume that the modulus of the order 
parameter p = |$| has a nonzero value, while its phase 6(x) defined by (5.3) is a random quantity. To 
be consistent with (5.3) in Sec. 5.1 the spin-charge variables (5.4), which separate the Nambu spinors 
into their modulus and phase, were simultaneously introduced. 

Applying (5.4) the Green's function (3.15) can be split into spin and charge parts 

Gaf>(x) = £ ^(^)((e lT3eW/2 )^(e- lr3e(0)/2 ) W ), (6.3) 

where 

Q a p{x) = (T Q (x)Tt (0)) (6.4) 

is the Green's function for neutral fermions. We note here that this equation implies that there is no 
interaction between the spin and charge degrees of freedom. Introducing the projectors P± = \{I ± r 3 ) 
one gets 

e ir 3 e/2 = p+el e/2 + p_ e -iff/2 > e ^r,e/2 = p _ e ie/2 + p +e -tf/ 2> (65 ) 



so that Eq. (6.3) may be rewritten as 

G(x)= PaS(x)P p (exp[iae(x)/2]exp[-if3e(0)/2\), (6.6) 

a,/3=± 

where a = (3 and a = —f3 correspond to the diagonal and non-diagonal parts of the Green's function 
respectively. 

For the frequency-momentum representation of (6.6) one has 

00 f d 2 p 

G(iu n ,k)=T W\2 ^ P a Q(iuJ m ,p)P/3D a f3(iuj n -iLj m ,k-p), (6.7) 

m=-oo J \^) a,0=± 



where 



rl/T r 

Q{iuj m ,'k) = J dr J d 2 r exp[iu m r — ikr]Q(r, r) (6.8) 



and 



rl/T r 

D a p{iVL n , q) = J dr J d 2 rexp(in n T -iqr) (exp[ia6(T,r)/2]ex-p[-i{36(0)/2}} . (6.9) 

There is good reason to believe (see Sec. 5.1) that, for T close to Tbkt, fluctuations of the modulus p 
of the order parameter^ are irrelevant. Thus one may safely use for the Green's function (6.8) of the 
neutral fermions the mean-field approximation given by Eq. (5.17) (compare with (3.15)). 

As we shall see, fluctuations of the #-phase transform the pole in the Green's function for the neutral 
fermions (5.17) into a branch cut in the Green's function for the charged particles (6.2) for T < T p . 
The CMWH theorem [126] concerning the absence of spontaneous breaking of a continuous symmetry 
in 2D means that symmetry-violating Green's functions which are described by terms with nonzero 
$ and $* must vanish. However, it says nothing about the gap in the spectrum of excitations, as 
is sometimes incorrectly stated. The correct explanation is that if the symmetry is unbroken, and the 
fermion excitation appears as a pole in the \P two-point function (6.2), then the fermion must be gapless. 
If the fermion does not have the same quantum numbers as ^ (e.g. the neutral fermion T) and so does 
not appear in the \l/ two-point function (6.2) as a one-particle state, then symmetry does not tell us 
whether the fermion (T) will be gapless or not. Thus there is no contradiction between the gapped 
form (5.17) and the CMWH theorem. We note that this formal consideration [210] closely resembles 
the more physical arguments presented by Abrikosov [64] about the possible existence of a gap in the 
one-particle excitation spectrum without the presence of long-range order. 



6.2 The correlation function for the phase fluctuations 



As stated above, phase fluctuations are expected to be responsible for the differences between the 
properties of the charged (physical) and the neutral fermions defined above. The latter are described by 
the Green's function (5.17) which coincides with the BCS Green's function (3.15)) only if one assumes 



These are so-called longitudinal fluctuations which correspond to fluctuations of the carrier density and 
undoubtedly must be taken into account in the very underdoped region. 



that the phase 9 of the order parameter <£> = pexp(i9) is a constant and can be chosen to be zero. This 
is not the case in the 2D model, where there is a decay of phase correlations and the Green's functions 
of charged and neutral fermions are non-trivially related via Eq. (6.7). To establish their relationship 
one needs to know the correlator for phase fluctuations. Its calculation is quite straightforward for 
T < Tbkt, while for T > Tbkt one may apply the results of the theory of the BKT transition (see e.g. 
[133] or the other textbooks mentioned). 



6.2.1 The correlator for T < Tbkt 

In the superconducting phase free vortex excitations are absent and the exponent correlator is easily 
expressed in terms of the Green's function 



D e (x) = (9(x)9(0)) 

(here as above x = r, r) via the Gaussian functional integral 



(6.10) 



D a p(x) = 
V9(x) exp 
1 rVT 



(6.11) 



l/T 



dr\ Id r\ 



-9(x 1 )Dg 1 (x 1 )9(x 1 ) + I(x 1 )9(x 1 ) 



exp 



I rl/T rl/T r r 

~2 Jo dTl Jo d7 ~ 2 J Tx Jd 2 ^ 1 ^!, r^Defa - r 2 , r x - r 2 )/(r 2 , r 2 ) 



where the source is given by 



l( Xl ) = -i^5( n - r)5(n - r) + 1^5(^)5^), (a, (3 = ±) . 



The Green's function 



D e \x) = - J(/i, T, p) V r 2 - K(p, T, p){d T f 



(6.12) 



(6.13) 



for the phase fluctuations has already been found in Sec. 5.1.1 and follows directly from Eq. (5.25). 
Recall that the superfluid stiffness J (see Eq. (5.26)) and the compressibility K (see Eq. (5.22)) are 
both functions of /i, T and p, and also that the Green's function (6.13) includes only the lowest order 
derivatives of the phase 9. Higher terms are also present in the expansion, but will be neglected. 



Substituting (6.13) into (6.12) one gets 

D aj3 {x) = exp 



rp CO 

- y 

4 „f-L 



d 2 q 1 — af3 cos(qr — f2 n r) 



(27T) 



Jq 2 + K^ n 



(6.14) 



It is easy to see that for zero frequency f2 n = the integral in Eq. (6.14) is divergent at q = unless 
a = p. Thus in the sums over a,j3 in Eq. (6.7) only two terms survive, namely 



P-G(iuj n , k)P_ + P + Q{iu n , k)P_ 



+ 



iuj + r 3 £(k) 
^ + £ 2 (k)+p2- 



(6.15) 



It is important that terms like P±G(iui n , k)P T which are proportional to T\ and thus violate the gauge 
symmetry, do not contribute in Eq. (6.7) due to vanishing of the _D + _ and correlators which follow 



them. This explicitly demonstrates that the non-diagonal part of the 2D Green's function is equal to 
zero at all finite temperatures^ in accordance with the CMWH theorem. 



For the non-zero correlators one has 



D(x) = D+Jx) = D(x) 



■ exp 



T ^ f qdqdp 1 — cos(gr cos <p) cos Q n r 
T 2-i 



4 n ^J (2vr)2 Jq' + KQl 



(6.16) 



In the following discussion only the static case r = will be considered. The restriction to this case is 
one of the few main assumptions used throughout this chapter. Recall also another essential assumption 
made in this chapter which is the absence of coupling between the spin and charge degrees of freedom. 

The summation over n and the integration over ip in (6.16) can readily be done yielding the following 
expression for the argument of the exponent in (6.16) 



J dqe- q/A [l - J (qr)} tanh (6.17) 



o 




16nVJK 
where we have introduced the scale 

r «=^jt> (6.18) 
which is a function of the variables used (in the simplest case r av ~ ^/nJT). 

In (6.17) we introduced a cutoff A by hand with the introduction of the corresponding exponential func- 
tion. This cutoff represents the maximal possible momentum in the theory - the Brillouin momentum. 

One can derive from (6.17) (see the details in [125]) the following asymptotic results 

8ttJ r > r av > A" 1 
£>(0,r) ~ l yTavJ T (6.19) 

[ (^)~ S^ 7 , r » A- 1 » r av . 
This long-distance behaviour then governs the ^-fluctuations physics to be studied below. 

Let us discuss the meaning of the value r av . Using again the phase stiffness J(T = 0) and compressibility 
K from Eqs. (5.26) and (5.22) one readily gets that r av = 2y* 'ep/m/T - the single-particle thermal de 

37 Strictly speaking, one can state this, based on Eq. (6.14), only for T < Tbkt- However, there are relatively 
simple arguments [247] why a correlator of the form (exp(i6(x) exp(i#(0)) must be identically zero both below 
and above Tbkt- Indeed, we know that the correlator (exp(i#(x) exp(— i9(0)} (see Eq. 6.23 below) tends to zero 
as r — > oo. On the other hand according to the principle of weakening correlations our correlator must reduce 
to the product (exp(i#(x))(exp(f#(0)) at large r, i.e.the vacuum expectation value (exp(i#)) must be identically 
zero. This of course simply means that one is in the symmetric phase, i.e. one has an exact symmetry with 
respect to the transformation 9 — > 9 + const. In turn this means that all correlators which violate this symmetry 
are identically zero. 



Broglie wavelength (e^ = nnf/m is the Fermi energy). Then, considering that T ~ Tbkt and taking 
T B kt as given by T B kt — £f/8 (see Sec. 5.1.2), one can estimate 



16 16V2 

r av ~ -— = = — — . (6.20) 

The value of the Fermi momentum kp for cuprates is less than the Brillouin momentum A. Thus the 
first line of (6.19) appears to be more relevant. 

There is a second way to estimate r av . One can rewrite it in terms of the value of 2A/T C in the form 

2A 

r av ~ V2n— £ pair , (6.21) 



- c 



where £ pa i r = vp/inA) is the BCS coherence length. (It is assumed, that for the carrier densities of 
interest, fi ~ €f, so that the coherence length £ co h and the pair size £ pa i r are the same. See Sec. 3.1.5). 
This shows that r av has the meaning of a coherence length, which appears to be rather natural since 
the minimal size of the region where there is the phase coherence should be of the order of £o- Since the 
coherence length in cuprates is larger than the lattice space A -1 , one again obtains that the first line 
of (6.19) should be applied. Therefore for T < Tbkt and for static fluctuations we have that 

T 

r 



av 



D(t) = [ — ) 87rJ , (6.22) 



where r av = 16/^/e^m. Note that all these estimations are correct for the clean system. For the dirty 
system one should replace them by the results from Sec. 5.4. 



6. 2. 2 The correlator for T > T BKT 

For T > Tbkt the expression for the static correlator (6.22) can be generalized using well-known results 
from the theory of the BKT transition as [130,131,133]: 

T 



where 



MT)=MT P )expJ J P J (6-24) 

V 1 ~ 1 BKT 

is the BKT coherence length and the value of £+(T p ) to be used in numerical calculations will be 
discussed later. 

Note that Eq. (6.24) represents a very specific temperature dependence of the coherence length which is 
completely different from the temperature dependence of the Ginzburg-Landau coherence length. The 
charge and spin susceptibilities for the negative-?/ Hubbard model (3.47) were estimated in [248] using 
the BKT and Ginzburg-Landau coherence lengths. It appears that the spin susceptibility calculated with 
the BKT coherence length better fits the Monte Carlo data. The fact that the BKT coherence length, 



with its specific temperature dependence, is used allows one to go beyond the T-matrix approximation, 
and to adequately describe the physics of the topological BKT superconducting fluctuations. 



There is also an analogy between the renormalized classical regime discovered in [135], when below 
the temperature Tx (see Sec. 4.5) the two-particle correlation length grows exponentially, and the 
exponential growth of the BKT coherence length below T p . 

One may consider Eq. (6.23) as a general representation for D{r) both for T > T B kt and T < T B kt 
if the coherence length £+(T) is considered to be infinite for T < Tbkt- The pre-exponential factor in 
Eq. (6.23) is related to the longitudinal (spin-wave) phase fluctuations, while the exponent is responsible 
for the transverse (vortex) excitations, which are present only above Tbkt- The pre-exponential factor 
appears to be important for the presence of non-Fermi liquid behaviour which we discuss below. Note, 
however, that the longitudinal phase fluctuations could be suppressed by the Coulomb interaction [123] 
which is not included in the present simple model. One further comment is that while the approximation 
used in [123] to study the vortex fluctuations is good for T well above T B kt, the form of the correlator 
D here is appropriate for T close to T B kt- 

The constant £ + (T p ) may be estimated from the condition that £+(T) cannot be much less than the 
parameter r av which appears in the theory as a natural cutoff, so we will use in our numerical calculations 
C+{T p ) = T av /4. In any case, for T > T B kt where the expression (6.24) is valid, the value £+(T) is large 
and not as sensitive to the initial value of £+(T p ). 

There is also a dynamical generalization of (6.23) 



proposed in [122] from a phenomenological background. Note that t is the real time and T is here the 
decay constant, so that (6.25) is the retarded Green's function. The temperature dependence of T(T) 
for the case of a classical 2D planar magnet was studied in [249] and showed a critical slowing down as 
T — > Tg^x due to the disappearance of mobile, free vortices. 

The more general case of the dynamical phase fluctuations (6.25) was considered in a related calculation 
[122] (see its discussion in Sec. 7.3) where a correlator (exp(i9(t, r) exp(—i9(Q)))\ ss ], which includes 
dynamical phase fluctuations, has been used in the numerical calculation of the self-energy of the 
fermions and the subsequent extraction of the spectral function from the fermion Green's function. 

6.3 The Fourier transform of D(r) 

For the Fourier transform (6.9) of (6.23) one has 




(6.25) 



D(ifl n , q) = J drj d 2 rexp (iQ n r - *qr) (r/r av )- T / 8wJ exp(-r/£+(T)) 

= 2n 5 -^r^ J rdrr^ J J (qr)eM-r/UT)). 
1 Jo 



(6.26) 



It differs from the correlator (6.25) only by a factor 1/2 multiplying the phase. 



The integral in (6.26) can be calculated (see, for example, [250]) leading to 

The hypergeometric function F(a, b; c; z) in (6.27) may be well approximated by a constant since it is 
slowly varying at all values of q. As such a constant we can take the value of hypergeometric function 
at q = oo. Thus, 

£(^,q) = %^ 2 + C 2 Cnr, (6-28) 

where 

47rr(a) / 2 x 2 ^- 1 ) T 

It should be stressed that for T > Tbkt the parameter a quickly deviates from unity as ep decreases; in 
other words, the underdoped region will reveal more non-standard properties than the overdoped one. 

Note that for f+^T) = (T < T BK t) Eq. (6.28) is the exact Fourier transform for the correlator (6.22). 

One should take into account that even for T < Tbkt the propagator (6.28) does not have the canonical 
behaviour ~ 1/q 2 which is typical, for example, for the Bogolyubov mode in dimensions d > 2. In 2D, 
modes with a propagator ~ 1/q 2 , would lead to severe infrared singularities [126], and to avoid this, 
the modes transform into the softer ones (~ l/q 2a , a < 1). 

Finally, substituting (6.15) and (6.28) into (6.7) one obtains 

n( . 1 \ A f d 2 q w w + T 3 £(q) 1 f . 

k) = (2^^ + {'(q) + ^l(k-q)' + GW <6 ' 30) 

The coincidence of the Matsubara frequency in the left and right sides of Eq. (6.30) is straightforwardly 
related to the static approximation. It is truly remarkable that the Green's function (6.30) may be 
evaluated exactly, as may the expressions for the spectral density and for the density of states. We note, 
however, that the static approximation can be justified only a posteriory when the dynamical case is 
also considered [260]. 



6.4 The derivation of the fermion Green's function in Matsubara representation and its analytical 
continuation 



The calculation of the fermion Green's function can proceed analytically along the same lines of calcu- 
lation as in the relativistic case of the Gross- Neveu model at T = [247] One splits the fermion part of 
(6.30) in the following manner 

«*/ + *(k) _ A, + m (631) 



where 



A 1 



1 



r 3 ~ 



U>r, 



Then , using the representations 



A, 



r 3 + 



UJ r 



1 ^ 

— = =p / dsexp [±is(a ± ib)], 

a±ib J 



(6.32) 



(6.33) 



and 



1 1 00 
c Q r(a) 7 



and taking into account (6.31), one may rewrite (6.30) as 



G (^ n) k ) = j^l r rfs r dtr -i e ~c + 2 (T)t- s ^? x 

rial Jo Jo 



I 



d 2 q 



A 1 exp 



A 2 exp 



is- ifxs — (k — q) 2 t 

g 2 

-is h ijis — (k — q) 2 t 

2m 



(6.34) 



(6.35) 



Note that the special form of the integral representation (6.33) (compare with the representation (6.34)) 
guarantees that the Gaussian integral over q is well-defined independently of the sign of £(q) = q 2 /2m — 
/i. Now the Gaussian integration over momenta q in (6.35) can be done explicitly. Subsequently changing 
the variables s — > 2ms and then t — > st one can integrate over s. Finally, making the substitution 
t — > — iu and expanding the quadratic polynomial in the denominator, gives (see the details in [125]): 



G(iu n , k) = 
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(6.36) 



where 



M 2 =?< 2 (T) 



' fc 2 £(T) + l 
. 2m^(T) 



/i + iW ^ + 



(6.37) 



with 



D 



. 2m^(T) 



P 2 ) + _., 2 



m£ 2 (T) 



(6.38) 



and 



u, 



(6.39) 



One can check from (6.37) that Re-Uj > for /i < 0, so that one can rotate the integration contour to 
the real axis: 
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du 
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(6.40) 



Then the integral representation (6.40) may be analytically continued to fi > 0. The change of variable 
z — u/(u + 1) allows Eq.(6.40) to be expressed in terms of Appell's function [251] 
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(6.42) 



For T < T BK t the BKT coherence length is infinite (^ + X (T) = 0) so that the first argument of the 
Appell's function [u\ — l)/ui = 1. This allows one to apply the reduction formula [251] 

F 1 (a, (3, (3>, 7; 1, x) = r(7)r(7 - «-^) ^ pl . i _ p . x) (g 43) 

r(7 - a)r(7-/3) 

and express the result via the hyp ergeo metric function 
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where the value of A from (6.29) has been substituted. 



(6.44) 



6.4-1 The retarded fermion Green's function 



To obtain the expression for spectral density, one firstly needs to obtain the retarded real-time Green's 
function from the temperature Green's function. One does this by means of the analytical continuation 



iu n — > uj + iO, so that ^Ju 2 + p 2 — > iy/uj T ^~p I . This results in the following rules (compare with (6.32), 
(6.37), (6.38)) 

Ai - * =5 + ' A2 A = K T3 - v^M ; (6 - 45) 



« 2 - ^ 2 = m£(T) ( ^[T^ 1 - » ~ ^ 2 ~ P 2 ~ ^) (6-46) 



with 
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and 

Vi = v^u 2 -p 2 - -y^-p 2 ). (6.48) 
Thus for the retarded Green's function one has 
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+ (V^ 2 -P 2 -> -V^ 2 -P 2 )}- (6-49) 

It is also easy to see that 

^2 = -2m£ 2 (T)(/i + ^/uj 2 - p 2 ). (6.50) 

Let us now discuss the conditions under which the imaginary part of G(cj + iO, k) is nonzero. 

For \u\ < p one can see that, v\ = v*, z>2 — v\ so that G(ui, k) is real and lmG(u + iO, k) = 0. The case 
\u\ > p is more complicated. It follows from the Appell's function transformation property [251] 

F 1 (a, (3, P', 7; x, y) = (l- x)- a F, (a, 7 - P - P', P', 7; ^77, j^f) • (6.51) 

that for real x and y the function F 1 becomes complex if x > 1 or/and y > 1. This implies that G(u;, k) 
has an imaginary part if v\ < or/and t>2 < 0. Looking at the expressions (6.46) for v\ and t>2 one can 
see that v\ is always positive, while t>2 may be negative. This means that G(u, k) has an imaginary part 
if V1V2 < 0. Now using (6.50) the condition for existence of nonzero imaginary part of G(u;,k) can be 
written in the following form p + ^u 2 — p 2 > 0. 



6.5 The branch cut structure ofG(u,k) and non-Fermi liquid behaviour 



Let us consider firstly the retarded fermion Green's function (6.44) for T < Tbkt- Applying the rules 
for analytical continuation from the previous subsection to Eq. (6.44) one gets 
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(6.52) 



Near the quasi-particle peaks when uo ~ ±£'(k) = ±y£ 2 (k) + p 2 the arguments of the hypergeometric 
function in (6.52) are close to 1. One can consider, for instance, the first hypergeometric function, so 
that 
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Using the following relation between the hypergeometric functions [251] 
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(6.54) 



one gets that near z\ ~ 1 



G(c,k)~-rW^) x 

A fr(l-2a) r(2a-l) 1 
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It is evident that the expression obtained for the Green's function is non-standard. In addition to 
containing a branch cut, it clearly displays its non-pole character. The latter in its turn corresponds 
to the non-Fermi liquid behaviour of the system as a whole. It must be underlined that the non-Fermi 
liquid peculiarities are strictly related to the charged (i.e. observable) fermions only - the Green's 
function (5.17) of the neutral ones has a typical (pole type) BCS form. It also follows from (6.55) that 
the new properties appear as a consequence of the ^-particle presence (leading to a ^ 1), and because 
the parameter a is a function of T (see (6.29)) non-Fermi liquid behaviour increases with increasing 
temperature and is preserved until p vanishes. 

It is interesting that, in Anderson's theory [252,253] (see also [59]), it was postulated that Fermi liquid 
theory breaks down in the normal state as a result of strong correlations. To introduce this breakdown 



(6.55) 



the following form of the single-particle fermion propagator for the normal state was assumed (see e.g. 
[253,254]): 



G(iJ,k) 
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(6.56) 



where uj c is a frequency cutoff which is introduced to make the Green's function dimensionally correct, 
< a < 1 is independent of T and g(a) is some function of 5 (see more detail in [254] where this form of 
Green's function was used as the starting point for the evaluation of the pair fluctuation propagator and 
the electronic self-energy). It follows from (6.56) that the associated spectral function A(ui,k) satisfies 
close to the Fermi surface a non-Fermi liquid homogeneity relation: 



A(Au,Ak) — A~ a A(u, k) , 



(6.57) 



where the limiting case of a Fermi liquid has a — 1. 

Here we started from Fermi liquid theory and found that it has broken down due to strong phase 
fluctuations, so that the final equation (6.55) reveals the same analytical structure as (6.56). The non- 
Fermi liquid behaviour may, as suggested in [252,59], lead to the suppression of coherent tunnelling 
between layers, that in turn confines carriers in the layers and leads to strong phase fluctuations. 
However, in contrast to [252-254] , the present model predicts the restoration of Fermi liquid behaviour 
as T decreases, since a — > 1 when T — > 0. 



The limit T = can be also obtained in the following way. Strictly speaking one cannot estimate the 
value of r av in the limit T — * in (6.52) via Eq. (6.20) because the substitution of the relationship 
Tbkt — ep/8 into (6.18) is not valid in this case. However, this is not essential because T/SirJ — > 
so that the correlator (6.22), D(r) — > 1 which clearly establishes long-range order in the system. 
Furthermore, the value of a in (6.29) goes to 1 for T — > so that the hypergeometric function in (6.52) 
directly reduces to the geometrical series: 
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(6.58) 



Therefore, substituting (6.58) into (6.52) one gets for the diagonal component Gu(u, k) of the Nambu- 
Gor'kov Green's function G(u, k) the ordinary BCS expression 



Gn(w,k) = 



u + ^k) 
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Clearly Eq. (6.59) results in the standard BCS spectral density [49] with two 5-function peaks 
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(6.60) 



To recover the non-diagonal components of G one has to restore the correlators D_ + (r) and D + _(r) 
that were omitted in Sec. 6.2.1. 



7 The spectral function in the modulus-phase representation and filling of the gap 



As is well known, [49,148,152,184], the spectral features of any system are entirely controlled by its 
spectral density. In terms of the matrix Greens' function (6.49) this is given by 



and, for example, in the cuprates is measured in ARPES experiments [4,7-9] (see also [10], where a 
transparent example shows how interaction, finite experimental resolution and other effects destroy a 
simple S function-like spectral function). The spectral function defines the spectrum anisotropy, the 
presence of a gap, the DOS, etc. Thus it is important to discuss its behaviour for the different models 
based on phase fluctuations. 



7.1 Absence of gap filling for the Green's function calculated for the static phase fluctuations in the 
absence of spin- charge coupling 

In this section we discuss A(u, k) for the Green's function obtained above in Sec. 6. This example is 
rather instructive since it allows one to obtain an analytical expression for the spectral function [125], 
check the sum rule and explicitly show the limitations of the static, non-interacting approximation used 
to calculate the Green's function and its spectral density. 



7.1.1 An analytical expression for the spectral density 

For v± > and v 2 < the retarded fermion Green's function (6.49) can be rewritten (see Appendix B) 
in the following form 
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Then, according to (7.1) the spectral density for the Green's function (7.2) has the form: 
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(7.3) 



Using the quadratic transformation for the hypergeometric function [251] 

2 F l{ a 7 fr t a-b + l;z) = (l- z)'\F x (j,-b+ I + a - b; -Jj^pj ( 7A ) 

the expression (6.29) for A, Eqs. (6.46) and (6.47) one finally obtains 




pa/2 

-(yju 2 - p 2 -> -yju 2 -f?) 



(7.5) 



where the chemical potential // can be, in principle, determined from the equation (5.30) which fixes 
the carrier density. Here, however as stated above, one assumes that the carrier density is sufficiently 
high that p = ep. 

In BCS theory A(u;,k) given by Eq. (6.60) consists of two pieces which are the spectral weights of 
adding and removing a quasi-particle from the system respectively. Note that our splitting of A(u>, k) 
is different since each term in (7.5) corresponds to both the addition and the removal of a fermion. 



7.1.2 The sum rule for the spectral density 

It is well-known that for the exact Green's function G(u, k) the spectral function (7.1) must satisfy the 
sum rule 

/oo 
duA{u,k) = 1. (7.6) 
-oo 

The Green's function (6.49) calculated here is, of course, approximate. One reason for this is the use of 
the long-distance asymptotic result (6.19) for the phase correlator (6.16). This means that its Fourier 
transformation (6.28) is, strictly speaking, valid for small k only, while the expressions have been 
integrated out to infinity. Another approximation that has been made is the restriction to static phase 
fluctuations. Thus it is important to check whether the sum rule (7.6) is satisfied with sufficient accuracy. 

It is remarkable that for (7.5) the sum rule (7.6) can be tested analytically (see the details in [125]). 
One can derive the following result, 

/°° T(a) 
duA(u,k) = —^- (7.7) 
-oo 1 (2 — a) 



The numerical value of the integral at the temperatures of interest may be estimated in the following 
way. At T = T B kt the phase stiffness is given by J = 2/nT BKT , so that for T near to T B kt the value a 
from (6.29) is 



T ~ Tbkt- 



(7.8) 



In particular, a(T = Tbkt) = 31/32 yields the following estimate for the right hand side of (7.7), 
r(a)/r(2 — a) ~ 1.037. This shows that for T ~ T BK t the spectral density (7.5) is reasonably good. 

The parameter a can however differ strongly from unity at T > Tbkt and in the underdoped regime. 

7.1.3 Main results for the spectral density. The problem of gap filling. Limitations of numerical ana- 
lytical continuation 

An example of the plots of the spectral density A(u, k) given by (7.5) at T > T B kt is presented in 
Fig. 24. To draw these plots the value of a from Eq. (7.8) and the mean-field value of p obtained from 
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Fig. 24. The plot of the spectral function A(u, k) as a function of uj in units of the zero temperature gap A 
for k> k F at T = 1.12T B kt as taken from [125]. 

Eq. (5.33) were used. 

The temperature dependence of the spectral function (7.5) is discussed in detail in [125]. It is shown, in 
particular, that below T B kt the quasi-particle peaks at uj = ±E(k) have a finite temperature dependent 
width which is, of course, related to the spin-wave (longitudinal) phase fluctuations. When T — > the 
width goes to zero, but this limit cannot be correctly derived from (7.5) because this is an ordinary 
function, while the BCS spectral density (6.60) is a distribution. This sharpening of the peaks with 
decreasing T in the superconducting state was experimentally observed [8] and represents a striking 
difference from the BCS "pile-up" (6.60) which is present for all T < T c . 

It was pointed out in [123] (see Sec. 7.2) that the broadening of the spectral function caused by these 
fluctuations can be greater than the experimental data permits. 

For T > T B kt the quasi-particle peaks become less pronounced as the temperature increases. Indeed, 
the value of A(u,k) at uj = ±E(k) is, in contrast to the case T < T BK t, already finite. This is caused 
by the fact that P ^ since £+(T) is now finite due to the influence of vortex fluctuations. As the 
temperature increases further, £+(T) decreases so that the quasi-particle peaks disappear (see Fig. 24). 
This behaviour qualitatively reproduces the ARPES studies of the cuprates for the anti-node direction 
[8] (see also [9]) which show that the quasi-particle spectral function broadens dramatically when passing 
from the superconducting to normal state. 

It is important to stress that due to the very smooth dependence of ^ 1 (T) on T (see Eq. (6.24)) as the 



temperature varies from T < Tbkt to T > Tbkt there is no sharp transition at the point T = Tbkt- In 
particular, there is a smooth evolution of the superconducting (excitation) gap A sc = p into the gap 
A PG which is also equal to p and in fact can be called a pseudogap because at T > T BKT the system is 
not superconducting. This qualitatively fits experiment [4,8,9,31] and appears to be completely different 
from BCS theory [49], where the gap vanishes at T = T c . 

We, however, would like to stress the main reason why further studies are essential. While the temper- 
ature behaviour of the quasi-particle peaks is quite similar to experiment [4,8], the gap in the spectrum 
remains unfilled, i.e. the spectral density is identically zero inside the gap (A(u;,k) = for \uj\ < p). 
Furthermore, there is also an excess of the spectral weight on the gap edges which is seen as the extra 
peaks [125]. This is obviously related to the facts that the classical phase fluctuations were treated in 
the static approximation and that the coupling between the spin and charge degrees of freedom was 
neglected. One therefore needs to consider how these effects could lead to the filling of the gap. 

The T- matrix approximation seems to result in the opposite result. For example, the classical fluctua- 
tions do lead to a pseudogap, as shown by the Hubbard model at half-filling [135]. We note, however, 
that our statement is applicable to the uncoupled static phase fluctuations only and neglects entirely 
modulus fluctuations. Further studies are necessary to make a comparison of this analytical result for 
the phase fluctuations with other complementary approaches, which we believe should reproduce our 
result when the same assumptions are made. Here we make only one further comment related to the 
numerical analytical continuation. 

It is very instructive to estimate the accuracy one would expect if a similar problem were to be in- 
vestigated numerically. Using the analytical expression (6.42) obtained for the Green's function on the 
imaginary axis one can perform its numerical analytical continuation by means of Pade approximants 
[194] (see also [183] and a recent comprehensive discussion of the method in [195]). Using Eq. (6.42) cal- 
culated to ten digits of precision we tried to recover the known spectral function (7.5) shown in Fig. 24, 
but even for 110 Matsubara frequencies the numerical result presented in Fig. 25 cannot indicate whether 
the gap is filled or not. As already mentioned in Sec. 4.1, numerical analytical continuation is often used 



A(co, k > k F ) 








-3-2-10 1 2 

cj/A 



Fig. 25. Plot of the spectral function A(u, k) as a function of oj in units of the zero temperature gap A for 
k > kF at T = 1.12Tbkt obtained by means of numerical analytical continuation of Eq. (6.42). 

for the subsequent extraction of the spectral function in the methods where the T-matrix equations 
are solved numerically on the imaginary axis. It is clear, however, that the precision of the numerical 
solution is even less than that used above for Eq. (6.42) and thus numerical results alone are insufficient 



for definitive conclusions about pseudogap filling. There is however a real-time technique [99,100,104] 
which allows one to avoid numerical analytical continuation, but demands some other assumptions. 

The results discussed in this section show that, while the physical picture based on non-interacting static 
phase fluctuations is sufficient to derive at least the qualitative phase diagram, additional consideration 
is necessary to understand the mechanisms responsible for gap filling. Thus the next two sections will 
be devoted to two possible ways that the pseudogap can be filled. 

7.2 Gap filling by static phase fluctuations due to quasi-particle vortex interactions. The phenomenol- 
ogy of ARPES 

As shown by Franz and Millis [123] (see also [124]) one way of filling the gap can be obtained for static 
phase fluctuations via the coupling of the spin and charge degrees of freedom which leads to a Doppler 
shift for the quasi-particles. 

As already discussed in Sec. 6.2, for this mechanism the disordered state above Tbkt can be thought of 
as a "soup" of fluctuating vortices with positive and negative topological charges and with total vorticity 
constrained to zero. Each of these vortices is surrounded by a circulating supercurrent which decays as 
1/r with the distance from the core. Such supercurrents, within a semiclassical approximation, lead to 
a Doppler-shifted local quasi-particle excitation spectrum of the form [238,255] 

£ k = 4 0) + ftk-v s (r), (7.9) 



where v s (r) is the local superfluid velocity and E^' = ye k + |A k | 2 is the usual BCS spectrum. The 
change in the local excitation spectrum will affect the spectral properties of the superconductor in that 
the physically relevant spectral function must be averaged over the positions of fluctuating vortices. 

This effect will be particularly pronounced in a rf-wave superconductor since Eq. (7.9) implies the 
formation of a region around a nodal point on the Fermi surface with E^ < for arbitrarily small v s (r). 
Physically this corresponds to a region of gapless excitations on the Fermi surface which leads to a finite 
density of states (DOS) at the Fermi level. As first discussed by Volovik [256], a similar situation arises 
in the mixed state of a rf-wave superconductor where the superflow around the field-induced vortices 
leads to a residual DOS proportional to \TR. 

7.2.1 Phenomenological and phase fluctuation coupled Green's functions 

In the mean field approximation (neglecting, among other things, phase fluctuations) the diagonal 
Green's function was written in [123] (compare with Eq. (5.17)) 

A 2 

g 1 (u,k) = u-e k + iT 1 (7.10) 

to + e k 

where the angular dependence of the gap function for rf-wave superconductor is A k = A^ cos 26* and 
following [257] the single particle scattering rate r x was added to the usual mean-field solution to 
describe the ARPES data in the overdoped samples where there are no strong pseudogap effects. It is 
instructive to use Dyson's equation (4.3) to extract from (7.10) the expression for self-energy as done 



in [257,9] 



S(o;,k) = -ir 1 + A * . (7.11) 
uj + tO + e k 

Note that the form of the scattering rate in Q (E) constitutes a non-trivial assumption. It is not pair- 
breaking, in the sense that it is ineffective at small uj and ek; i.e. in the region uj < By contrast in 
a <i-wave superconductor, a conventional scattering rate enters via the replacement uj — > uj + iT, leading 
to a broadening which is effective even at low uj and ek. As shown by Norman et al. [257] the form given 
in Eqs. (7.10) and (7.11) agrees with the ARPES data at T < T c . It is demonstrated in [123] that it 
also agrees with STS. 

At T > T c Norman et al. [257] showed that additional pair-breaking scattering is needed to account 
for the ARPES data in the underdoped case where the gap is filled well above T c , They modelled 
this phenomenologically by introducing another scattering rate F ^ F 1 , and making the replacement 
uj — > uj + iT Q in the last term of Eqs. (7.10) and (7.11): 

E(w, k) = -iT, + A ^ . (7.12) 
uj + e k + il 

They suggested that To could arise from exchange of pair fluctuations and precisely this kind of effect 
has been intensively studied within T-matrix approximation (see Sec. 4.1), but as argued in [123,135] 
this proposed mechanism does not account for the observed magnitude of To- 

As mentioned above another likely source of the pair-breaking scattering leading to gap filling has been 
suggested in [123], namely supercurrents induced by phase fluctuations. In order to determine how Q is 
changed in the presence of superflow it is useful to recall the origin of the energy shift in Eq. (7.9). This 
can be derived [238,255] by assuming a state of uniform superflow with v s = hq/m induced by an order 
parameter of the form e 2iqr Ak. By solving the appropriate set of Bogolyubov-de Gennes equations and 
retaining only terms to linear order in q, one finds that the energy is modified as indicated in (7.9) while 
the coherence factors are to the same order unchanged. This result is then semiclassically extended to 
non-uniform situations by assuming slow spatial variations of v s (r). 

One can follow this exact procedure and solve the appropriate Gor'kov equations for (? q in the presence 
of superflow. One finds the following intuitively plausible result which is exact for uniform flow up to 
terms linear in q: 

g n (uj,k) =0 o (w-77,k-q), (7.13) 

where rj = ftvp(k) q ~ Kk-v s . Here v^(k) = (<9ek/<9k) fe=fcF ~ ftk F /m is the Fermi velocity and the last 
equality holds when the Fermi surface is approximately isotropic. In the following it is assumed that 
Eq. (7.13) can be applied locally when v s (r) varies slowly in space. Applying the above prescription to 
(7.10) one finds, again to the leading order in q, 

Q-\u, k) = uj - e k + iT, - (Ak " C) 9 2 , (7.14) 

uj + e k — 2,7] 

where ( = v A (k) -q with v A (k) = ((9A k /i9k) fe=fcF . One can easily estimate va/v f ~ (Ccoh&F)^ 1 ~ A rf /e F 
which is typically a small number in superconductor. One therefore expects that ( rj. A more detailed 
numerical analysis indeed shows that, as long as A^/ei? is small compared to unity, the effect of ( on 
the spectral lineshape is negligible compared to that of r/, and will be dropped in the following. 



A typical experimentally measured quantity, such as the ARPES or STS lineshape, will provide infor- 
mation on the spectral function associated with £ q averaged over the phase fluctuations (see also [258]). 
Thus, one needs to evaluate 



<? q (w, k) = J dr}P(rj)Gq(u, k) , 



(7.15) 



where P is the probability distribution of r] given by 



P(r ] ) = (5[r ] -hk-v s (r)]). 



(7.16) 



The angular brackets indicate thermodynamic averaging over the phase fluctuations in the ensemble 
specified by the 2D XY Hamiltonian (5.14). 

Further details of the precise derivation and the comparison with the ARPES and STS experimental 
data can be found in [123], but will not be discussed here since our main goal is simply to discuss 
possible mechanisms for gap filling. We note only that since the calculations reported in [123], new and 
controversial ARPES data has been obtained [10,11]. 

As already discussed the last term in Eq. (7.10) can be thought of as a superconducting self-energy given 
by Eq. (7.11). Eqs. (7.14) and (7.15) then imply that the primary effect of the phase fluctuations is to 
smear the functional dependence of £(a>, k) on the energy variable, broadening the spectral lineshape. 
The detailed analysis performed in [123] shows that n acts primarily to fill in the gap, in a way similar to 
the inverse pair lifetime To introduced by phenomenological considerations in Ref.[257]. IY on the other 
hand, does not affect the lineshape at low energies: notice that since the self-energy given by Eq. (7.11) 
diverges at the point k = kp, u = if A k ^ 0, the corresponding Green's function Q (k F ,uj = 0) = 
for any IY 

7.2.2 Discussion 

The qualitative behaviour of ARPES and STS lineshapes (see [4,8,31,257]) in underdoped BiSCCO 
clearly establishes the existence of a scattering mechanism which becomes operative at T > T c and 
which acts primarily to fill in the gap at low energies. 

It is shown by Franz and Millis [123] that vortex phase fluctuations associated with proliferation of 
unbound vortex-antivortex pairs in the system interacting with the quasi-particles provide a reasonable 
explanation for this scattering. 

Analysis of [123] also indicates that longitudinal (spin wave) phase fluctuations are almost completely 
suppressed, above and below T c . As we mentioned in Sec. 5.3.2 (see also Sec. 5.4.3) it has been proposed 
in [205] that in high-T c materials, longitudinal phase fluctuations governed by the XY Hamiltonian 
(5.14) are important in that they significantly contribute to the observed temperature dependence of 
the magnetic penetration depth. The broadening of the spectral function which would be caused by 
these fluctuations was calculated in [123] , and it was found that it is much greater than the experimental 
data would permit. They therefore conclude that longitudinal fluctuations are suppressed, perhaps by 
the Coulomb interaction. 

Quantitatively there exists considerable discrepancy between the parameters describing the ARPES and 
STS lineshapes, in particular the single particle scattering rate IY It is possible that this discrepancy 
can be resolved using new ARPES data [10,11]. 



Sizable transverse phase fluctuations implied by this work will also affect other properties of the un- 
derdoped systems, such as the electronic specific heat, fluctuation diamagnetism and transport [258]. 
Vortices existing above T c should also generate local magnetic fields which are zero on average but have 
a non vanishing variance. If such fields could be detected, e.g. by muon spin rotation experiment, this 
would constitute direct evidence for the phase fluctuation model of the pseudogap phase. 

Finally we mention the recent paper [259] where the effect of phase fluctuations on inelastic neutron 
scattering and NMR experiments was studied using an approach similar to that of [123]. 



7. 3 Gap filling due to dynamical phase fluctuations without quasi-particle vortex interactions 



Another mechanism of the pseudogap filling due to vortex-vortex dynamics has been suggested by 
Capezzali and Beck [122] where the same s-wave pairing model as in Sec. 6 was considered. This model 
is treated by the Stratonovich-Hubbard transformation (see Sec. 3.1.1), decoupling the interaction term 
by the complex pairing field $. The one-electron Green's function is then approximately given by 
[97,226]: 

-j ( $ ) 2 
G(iu n ,k) = iu n - e(k) + n - a(iu n ,k) - - ■ — — ■ — - — : — - . (7.17) 

iuj n + e(k) - pL + a(-iuj n , k) 
The expression for the self-energy in the next to leading order approximation 

00 r d 2 q 

a(iLU n ,k) = - ]T / — d -(mia m ,q)\ 2 )G(iu; n -iu m ,k-q) (7.18) 

m=-oo J \ Z7T ) 

involves the dynamic correlation function of the pairing field. Note that Eq. (7.18) does not contain 
the leading order of the self-energy which is included into the Green's function (7.17) via the term 
with pairing field (<&) (compare with Eqs. (7.11) and (7.12)). In diagrammatic language Eq. (7.18) 
corresponds to a two- vertex bubble, one of the lines of which represents the propagator of the complex 
bosonic Hubbard-Stratonovich field, while the other line corresponds to the dressed propagator of the 
fermions. 

In [122] rather than aiming at a self-consistent solution as described in Sec. 4, a simple form for 
the pairing correlation function was assumed and its influence on the one-electron properties were 
investigated. Introducing the amplitude and phase in the pairing field (see Eq. (5.3)) but without the 
corresponding replacement for the Fermi field (see Eq. (5.4), Capezzali and Beck made the following 
assumptions: 

(i) The mean-field transition temperature T C MF is, following [32], identified with the temperature T*. 
Below this temperature the amplitude (modulus) fluctuations are assumed to be space- and time- 
independent, i.e. they are treated in the mean-field approximation in a way similar to that discussed in 
Sees. 5.1 and 6.1. They are thus approximated by a typical mean-field (BCS-like) analytical form (see 
Eq. (5.36)), 

(|$(t,r)| 2 )=A 2 (T)^A^l-^ (7.19) 

where Aq is treated as a free parameter. For BCS theory it is given by A 2 , = 3.02A 2 (T = 0). Note that 
once more the average ($(£, r)) is zero. Using the definition of the Fourier transform (3.13), and space 



and time independence of the amplitude, one can obtain that 



(|$(^ m ,q) 2 |> = / dr /rf 2 rexp(^ m r-2qr)($*(r,r)$(0,0)) 

p Jo J (7.20) 
= A 2 (T) J dr J d 2 r exp(^ m r - iqr) (exp(^(r, r) - 6(0, 0))) . 

The Fourier transform of the phase correlator in (7.20) was taken in the form given by Eq. (6.25) (see the 
details in Sec. 6.2, in particular the last paragraph of the section). Thus the main difference between the 
first assumption of [122] and the assumptions used in Chap. 6 is the consideration of the more general 
case of the dynamical phase fluctuations. As will soon become apparent it is their presence that leads 
to gap filling. 

(ii) Below the critical temperature T c one keeps the same form (6.25) for the phase correlations, but with 
^(T) = 0, corresponding to the algebraic decay of correlations. However, in contrast to the calculations 
in Chap. 6, a non-zero coupling between the planes (in the third dimension) is taken into account in 
[122]. This is achieved by introducing a non-zero value for the average of 3>(x), ($(a;)) 2 = A(|$(a;)| 2 ) 
with a variable parameter A < 1. Due to the presence of the third direction the nonzero value of (<&) is 
not forbidden by the CMWH [126] theorem. 



7.3.1 Discussion 

For numerical computation the value of a(iu n , k) from (7.18) was evaluated in [122,226] to lowest order 
by using the noninteracting G(ia> n ,k) with an isotropic spectrum e(k) = k 2 /2m. 

In Fig. 26 we show the spectral function from [122] obtained for the decay constant (see Eq. 6.25) 
T = 0.5e,F and the modulus strength (see Eq. 7.19) A 2 , = 0.76ef7. The following observations can be 
made on the basis of this figure and the other results presented in [122] : 

(i) Above T c , a pseudogap opens around ji ~ e F . Its effective width is almost T-independent, in spite 
of the T-dependence of (l^l 2 ) in the self-energy. 

Below T* , but well above T c , one sees a one-peak structure in the spectral function for all wave- vectors 
k, centered below (above) the chemical potential for k < (>)k F - As one approaches T c , a two peak 
structure emerges for wavevectors ~ k F where the peak below (above) the chemical potential has a 
greater spectral weight for k < (>)k F . This structure becomes increasingly visible near T c where the 
density of states at the chemical potential, N (//), becomes practically zero for all k ^ k F . Thus in this 
limit the pseudogap is delimited by two rather pronounced peaks. 

It is important that, in contrast to the previous section, the filling is achieved without any vortex-quasi- 
particle interactions. This, of course, does not exclude the fact that there are phase fluctuations. If one 
could divide the total contribution of the fluctuations in (7.20) into the classical (m = 0) and other 
(in some sense non-classical with rn^O) fluctuations, it can be shown that the latter contribution also 
results in some gap filling [260]. 

(ii) Below T c , the presence of ($) ^ produces two new peaks: a true "superconducting gap" emerges 
out of the pseudogap. Its width is given by the geometrical superposition of average and fluctuating 
part of $(x). The fluctuations of the latter remain visible in the form of secondary shoulders inside the 
superconducting gap which approach each other about in the same proportion as the main peaks move 



Fig. 26. The spectral function A{oj,\l) taken from [122] for k = 0.5/cf (dashed line), k = 0.7kp (dotted line), 
k = kp (full line), k = \.2kp (dash-dotted line), k = l.5kp (heavy full line), for T = 1.01T C and (inset) 
T = 1.5T C . 



away from each other. 

(Hi) The value for the decay constant T ~ 0.5/x is relatively large. For a smaller T, and in particular 
for the case of critical slowing down of phase fluctuations (T (T c + ) = 0) [249] the two shoulders would 
become secondary peaks inside the superconducting gap. 

(iv) Near T c , as seen above, the spectral functions are doubly peaked in some wave- vector domain around 



k,F- In this case the width of the pseudogap is given by the separation between these two peaks, which is 
essentially determined by (|$| 2 ). At higher temperatures, approaching T*, the spectral functions have 
only one peak. However their width is enhanced over essentially the same /c-domain. Due to this fact, 
the pseudogap fills up gradually, as T* is approached (from below), but with little change in its width. 

Thus the mechanism for gap filling due to vortex- vortex interactions which are described by dynamical 
phase fluctuations can be also regarded as a strong candidate. Both mechanisms are based on the same 
physical premise, namely the existence of a "soup" of vortices above T B kt- Indeed both mechanisms 
ought to be present simultaneously, but it is not yet clear which is the dominant mechanism for gap 
filling, and more studies are certainly needed. 

On the other hand, it must be pointed out that the existence of fluctuating vortices is not a necessary 
condition for the appearance of a pseudogap. Indeed, systems with SO(n > 2) symmetry (where n is 
the number of the order parameter components) do not have topologically stable vortices, but still may 
possess a pseudogap- like regime. As discussed in [261] there is a situation when the second order phase 
transition in such a system is driven by long- wavelength transverse (or directional in terms of [261]) 
fluctuations, and not by longitudinal (or size) fluctuations as happens in more conventional scenarios. 
During this sort of transition the ordering field acquires firstly as the temperature decreases its size 
(modulus) while the direction is still strongly fluctuating and only a further temperature decrease 
leads to the formation of the ordered state. Thus one can concede this scenario of a second order 
phase transition leaves room for pseudogap-like behaviour even in more general theories where the 
stable vortex excitations are absent. Nevertheless, their presence in the particular case of 5*0(2) theory 
considered here suggests that vortices can play an essential role and their specific contribution may even 
be dominant in this case. 



8 Concluding remarks 

We have concluded the description of the physical properties of 2D metals in which conductivity is the 
result of doping. We paid particular attention to the normal and superconducting properties directly 
related to phase fluctuations of the complex ordering field which increase significantly as the carrier 
density decreases. The vital role of phase fluctuations for phase transitions in low dimensional systems, 
as mentioned above, has been clearly outlined in the classical papers of Berezinskii, Kosterlitz and 
Thouless and the many other physicists who followed. However for many years physicists have ignored 
the important fact that even above T c (Tbkt), in the so called normal phase, there is a remaining non- 
zero parameter, which drastically changes the spectrum of the phase and leads to pseudogap behaviour. 
While in theoretical papers this feature was simply overlooked, experimental techniques probably did not 
have sufficient precision to observe the pseudogap. During the last 3-5 years the situation has drastically 
changed and the resolution of ARPES and thermodynamical experiments permit one to measure quite 
subtle spectral features, and thus leading to the discovery of the pseudogap. According to Abrikosov 
"the pseudogap in the quasi-particle spectrum above T c and its amazing stability in strongly underdoped 
samples was the most spectacular among phenomena" observed in high-T c compounds [64] . The theory 
discussed above, however, does not relate the pseudogap appearance directly to the superconducting 
properties of oxides and stresses the possible main role of the low dimensionality of their magnetic 
and electron spectra. And indeed, the experimental data is not in contradiction with this conclusion. 
Pseudogap features have already been found in quasi-lD superconductors [262,263] as well as in the 
HTSC systems mentioned above. It is interesting that pseudogap behaviour close to the Fermi level 



is also seen in 3D manganite systems [264]. However, in the last case its connection with some kind 
of fluctuation is under question [265]. Nevertheless even if the pseudogap in manganites is caused by 
different effects, it is remarkable, that in all these cases the systems are bad metals, i.e. for the systems 
under consideration, the undoped state is an insulator and conductivity only results from doping. 

Throughout this review we discussed mainly the simplest 2D fermi-system with attraction. There is, 
however, a general theoretical question as to what are the generic effects that lead to a pseudogap. It 
is clear that to gain greater insight into this problem a wider class of the models should be considered. 
Nevertheless, some general arguments can be given already. First of all, the lowered dimensionality 
of the space is an essential component. This results in the CMWH theorem that allows one to state 
that there is a new region between the mean-field transition temperature T C MF and a true transition 
temperature which is in d = 2 for the systems with two-component (SO (2)) order parameter may be 
well approximated by the of topological transition temperature, Tbkt- This case of SO (2) symmetry, 
which is relevant for superfluids and superconductors, led us to the rich BKT physics of vortices that 
undoubtedly contributes to the pseudogap properties of these systems. However, as we saw even in 
these 2D systems the pseudogap region appears to be small if the carrier density is too high, or in other 
words the coherence length is too big and size fluctuations dominate over directional [261]. In d = 3 
this coherence length should be even smaller [261] and therefore the conditions for the pseudogap in 3D 
is difficult if not impossible to satisfy. 

What would happen if one considers a multi-component (SO(n > 2)) order parameter? Since there is no 
stable vortices in this case, the topological transition is also absent. Thus as mentioned, for example, in 
[112] the higher symmetry of the order parameter space leads to an even larger region with directional 
fluctuations than in the case of SO (2) symmetry. Therefore the larger dimensionality of the space of 
the order parameter also favours pseudogap formation. Further study of these types of models would 
be helpful for a deeper understanding of the pseudogap origin. 

A series of experiments, in particular, NMR, neutron scattering and spin relaxation demand the intro- 
duction of the concept of a spin gap. This is not surprising because the undoped magnetic subsystems, 
which are antiferromagnetic insulators have spins localized on the copper ions and are characterized 
by the magnetic anisotropy of "easy plane" type. This shows that the average magnetic moment of 
every site is nonzero, and therefore the average antiferromagnetic vector in the basic CuC>2 plane is 
also nonzero.^] Doping which suppresses the long-range magnetic order, does not, however, suppress 
the average value of the modulus of the spin on the lattice point. This spin can also be parametrized 
via its modulus and phase. Then the physical picture in the spin spectrum will be very similar to that 
for superconductivity, when the gap in the spectrum is related to the average spin (modulus of the 
magnetisation) and the absence of the long-range order is the result of "disorder" by angular (direc- 
tional) fluctuations which can also be described by a correlation function. Thus the system may have 
an additional crossover temperature related to modulus formation for another - spin - physical value. 
The consideration of this system, which will be close to that described above, is given either by the 2D 
sigma model or by the model of nearly antiferromagnetic liquid. We did not however describe them in 
detail because our goal was to present simple analytical (but at the same time, we hope, non-trivial) 
models which can lead to pseudogap character in the one-particle spectrum. The spin aspect inevitably 
demands in the absence of long-range order a separate treatment and the explanation is therefore less 
clear, and was avoided for simplicity. 



39 It is remarkable that in the paper [266] the role of these planes in the physics of HTSC is compared with the 
role of the hydrogen atom in atomic physics. 



Another important and new feature we wanted to stress is the role of impurities. All bad metals are 
such that each carrier has its "own" dopant. Thus it is necessary to show that impurities do not lead to 
pseudogap suppression. It is clear, however, that this direction needs substantial development to include 
more fully the physically relevant case of strong disorder and the effects of localization. As seen above, 
these effects are probably in fact responsible for the development of an inhomogeneous state. 

Finally, it is very important to develop the theory for the case of anisotropic pairing. Although, as we 
saw, the BCS-Bose crossover theory has already been generalized to the case of (i-wave pairing, the 
formation of an anisotropic pseudogap and an anisotropic BKT transition remain serious problems to 
be addressed.^ Unfortunately, the analytical treatment of the problem is very limited and, as we saw 
in Sec. 7.1.3, numerical methods have their own limitations. 

Therefore, there is no doubt that physicists studying pseudogap properties have many problems to solve 
and the results described above are only one step in this process. 



A Calculation of the effective potential 



Here we sketch the derivation of the effective potential. To obtain it one must write Eq. (5.12) in the 
momentum representation: 
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where 
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are the inverse Green's functions. Note that the derivation of the effective potential (3.16) in the "old" 
$ and <£>* variables differs only in the replacement of the T\p term by r + $ + r_$*. 

The exponential factor e iSuj " T3 is added to (A.l) (see also Eqs. (3.63), (4.6), (4.11), (4.21) and (5.81)) 
to provide the correct regularization which is necessary to perform the calculation with the Green's 
functions [152]. For instance, one obtains 



40 Here we refer a reader to a very recent paper [267] which exploits a discrete version of the variable transfor- 
mation (5.3), (5.4) to study the low temperature properties even for d-wave pairing. 
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To calculate the sum in (A.3), one must first use the identity TrlnA = lndet A, so that (A.l) takes the 
form 
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Calculating the determinants of the Green's functions (A. 2) one obtains 
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where the role of G (iui n , k) in the regularization of f2 pot is now evident . The summation in (A.5) can 
be done if one uses the representation 
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Integrating (A. 8) over x, one thus obtains 
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Finally, substituting (A. 9) into (A. 5), 
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It is easy to show that at T = 0, the expression (A. 10) reduces to that obtained in Sec. 3.1.2. 
Finally, we give formulas for the summation over the Matsubara frequencies used in Sees. 5.1.1: 
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where the Green's function Q(iu n , k) is given by (A. 2). Both formulas can easily be calculated using 
Eq. (A. 7) and its derivative with respect to c. 



B Another representation for the retarded Green's function 



Here we shall obtain another representation for the retarded fermion Green's function which is more 
convenient for the derivation of the spectral density. Recall that when the imaginary part of G(u, k) is 
nonzero p + ^ui 2 — p' 2 > and v± > 0, v-i < 0. This allows one to transform the analytically continued 
(by means of Eq. (6.46)) integral 
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from Eq. (6.42) in the following manner (a < 1) 
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The first Appell function in (B.2) can be reduced to the hypergeometric function using the identity 
[251] which is valid for 7 = (3 + (3' 
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Thus one gets 
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This finishes the derivation of Eq. (7.2). 
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